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Abstract 



In this thesis three separate problems relevant to general relativity are consid- 
ered. Methods for algorithmically producing all the solutions of isotropic fluid 
spheres have been developed over the last five years. A different and some- 
what simpler algorithm is discussed here, as well as algorithms for anisotropic 
fluid spheres. The second and third problems are somewhat more specula- 
tive in nature and address the nature of black hole entropy. Specifically, the 
second problem looks at the genericity of the so-called quasinormal mode 
conjecture introduced by Hod, while the third problem looks at the near- 
horizon structure of a black hole in hope of gaining an understanding of why 
so many different approaches yield the same entropy. A method of find- 
ing the asymptotic QNM structure is found based on the Born series, and 
serious problems for the QNM conjecture are discussed. The work in this 
thesis does not completely discount the possibility that the QNM conjecture 
is true. New results released weeks before this thesis was finished showed 
that the QNM conjecture was flawed. Finally, the near-horizon structure of 
a black hole is found to be very restricted, adding credence to the ideas put 
forward by Carlip and Solodukhin that the black hole entropy is related to 
an inherited symmetry from the classical theory. 
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Introduction 



This thesis looks at three problems in general relativity: 

• The construction of "physically reasonable" static spherically symmet- 
ric solutions of Einstein's field equations. 

• The determination of the asymptotic spectrum of quasinormal modes 
of dirty black holes. 

• Investigation of a conformal symmetry at the event horizon of arbitrary 
stationary black holes. 

The first problem is the odd man out. It looks at the bounds that general rel- 
ativity alone can place on a star in equilibrium without imposing an equation 
of state. It is shown that no statement about a maximum mass can be made, 
although various statements about how compact an object is, can be made. 
It is also shown how to take the gravity profile of an isotropic fiuid (that is, 
one where the tangential and radial pressures are equal) and from this deduce 
the pressure profile, the density profile and (hence trivially) the mass pro- 
file. Bounds are also discussed that ensure physically reasonable anisotropic 
fluids. Three separate algorithms are presented for the anisotropic case, but 
all require the specification of a measure of the anisotropy - a hard thing to 
determine either theoretically or observationally! 

The other two problems are both "windows to quantum gravity" ; they look 
at the properties of classical black holes in general relativity and see what 
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this can tell us about quantum black holes and ultimately quantum gravity. 
There are now many approaches on how to quantise gravity (particularly 
"toy models" in 2+1 dimensions) and they all seem to reproduce the semi- 
classical formulae for black hole entropy and temperature. We are left with 
three possibilities: 

• Black hole entropy and temperature are not generic. It is a coincidence 
that every proposed model so far of quantum gravity reproduces these 
results. 

• Black hole entropy and temperature are not generic. Loop quantum 
gravity, strings and other models of quantum gravity are somehow the 
"same" theory described differently (M-theory). 

• Black hole entropy and temperature are generic. Something that all 
the approaches take for granted (such as diffeomorphism invariance) 
ultimately requires these results. As a result calculations of a particular 
model will give these results and hence cannot be used to discriminate 
between theories. 

The view taken in this thesis is that the semi-classical results are generic. 
To state this precisely one needs to be a lot more careful - extending New- 
tonian gravity to the quantum realm certainly would not reproduce black 
hole horizons, let alone black hole entropy or temperature. A less extreme 
example would be to investigate the statistics of horizon states as carried 
out by Polychronakos [^2] where the entropy was found to be proportional to 
the area provided that the individual states where somewhat distinguishable. 
Nevertheless it is believed that under suitable restrictions that the black hole 
entropy and temperature are generic in the sense that a large class of wildly 
different theories of quantum gravity will give results compatible with the 
semi-classical calculation. 

The quasinormal mode conjecture, or Hod conjecture, speculated that we 
could find information out about details of quantum gravity by looking at 
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Fi gure 1: The area of a black hole horizon is calculated in LQG by the spin links 
intersecting the horizon, indicated by black dots. On the right a link breaks off the 
horizon, and goes off toward infinity. Hod's conjecture is that highly damped modes have 
energies ~ fku, and hence relate this loss of energy to the area of horizon lost 

the resonances of a black hole relaxing back into a stationary state after a 
perturbation. These resonances, or quasinormal modes (QNMs), dictate the 
late time response of a black hole to an external perturbation. The QNM 
are harmonic functions with complex frequency, so as well as oscillating they 
are exponentially damped. This damping occurs because energy is lost both 
toward infinity and over the black hole horizon. The QNM conjecture related 
properties of these purely classical QNM to an ambiguity in loop quantum 
gravity (LQG). 

To understand how QNM may be related to LQG, an understanding of how 
LQG calculates the area of a black hole is needed. We construct a trivalent 
network, called a spin network, that underlies the geometry of space. Each 
link on the network is labelled by a representation of the gauge group of 
LQG, usually taken to be SU{2). In the case of SU{2) the representations 
are the same as the representations of particles of spin j, hence the name 
spin network. The area of a black hole horizon is calculated as a function 
of the links that intersect the horizon (see the left hand side of figure Q). If 
a particle is emitted from the horizon, as on the right hand side of figure ^ 
then the area decreases. But a decrease in area implies a decrease in mass. 
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the exact amount of mass decrease depends on the angular momentum of 
the hnk carried away and the surface gravity of the original black hole. The 
mass m, or energy, that the spin link carries off gives rise to a frequency 
uj = mc^/h. As the spin links can only by assigned half integer spins j, the 
possible values for uj for a Schwarzschild black hole are restricted to 



On the other hand, the QNMs of a Schwarzschild black hole appeared from 
numerical calculations to have the asymptotic form 



If one assumed that Hawking radiation was due to primarily spin 1 links, 
then the (real) frequency of a QNM would act like the emission of a parti- 
cle with energy E = huj\ To get the Schwarzschild mass-area relationship 
correct required the value of the Barbero-Immirzi parameter 7, an arbitrary 
parameter in LQG, to be fixed to a specific value. The Barbero-Immirzi pa- 
rameter chooses the valid sector of LQG, and so cannot take different values 
for different problems. So it seemed that assuming that treating QNMs as 
the emission of particles allowed LQG to be specified unambiguously. An 
excited flurry of work began to verify if Q could be established analytically, 
not that it was just numerically close to the "In 3" form. It was recently 
shown that the real part of the quasinormal modes was not of the "In 3" 
form by Domagala and Lewandowski and hence the QNM conjecture 
did not hold for the Schwarzschild black hole, let alone more complicated 
examples! However, the work we did in chapter |21 and chapter El was on the 
basis of extending the QNM conjecture to other black hole spacetimes, as it 
was thought to hold at least on Schwarzschild spacetime. 




(as n 00). 



(2) 



The final problem discussed in the thesis looks at why different theories 
give the same result for the black hole entropy. It may be that the entropy 
calculation only checks that the theory of gravity can be expressed as a 2D 
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conformal field theory on the horizon. The conformal symmetry in 2D al- 
lows an almost complete determination of the density of states, and hence 
the entropy. There is speculation, initiated by Cailip and Solodukhin, that 
classical general relativity can be expressed as a conformal field theory on a 
2D section near the horizon. As a Killing horizon is a null surface, a Killing 
horizon cannot be modelled by a massive field theory, hence breaking con- 
formal invariance. The research carried out by Dr. Visser, Dr. Medved and 
myself is much more modest than constructing a full conformal field theory 
modelling classical general relativity. Instead, we focus on the near horizon 
geometry of a black hole under the restriction that each of the curvature 
invariants remains finite. 

Structure of the thesis 

When writing, choices must be made about the point at which you refer to pa- 
pers for background and what you include. This thesis has been written with 
the goal of being accessible to people with a basic background in quantum 
field theory and general relativity, and much of the background material that 
could have been referenced has been included to enable this thesis to show 
why the topics included are of general interest. Although complete separa- 
tion between research and literature review is not possible, I have attempted 
to use even chapters as an introduction and motivation while expanding on 
original research in odd chapters. The research chapters give more explicit 
calculations and details to the papers included in the appendices. 

This thesis is made up of five chapters and seven appendices. In the version 
of this thesis that was handed in, the appendices are full copies of the papers 
published, and were produced as a collaboration with Dr. Matt Visser and 
Dr. Allan Joseph Medved. Because these papers are available on the arXiv, 
this is inappropriate for the thesis on the arXiv. Instead the appendices have 
been truncated to just the abstract, journal and arXiv references so that 
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people reading this can pursue the papers they are interested in. The only 
other differences between the final version and this one are this paragraph 
and parts of the conclusions section. 

Conventions 

The original source is unknown, but it has often been said that the beauty of 
standards is that there are so many to choose from. In this section I outline 
the basic conventions that I use throughout the thesis. 

Signature 

The signature of spacetime is taken to be ( — h ++) throughout. 

Differential geometry 

The Riemann tensor is defined to be 

Rabc'^OJd = (VaVb - VfcVa) UJc = [V a, V b] (3) 

for any one-form lOc- The Ricci tensor is defined to be the contraction over 
the first and third indexes: 

-Rafe = R^acb- (4) 

These conventions agree with the definitions given in chapter 3 of Wald j79j . 
The Ricci rotation coefficients are defined to be 

^iiva = e(^)^Vae(y);,, (5) 

where 6(0) are the basis vectors of a tetrad. The Greek letters are used to 
emphasise the first two indexes are labels and the last index is a coordinate. 
Technically I should write {uj^i,)a, as there are 6 one- forms ci;^,^ each with four 
components. The brackets are overly cumbersome and have been dropped. 
This differs from the convention used by Wald. 
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Similarly, the spin coefficients are given by 

i^fiua = i^u.,^be^''\- (6) 

Symmetrisation 

• The Lie bracket is defined in the standard way. For example, for vector 
fields the Lie bracket is identical to the commutator: 

[c^i^r^c'd^r-^'d^c, (7) 

• Symmetrisation and antisymmetrisation on indexes include a normali- 
sation factor: 

Aaia2...a„lbi...bm] ^ 37 signum((7) Aa^„2...a„6^(i) ...6,(^j , (8) 

III. 

■^aia2—an{bi...bm) ~ "ZT ^ ^ ^aia2—anb„(i)...b„^^y (9) 
III. 

S„i is the group of permutations on m objects, and signum(o") is 1 (—1) 
if the permutation is even (odd). Any indices that lie between vertical 
bars are exempt from (anti)symmetrisation. For example: 

Aiab) = ^{Aab + Aba) (lOa) 
Alab] = ^iAab - Aba) (lOb) 
B[a\b\c] = -^{Babc - Bcba)- (lOc) 



The nature of symmetrisation is quite different and there is no possibility 
of confusion. While it would be slightly more consistent to adopt a normal- 
isation for Lie brackets or drop normalisation for symmetrisation, a higher 
value has been placed on the ability to compare directly with the literature. 
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Veiriable definitions 



In this section I try and establish the notation that will be used throughout 
the thesis. This table is not meant to be exhaustive, but is meant to cover 
common definitions that occur throughout portions of the thesis. These are 
ordered roughly by order of appearance in the thesis, but some exceptions 
are made to keep related variables together. 
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Symbol 


C'liapters 


Moaning 


T 


All 


The stress-energy tensor. 


M 


All 


The total mass of a star or black hole, as measured from 
infinity. 


m(r) 


1-2 


The total mass which is less than a coordinate distance r 
from the origin. Also known as the Hawking quasi-local 
mass. Note the consistency condition M — m{oo). 


R 


All 


The coordinate distance that marks the edge of a star 
or the event horizon of a black hole or the Ricci scalar. 


r 


All 


A radial coordinate. Note that this coordinate is not 
the radial distance, except in flat space. If there is an 
origin, r is chosen so that the area of a set of points 
equidistant from the origin on a constant time surface 
has area 47rr^. 


P 


All 


The energy density as measured by a particular ob- 
server. 


a 


1-2 


The specific angular momentum a — J/M. 


P 


1-2 


The pressure (isotropic spacetimes only). 


Pr 


1-2 


The radial pressure (spherically symmetric spacetimes 
only). 


Pt 


1-2 


The tangential pressure (spherically symmetric space- 
times only). 


A 


All 


The area of a black hole horizon or a general area. 




2-3 


The eigenvalues of the area operator. The n acts as 
an index for different eigenvalues, and depending on the 
theory may mean a collection of indices. 


dn 


2-3 


The contribution to the area by a single spin link. 




2-3 


The (complex) angular frequency of the quasi-normal 
mode. 


N 


4-5 


The lapse function. 


n 


4-5 


The direction that is normal (and not tangent) to the 
horizon, n = on the horizon. 




4-5 


A direction that is axisymmetric and is a free coordinate 
on the horizon. 


z 


4-5 


The other coordinate that becomes a free coordinate on 
the horizon. No particular symmetry. 
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Energy conditions 

The energy conditions adopted in this thesis are summarised in tabled and 
discussed in more detail below: 

• The null energy condition (NEC): 
For all future pointing null vectors /c" 

Tabk'^k' > 0. (11) 

In terms of the principal pressures this implies that 

P + Pi>0, i = 1,2,3. (12) 

• The weak energy condition (WEC): 

Traditionally stated, the WEC states that for all timelike vectors t"' 
that 

Tabf't' > 0. (13) 

Any (unit) timelike vector can be a tangent to an observer's world 
line. The WEC condition states that the energy density measured by 
any observer is non-negative. Note that we do not need to assume the 
existence of principal pressures to make this statement. 
The WEC =^ NEC as Tab is a multilinear map and hence continuous in 
both arguments. As any future directed timelike vector satisfies fll3p 
any null vector will as well.^ 

• The dominant energy condition (DEC): The DEC assumes that the 
WEC holds, and that for all future directed timelike vectors t°- that 
—Taht'^ is a future directed non-spacelike vector. This ensures that the 

^This argument can be made because we are comparing all vectors in a particular 
tangent space, and are then unaffected by singularities in other points in spacetimc. Con- 
tinuity arguments would not hold when studying global structure, such as extending the 
chronological future of a point to the causal future. 
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Energy condition 


Definition 


Meaning (for any observer) 


Null Energy Condition (NEC) 


Tabk'^k^ > 





Weak Energy Condition (WEC) 


Tabf't'' > 


Non-negative energy densities. 


Strong Energy Condition (SEC) 


Tabf't'' > -T/2' 




Dominant Energy (DEC) 


1. WEC and 


Energy densities non-negative 


Condition 


2. -T'^bt'' is a future di- 


and the net flow of energy is less 




rected non-spacelike vector 


than the speed of light. 



Table 1: A summary of the energy conditions. In the above k"" is any null 
vector and is any timelike or null vector. The last column gives a simple 
physical interpretation for generic matter where applicable. 

net energy flow does not exceed the speed of light. By construction 
DEC =^ WEC. In terms of the principal pressures: 

p>|p,|, (z = l,2,3) (14) 

• The strong energy condition (SEC): The strong energy condition states 
that for all timelike vectors that 

Tabft' > -^T, (15) 

which by the same continuity argument used in the WEC case. In 
terms of principal pressures we have 

3 

P + J2pi>0, P+Pi>0, (^= 1,2,3). (16) 

1=1 

Note that the SEC implies the NEC, it does not imply the WEC. For 
example, matter with a negative energy density but sufficiently high 
pressures could satisfy the SEC but would violate the WEC. 

Hawking' s area theorem for black hole horizon relies on the NEC, and 
hence evaporation of a black hole must violate the NEC. Violating the NEC 
implies violating the DEC, SEC and WEC as well [THj. In fact, it is well 
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known that quantum fields violate many of the traditional energy conditions. 
In chapter ^ it is assumed that for bulk classical matter that violations of 
the energy conditions are insignificant. For a contrary point of view, see the 
essay by Barcelo and Visser [?]. Subsequent chapters look at near- horizon 
structure and Hawking radiation, so the energy conditions are not applicable. 

While the energy conditions are used to bound the types of matter we 
have available, the theorems in general relativity (such as the singularity 
theorems) usually rely on statements about the geometry. By assuming the 
Einstein equations the energy conditions can be transformed into statements 
about the geometries we are interested in. The NEC becomes the condition 
the Rabk"'k^ = Gabk"'k^ > 0, the null Ricci convergence condition. The WEC 
becomes the Einstein convergence condition Gabt""^^ and the SEC becomes 
the Ricci convergence condition Rabt°'t'^. 

The definition of the energy conditions used conforms with Hawking and 

Ellis m]. 



Chapter 1 



Constructing perfect fluid 
spheres 



1.1 Introduction 



The most basic model for a star is one in hydrostatic equilibrium, a condition 
which connects the rate of change of pressure throughout the star with the 
density profile. When supplemented with a relation between p and P (a 
so-called equation of state) we can then integrate out the equations and find 
both the density and pressure profile. The simplest example of this method is 
to assume a fluid with P = p'^, an example which is treated in most textbooks 
|32l EOj- The disadvantage of such a method is that we need to know not 
only general relativity but also a great deal of nuclear physics at extremely 
high densities to get realistic profiles out for compact objects such as neutron 
stars. If more than one type of matter is present or different phases exist 
then the "equation of state" used to integrate the density profile is really an 
equation that ties the density at a specific radius to the pressure. Hence the 
function P{p) may not be invertible. Even in the cases that such a relation is 
invertible it will generally be very different from the equations of state of the 
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matter making up the star.^ Beyond the purely technical problems of finding 
equilibrium pressure and density relations there is also the problem that we 
do not have direct experimental evidence for how matter behaves under such 
pressures. This is refiected in the difficulty of finding a strict mass limit for 
neutron stars and other superdense objects such as quark stars [32] . 

An alternative algorithmic approach can be used to generate any spheri- 
cally symmetric static solution of Einstein's equation without recourse to the 
details of the material that makes it up. For Newtonian gravity the condition 
of hydrostatic equilibrium is quite simple: 



where g{r) is the local acceleration due to gravity, m(r) is the mass contained 
within a radius r and p(r) is the density. The mass within a radius r is simply 
related to the density: 



^0 

If we tried to use p(r) as a known function we have a second order differential 
equation to solve: 



We are relatively fortunate in this the differential equation has a nice 

solution, albeit in terms of an integral of the density: 



"'^An extreme example is a constant density star. By combining different phases of 
matter it is possible to have a fluid of constant density, but it would still be compressible 
(i.e. at other densities the relation p = Const would not be valid), whereas if the relation 
p = Const was assumed to hold true at all pressures then we would have a incompressible 
star and infinite speed of sound. |57| 




(1.1) 






(1.2) 
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This is what is meant by an algorithmic procedure: given an arbitrary func- 
tion (in this case p(r)) we can then use that to generate an exact solution 
for the pressure profile ()1.2j) and mass distribution m(r) (trivial). Because 
the pressure does not act as a source in Newtonian gravity we have no way 
of knowing what the tangential pressure is doing without some assumption 
of the nature of matter making up the star. 

The above algorithm is effective and simple and it will allow us to find 
the (radial) pressure profile for any time-independent spherically symmetric 
Newtonian fluid with no external forces. By allowing arbitrary functions for 
p{r) we could have solutions that are clearly unphysical, such as any solution 
with p(r) negative. The advantage of this particular algorithm is that ensur- 
ing that p(r) is positive is trivial, as this function is freely specifiable. The 
only non-trivial constraint that we can sensibly impose on the Newtonian 
case is that the pressure remains finite everywhere. That is, we demand: 

• The exact solution can be integrated out for radial pressure. 

• Any reasonable conditions {p > 0, Pr finite) are easily imposed. 

• The "generating function" is not an esoteric function but has direct 
physical relevance. 

Unfortunately attempts to generalise this equation to relativistic stars fail as 
the generalisation of ()1.2|) is a non-linear second order differential equation. 
The non-linearity arises for a simple physical reason; now the pressure also 
acts as a source of gravity as well. 

In contrast, let us return to and now let g{r) be the freely specifiable 
function. Simple algebra allows us to find m{r), and differentiation allows us 
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to find p(r) and then Pr{r): 

m{r) = r'^g{r) = An j r^p{f) df 

^ p(r) 







dr 



47rr 



2 
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The differential equation has roughly the same form for relativistic stars, 
but now pressure and density are sources. Despite being non-linear, it is still 
a first order differential equation and specifying the function g{r) allows one 
to construct the density and pressure profiles. In both the Newtonian and 
relativistic cases we would like to ensure that: 

• The energy density is non-negative, and 

• The pressure is finite and non-negative. 

It is comparatively more difficult to ensure the non-negativity of energy and 
pressure in the second algorithm, as it involves choosing a function g and 
then computing the density. In the relativistic case we also want to ensure: 

• The compactness x(^) = 2m(r)/r is always less than one (no black 
holes), and 

• That the various energy conditions are satisfied. 

These conditions have no sensible analogue for perfect fluid spheres in New- 
tonian gravity. 

The rest of this chapter is devoted to properties of perfect fluid spheres in 
general relativity. First an algorithm for perfect fluid spheres analogous to 
the one presented here for Newtonian stars is developed. It is then interesting 
to ask questions like 

• Can we place bounds on the generating functions g{r) so that we get 
reasonable solutions? 
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• Can we bound the internal compactness 

• Can we find more general bounds than the Buchdahl-Bondi bound 

X = X{R) < 8/9? 

A few comments are in order on why we have restricted ourselves to looking 
at time independent non-rotating fluids, as most stars are radiating (destroy- 
ing the time independent condition^) and rotating (which leads to non-static 
metrics and equatorial bulges breaking spherical symmetry). The effect of 
radiation can be treated as a perturbation in almost all of the astrophysically 
relevant examples, and on small timescales is swamped by the effect of rota- 
tion. The problems facing an algorithmic procedure for a stationary rotating 
star or planet are pragmatic ones. The first is that there is no equivalent of 
Birkhoff's theorem for rotating systems, so we have no reason for thinking 
that the exterior solutions to rotating bodies are described by the Kerr met- 
ric [Zni particularly §7.1, §12. 3]. That does not stop us from restricting our 
attention to only those stationary solutions that do join onto the Kerr met- 
ric as a first attack on the problem. However, despite numerous attempts no 
solutions of the Einstein equation in reasonable matter have been found that 
join smoothly onto the Kerr vacuum, although various perturbation tech- 
niques have been developed [2011221 • Thus the idea of being able to generate 
all physical solutions to describe rotating bodies "algorithmically" would be 
far-fetched at the current time, as not even a single explicit solution is known. 
To give a rough estimate of the order of the effects being neglected, a table 
of values has been constructed (table ITTT|) . In a asymptotically fiat, spheri- 
cally symmetric system Birkhoff's theorem is the statement that the solution 
in vacuum is a piece of the Schwarzschild solution (see §1.2.11 for a proof). 
The number of degrees of freedom inside matter is also tightly constrained, 
making stationary spherically symmetric metrics suitable for solving algorith- 
mically. We shall see that the algorithm becomes complicated even under 

^For all reasonable systems involving finite energy. 
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these restrictions, and that a useful algorithm only exists for matter models 
that are also isotropic. 



Name 


Mass 


Mean radius 


Ai?/i?eq 


M/R 




a 


/R 




a/M 


Earth 


4 mm 


6.4 X 10^ km 


3 X 10^^ 


7 X 10-1° 


3 


X 


10- 


-7 


429 


Jupiter 


1.4 m 


71500 km 


7 X 10-2 


2 X 10-* 


1 


X 


10- 


-5 


524 


Sun 


1.5 km 


695 500 km 


9.1 X 10^*^ 


2 X 10-6 


5 


X 


10- 


-7 


0.234 


Neutron star 




M~ 


1.4M0, i?- 


20 km, period 




1 


s 








2.1 km 


20 km 


< 1 


0.105 


2 


X 


10- 


-4 


1.9 X 10-3 



Table 1.1: A table showing the typical size of rotation effects on astrophysical bodies 
in units in which c = G — 1. Notice that for bodies with low compactness {2M/R) it is 
possible for the parameter a/M ^ 1. The oblateness, defined as the difference between the 
polar and equatorial radius over the equatorial radius, is also given to indicate departure 
from spherical symmetry. 



As seen in table the specific angular momentum may exceed the mass 
of the star or planet, unlike a black hole. An extreme example of such a 
violation comes about from considering a planet orbiting a star. The (orbital) 
angular momentum can be approximated by Newtonian gravity: 

• - -ip^/M^, (1.3) 



nipr u = TTir, 



where rup is the mass of the planet and Mq is the mass of the star, and 
rrip <^ Mq. The specific angular momentum a is then 

L m„ t— — m^) 



M^^mp MQ^mp^ Mq ^ " ^ ' 

Expressing this in terms of the Schwarzschild radius of the planet and star 
gives 

Rsp ^ jRsQr ^^^^^ 
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RsQ V 2 

or that the ratio of specific angular momentum to mass is given by 
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The orbital radius r is of the order of astronomical units, making this ratio 
incredibly large. Taking Jupiter and the Sun gives 

(^) "21. (1.7) 

\M / Jupiter + Sun 

Looking at the ratio a/M for planets alone typically yields enormous values. 



1.2 Setting up the problem 
1.2.1 Spherical symmetry 

We start with a spherically symmetric metric and try and see what con- 
straints this places on the geometry. The presentation given here closely 
follows Misner, Thorne and Wheeler jSTj (§23.3). By spherical symmetry it 
is meant that the group of isometrics of spacetime possess SO (3) as a sub- 
group. It is also assumed that for every point x except the origin that the 
action of this subgroup is a two-dimensional surface, denoted S{x). Natu- 
rally the origin, if one exists, will be mapped to a zero dimensional surface 
(namely the origin itself for any rotation). For nonstandard topologies (such 
as wormholes) it is not necessary for an origin to exist. Each of these two 
dimensional surfaces can be given an induced metric of a sphere 

ds\^^)=R{xfdn\ (1.8) 

where df2^ = dO"^ + sin^ Odcj)^ is the metric on a (unit) two-sphere and R{x) 
is a "scaling" function that depends on the point x of the r — t submanifold. 
The spheres are a submanifold and are also maximally symmetric implying 
that we can always separate the coordinates on a submanifold 

ds^ = gQQ{r,t)dt^ + 2gQi{r,t)drdt + gii{r,t)dr'^ + R{r,tfdVt^ (1.9) 

= ^00 (dt + ^dr) ' + - f \ dr' + i?(r, tf dU' (1.10) 
V ^00 / V V^oo// 
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In Misner et al j^Z] they carry out this construction in considerably more 
detail but at the cost of assuming that the spacetime is stationary. The time 
coordinate can be redefined as i in the following way to ensure the cross-term 
vanishes. Introduce the (as yet undetermined) function f{i,r): 

t = t + f(t,r) (1.11) 
We then ensure that the cross term vanishes by defining fit, r) as 

We have now got a definition of t in terms of t, up to an arbitrary constant. 
The metric now takes the desired form: 

ds^ = goodt" +(gn- —] dr^ + R(r, tf dfil (1.14) 
V 9oqJ 

Finally the radial coordinate is redefined so that R{x) = f. This means we 
can now write the metric in the form (dropping the bars): 

ds^ = -C(r, t)^ dt^ + 7]{r, tf dr^ + dfi^ (1.15) 

for some functions rj and C^. The squares of the functions have been introduced 
to ensure that spacetime has a Lorentzian structure. 

Mathematically there are lots of things that can go wrong with this con- 
struction globally. The "timeshift" function f{i,r) or the "radial function" 
R{x) may not be invertible, but the implicit function theorem (almost) guar- 
antees that locally we can put the metric in this form. Notice that the only 
condition so far is spherical symmetry, the field equations have not yet been 
imposed. Doing so in vacuum leads us to Birkhoff's theorem: 
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Theorem 1. (Birkhoff's theorem) In a spherically symmetric spacetime the 
only solution to the Einstein equations (without cosmological constant) in 
vacuum is a piece of the Schwarzschild solution. 

Proof. The field equations state that the Einstein tensor vanishes identically 
in vacuum. Calculating the Einstein tensor from the metric in an orthonor- 
mal frame gives the following equations: 



W -C-2rg 



Gff = " 2 = (1-17) 



Gi, = = (1.18) 

Gee = (1-19) 
% = Geh (1-20) 

where the last expression follows from spherical symmetry, not from the 
requirement that the Einstein tensor vanishes. The angular expressions have 
not (yet) been given as they are incredibly ugly in their current form; huge 
simplifications are obtained by solving the other equations first. As a point 
of notation the Einstein tensor has been evaluated in an orthonormal frame, 
which will be denoted here and throughout by placing "hats" on the indices. 

The Gff equation tells us that the function rj is independent of time. The 
Gff equation can then be explicitly solved for 77 (r): 



^ r-2M 

where M is currently just a (fortuitously named) constant of integration. 
Substituting 77 (r) back into the G^f equation gives 

2M^ 



ar,ty^ary^A\^i 

where A is another constant of integration. 
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The next logical step is introduce the angular components of the Einstein 
tensor and check that they do indeed vanish: 



dr^ dr ' \ dr dr dr^ ' 



or ^ or I I 

= 0; 

as it had to be, by virtue of the Bianchi identity. By rescaling the time 
coordinate the metric can be put into standard Schwarzschild form: 

ds^^-(l-—] dt' + —^dr' + r'dn'. (1.21) 



□ 



It should be noted that before applying the Einstein equation we were 
able to transform into coordinates that ensured the metric decomposed into 
spacelike and timelike parts i.e. gu — 0, i — r,9, (f). Hence if the spacetime is 
stationary (i.e. there is a timelike Killing vector) and spherically symmetric 
automatically imply that it is static (i.e. that there is a spacelike hypersurface 
orthogonal to the orbits of the timelike Killing vector field) . An example of 
a stationary but non-static case is the Kerr black hole, as the rotation rate is 
not changing but the surfaces orthogonal to the timehke Killing field are not 
spacehke. In looking at stars we shall always be looking at stationary cases 
(and hence static by spherical symmetry) so no use will be made of this fact 
in this thesis. 



The final comment about Birkhoff 's theorem is that we know that spher- 
ically symmetric spacetimes do not radiate - even if the entire spacetime is 
not static! More specifically the vacuum region is forced to be static even if 
the matter regions are undergoing (spherically symmetric) pulsations. 
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1.2.2 Spherically symmetric matter 

So far we have only looked at the field equations in vacuum. The full Einstein 
field equations are 

G^, = SvrT^, (1.22) 

where T^j^ is the energy-momentum tensor. Performing the necessary incan- 
tations of the previous section we take the metric to be of the form 

ds' = -C(r)2dt2 + ^— -dr^ + rMfi^ (1.23) 

1 — 2m(r) jr 

I am putting in an additional assumption now that the spacetime is station- 
ary, which is not guaranteed inside matter. I have also replaced ri{r) with 
the function m{r), a label that will be justified as we look at the field equa- 
tions. Calculating the Einstein tensor directly from the metric ()1.23|) in an 
orthonormal frame to find (where primes denote differentiation with respect 
to r): 

2 d?Tz 

G,,- = 87rp=-— (1.24a) 
G,, = SttP, = 4 fl - — 1 C' - ^ (l-24b) 




G,, = 8nP, = -Wl - — I C'a/1 - — 1 + ^ ( 1 - — 1 - - (-)' 

r \ r / 

(1.24c) 

% = G^§, (1.24d) 

where we have used the Einstein field equations only to relate G to the prin- 
cipal pressures. Again we see that (jl.24djl holds and is a direct consequence 
of spherical symmetry. The label m(r) is justified as the "mass inside radius 
r" by integrating from an arbitrary but constant Tq to p.24a|) 

m{r) = 4:71 / f^p(f)df + mo, 



■^Thc mass "inside" a radius r only makes sense if we can define a point r = 0. In cases 
that we cannot define r = (wormhole throats, for example) it would be better to say the 
mass between vq and r). 
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which is the Hawking quasi-local mass in spherical symmetry. This is pre- 
cisely the equation you would expect for a spherically symmetric mass in a 
Euclidean space. In a curved space one would anticipate that the volume 



the proper mass jZHI, but as gravity itself gravitates we also have a "grav- 
itational mass". In spherically symmetric spacetime the gravitational mass 
exactly compensates for the smaller measure being used leading to the naive 
formula one would expect in flat space. This result is not generally true, and 
unless the spacetime is asymptotically simple there does not appear to be a 
meaningful way of defining a mass. The asymptotically flat case is discussed 
in chapter 11 of Wald's book [79J. Note that if we require the geometry to 
be regular at r = (which includes the case where we do not want a black 
hole) then mo = 0. For future reference the local acceleration is given: 

Theorem 2. The local acceleration measured by a (timelike) observer in a 
static spherically symmetric spacetime is purely radial with magnitude a{r) = 



Proof. We use the fact that the spacetime is static, so that we have a Killing 
vector field ka = d/dt. Let us have an observer who is "at rest", meaning 
that the four velocity is parallel to the Killing field (i.e. (r, 6, (p) is constant). 
We choose to normalising the four velocity by 



form would have to be replaced 




This would define 



y/1 - 2m{r)/r C'(r) /C(r) . 






which is possible as the motion of the observer is timelike. These conditions 
imply that V = (C"S 0, 0, 0) = C"^ k". 
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We can then calculate the four acceleration as follows: 



a« = V^VbV (1.26) 
= h^Vb (hA (1.27) 



C "VC 

= + ^^'V,r . (1.29) 

In this case we know that C is a function of r alone, so VC~^ will be pointing 
in the purely radial direction. As the Killing field k"" is pointing along (and 
indeed defining!) the timelike direction, the first term must be zero. 
1 



a" = -^k^{'Vbkc)g"''^ (as we have a metric tensor) 
1 



— —^k^{yckb)g°''^ (Killing's equations) 
'l ... . 1 



2\ „ac 



,(yck'kb)g^' = -—jycC)9 



2(2^ '' "'^ 2(2 
9 ^ 

which is zero unless a = r. For the radial acceleration we have 

a"- = (1.30) 

which completes the proof. □ 

Instead of using the acceleration due to gravity directly it is often more 
convenient to express things in terms of g{r): 

Definition 1. Let g{r) be related to the gu component of the metric as fol- 
lows: 

g{r) ^ \%M-9u) (1-31) 

Cir)' 

- W) 

= grr {no sum) (1.33) 
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Alternatively we can write the gu component of the metric in terms of 
g{r) by a simple rearrangement of the definition: 

gu = -exp(^2 j g{f) df^ . (1.34) 

We can eliminate C,{r) from the Einstein tensor in favour of the g{r), a 
quantity that has a direct and simple physical interpretation. The Gfr equa- 
tion can be rearranged to give 

m{r) +A-Kr^Pr{r) 

9{r) = ^ . ■ (1.35) 

r(r — 2m[r)) 

The Gfjg equation takes a little more work. By substituting C('^)' = C(^)fi'('^) 
we find 



+ ^ 1 - 

r \ r J r \ r 

i(.'C..O.f)(i-^)-(..i)(^ 

(7)' 

The algorithm for isotropic stars is obtained by taking Pt = P^.. When 
coupled with the equation of hydrostatic equilibrium (discussed below) the 
solution can be found in terms of g{r) quite simply. 



1.2.3 Hydrostatic equilibrium 

For Newtonian self-gravitating fluids the equations of hydrostatic equilib- 
rium are obtained easily and are physically motivated (see figure IT!T|) . For 
anisotropic relativistic fluids a nicely motivated approach does not appear to 



15 



P(r)A 



pg + P(r+dr)A 

Figure 1.1: The forces acting on a fluid parcel element subject to Newtonian 
gravity. The tangential forces are ignored as they do not play a role in the 
equilibrium condition. Here P denotes radial pressure and A denotes the 
area of the fluid parcel element 




be possible. We start with the expression for g{r) p.35p and differentiate: 

m' + 47r(r3p,)' m + A-Kr^Pr 

g = — 2{r -m-mr)— — — (1.37) 

r[r — 2m) r^(r — zmj"^ 

m' (r^P.)' , , X 9 . X 

+ 47r-^ ^ - 2(r - m - mr) , \ . (1.38) 



r(r — 2m) r(r — 2m) r(r — 2m) 

1 



;i + 2^r)m' + Vl'nr'^Pr + ^■nr'^P[ 



r{r — 2m) 

-2{r-m)g}. (1.39) 



We can solve the G^g equation ()1.36p to find an expression for g' as well: 

, (SirPAr or + 1 /m\' n Q , ^ 



r — 2m r — 2m V r / r 
Eliminating g' we solve ffO^ and ([00)) for P,': 

Anr^P'. = SnPtr'^ + (1 + 9r) (^(— ) i" - "^'^ " i^'9i" (1-41) 
— g'^r(r — 2m) + gr — 12nr^Pr. 
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The derivatives of m(r) are eliminated via the equation ()1.24ap : 

47rr=^P; =Mr\Pt - Pr) - ""^^ ^ ^"^^ - Aur^gp (1.42) 

r 

— g^rir — 2m) + gr — Arrr'^Pr. 

We expand all but one of the factors of g{r) and collect terms: 

47rr^p; = - 4nr^P,. - Aur^gp + Sr^r'^^Pt - P^) (1.43) 

+ — — -| (m + Anr^Pry - (m + Anr^Prf 

r[r — 2m) y 

— m{r — m + Anr^Pr 

Some remarkable cancellations in the bracket lead to 

47rr^P; = - Anr^Pr - Anr^gp + 8nr\Pt - Pr) (1.44) 
Aixr^ P 

+ \{r - 2m) - (m + A-nr^P,)] , 

r[r — 2m) 

or, in its final form 

^._,(,,P,),H(iWl). (1.45, 
dr r 

This equation is known as the anisotropic Tolman Oppenheimer Volkoff 
(TOV) equation. In the case of an isotropic fluid we have the more familiar 
TOV equation. 

(1.46) 

1.3 Algorithmic construction 
1.3.1 Isotropic fluid 

The most difficult thing about the algorithmic constructions is deciding which 
functions we are going to eliminate and which we are going to keep. We have 
a plethora of functions to choose from : P = P^ = Pt, g, p, m and (. By 
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referring to p.lSj) . we see that spherical symmetry allows us to reduce our 
problem to at most two free functions . As we are interested in isotropic 
fluids the extra condition Pr = Pt restricts us to one arbitrary function, the 
so-called "generating function" . 

While in principle we have one freely specifiable function we are not guar- 
anteed at this stage that we will be able to find a closed form solution for the 
resulting differential equation. If we were to take the mass as a "generating 
function" then finding p is trivial: 

Then we try and construct the pressure profile: 
dP 



dr 



g{p + P) (1.48) 

m + 47rr^P\ , , , , 

iP + P) (1-49) 



r(r — 2m) 

mp f'm + ATTr^p\ Anr 



. , , sP-—, -.P ■ (1-50) 

r(r — 2m) \ r(r — 2m) / r(r — 2m) 

This is a non-linear differential equation although it will in general have a 
solution, and for some choices m(r) there will be a closed form solution, there 
is no closed form solution for arbitrary functions mir). The culprit is the 
last term which makes the equations non-linear. Note that in the Newtonian 
case the first term would not have m in the denominator and the other two 
terms would not exist at all; hence in cases where the star is not too compact 
one can neglect the last term, integrate the resulting equation and have an 
approximate solution significantly better than the Newtonian solution. This 
is not of much use to us as we are looking for exact solutions and bounding 
the compactness by the condition that the pressure remains finite; in such 
limiting cases the last term becomes dominant rather than negligible! 



Attempts to generate solutions from the density p(r^ leads to a very sim- 
ilar conclusion. To eliminate m(r) we need to solve ()1.35|) and differentiate 
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now leading to a second-order non-linear differential equation. After these 
attempts it may be a surprise that any closed form algorithm exists at all. 



To actually construct the algorithm we start with the isotropic version of 

dnni): 

9 = TT- + ( - ) -9'--- (1.51) 



r — 2m r — 2m V r / r 

We choose to write our function in terms of mass and the gravity profiles. 
We eliminate P by using the Gfr Einstein equation ()1.24b|l : 



, q 2m qr+l /mx' r, 

9 =--—f ^ + ^^ - '9'- 1.52 

r r''[r — im) r — im \ r / 



This is a Ricatti differential equation in g, so again we cannot specify m(r) 
and find a general closed form solution. However, when considered as a 
differential equation in m(r) we see that it is actually a linear differential 
equation: 



m' = ^ [g' + 9'-^)+ (3(1 + 9r) - 2r\el + /)) . (1.53) 



As a first order linear differential equation we are able to integrate it directly 
by introducing an appropriate integrating factor. 



The integrating factor (denoted by A) can be found by solving the homo- 
geneous differential equation: 

y . Af?_2.to:±i!)V (1.54) 

\r l + rg ) ^ ' 
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The solution takes the rather simple final form: 



= exp ^3 In \r\-2 J ^'[^^J dr^ (1.56) 
= r3exp(-2|M^^^±^dr) (1.57) 
= r^exp ( -21n|l +r5r| + 2 / g- — — dr] (1.58) 



1 -\- gr 

exp {2 [ g ]~^^[ dr ) , (1.59) 

ro 



[1 + rg)^ V Jro ^ + 9r 



where in the last integral we have explicitly introduced the constant of in- 
tegration tq. It is now a simple matter to obtain a closed form solution for 
m(r) (with the constant of integration absorbed into the indefinite integral): 



T'^ ( 1 — gv 
'>n(r) = — rrrexp 2 / g- dr 1 (1.60) 



(1 + rg) 

where we have introduced the function 



I ^e-^'^C^) {{gr)' - 2g + g\) dr, (1.61) 



Formally wc have done what we have set out to achieve; given a generating 
function g{r) we can construct the mass profile and from the Grf equation 
we have the pressure profile. 
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We can improve it by eliminating the function g' by an integration by 
parts. Performing integration by parts on the offending term: 

i±^(^?r)'e-2''W dr = y {gr)''-^ dr + ^ j (A')'^^ dr (1.63) 



After a httle bit of algebra we find that m(r) is given by 
(l + e^'^ / 2^^(^V - 3) 



gj-2 , ^3 

m(r) 



{l + rgy 

where there is still a constant of integration to be determined. 



For stars with finite central pressure the integrand in p.(j7|l is well be- 
haved. Sufficiently close to the centre of the star we can make the following 
approximations: 

4 

m(r) ^ -vrrVc, P{r) ~ Pc (1-68) 
3 

which implies 

g{r) ^ |vrr-^^^4^ ^ |vrr(p, + 3P,) (1.69) 
3 1 — Snr'^pc/S 3 

i.e. the gravity in well behaved stars is linear at small r, making the integrand 
well behaved. 
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Another "well behaved" algorithm comes from considering a different in- 
tegration by parts in ()l.(il|l : 



1 + rg 



{{gr)' - 2g + g^r) e'^^ dr 



(l + r<7)(^-4 + ^)e-2'^dr 

^-^ ^99'-^-^{\^rgf^ e'^^dr 
rJ r r 



r 2 



(1.70) 
(1.71) 
(1.72) 
(1.73) 
(1.74) 



+2 



1 + gr J \r 2 



dr 



9_^^2^ 

r 



(^^f)^^/f(r^){(^^f)(--) 



9^ , 



V-(l + f: )+2 



r(l + rg) 



(1.75) 
(1.76) 



which leads to the algorithm [HH] : 



( \ f. . 9r 



2^3 



9r 



{l + rgy 



1 



gr 



(l + r^)2 
(1+^ 



e-2^^dr 



r(l + rg) 



3 r{l + rg) 



(1.77) 



dr , {l.n 



where the constant of integration has been put in explicitly. This algorithm 
is slightly better than ()1.67|1 insofar as g only enters as a quadratic rather 
than a quartic. To compute the pressure we simply use the Gff equation: 



p{r) 



1 



87r[l + r g{r)Y 



-g^-2[l + 2rg]e'UCr,+ 



19' 



r(l + rg) 



-2^ dr 



(1.79) 
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1.3.2 Anisotropic fluid 

We now consider the more general case of an anisotropic fluid: the radial 
pressure Pr is not required to be the same as the tangential pressure Pt- This 
case is a little more artificial as we now have the freedom of both arbitrary 
functions of r allowed in the metric by the constraints of a static spherically 
symmetric spacetime. This tends to make the results of this section somewhat 
"artificial" and so the derivations are done in less detail. 

We are able to tell about the gross features of anisotropy by looking back 
at the anisotropic TOV equation ()1.45|) : 



The first term is non-positive by the NEC {\Pr\ < p, p > 0).^ The radial 
pressure is zero at the surface of the star by definition, and in the isotropic 
case the derivative shows that Pr comes to zero from above so Pr > 0. The 
more negative dPr/dr is the larger the central pressure, and hence we have: 

• Pt > Pr- One would expect this condition to hold in rigid solids. The 
tangential pressure would help the solid retain its shape. Thus less 
radial pressure would be required to stop collapse. The fall off of Pr 
as the surface is approached will be gentler, or equivalently the central 
pressure will be lower, than an isotropic star of the same radius and 
mass. We can see this mathematically by looking at ()1.45|1 and seeing 
that the second term is positive. 

One may have Pt slightly larger than Pr in "rocky" planets such as 
the Earth. Astrophysically compact objects such as neutron stars and 
white dwarfs may form crusts that would allow the tangential pressure 

^We cannot assume that Pr > 0. For example, we can have an empty soUd sphere 
immersed m an atmosphere. The radial pressure is pushing inward {Pr < 0), but the 
rigidity of the surface creates a tangential pressure which can support the object against 
collapse. 




(1.80) 
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to exceed the radial pressure. The most extreme example of this effect 
would be a Dyson sphere, which relies on the tangential pressure over 
the surface of the sphere for its stability. As the central pressure is 
lowered these stars are able to support a higher maximum mass than 
their isotropic counterparts. 

• < Pt < Pr'. This condition would occur in fluids with weak surface 
tension; gravity is still needed to hold the fluid together. Most "normal" 
stars would satisfy this condition, although the effect will typically be 
negligible especially compared to rotation. As a consequence of the 
surface tension the radial pressure profile must be steeper and hence 
the central pressure is greater than its isotropic counterpart. 

• Pt < 0: This is the case for fluids with strong surface tension, or so- 
called self-bound matter (see chapters 8 and 12 of |S2])- These are forms 
of matter that do not require the gravitational force to be bound, such 
as small water droplets. The radial pressure in such cases does not 
have to go smoothly to zero at the surface. The radial pressure profile 
will be pushed up further in comparison to the isotropic case, so the 
central pressure will also be higher. 

We have experimental evidence of objects satisfying this condition but 
they are all small such as raindrops. There is some speculation that 
hybrid or quark stars may satisfy this condition; this is an assertion 
that the nuclear force is capable of binding the star even if gravitational 
forces were absent. However there is no observational evidence of any 
astrophysically relevant body satisfying this condition. 

These cases show the qualitative difference that introducing anisotropy makes. 
One needs to be careful as it is possible to construct solutions where neither 
the radial pressure or the tangential pressure is dominant throughout the 
star. The requirement of spherical symmetry requires that Pr = Pt = Pc at 
the centre of the star as well, so in particular we cannot have a static solution 
where Pt < all the way through the star. 
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Although the effect of anisotropy will typically be much less than that of 
rotation in stars, for completeness the problem is addressed. Three separate 
algorithms are discussed: 

• By using g{r) and the anisotropy function {Pj. — Pt) to generate the 
mass and pressure profiles. 

• By using g{r) and the tangential pressure to generate new solutions. 

• By using g{r) and 5{r) = Pt/ Pt- 



1.3.3 The anisotropy function solution 

This solution is perhaps the most "direct" attack as we simply integrate the 
first order linear differential equation. With the integrating factors somewhat 
simplified we obtain 

^{l + \rg{r)) 



Mr) = T^^rTTTM^e 



2i?(r) 



+ < q{r)r , , 
exp(27?(r))r^ / /" g{r)'^ 



+ ,r /\ \9 C + 2 / ,, exp(-2^(r)) 
= m^mso + friiso- (1-82) 

In the last line we have written the mass as the sum of an isotropic mass and 
an anisotropic mass. The isotropic mass is defined to be the mass of a star 
with generating function g{r) but with P^ = Pt throughout. The anisotropic 
mass is defined as whatever is left over. By comparison with our algorithm 
for isotropic fluids p.78|) we see that the first term in ()1.81|1 is the anisotropic 
mass and the second term is the isotropic mass. 



An alternative definition of an isotropic mass would be the mass profile that 
generated Pr, under the constraint that Pt = Pr- While this is a legitimate 
construction [although not able to be solved in closed form in general] it is 
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distinct from the one I have chosen to use. It is important to note that the 
pressure profile corresponding to the isotropic mass profile the way I have 
defined it will not give the same Pr- 



Some general comments about this solution: 

• It is very easy to find the isotropic limit as this function decomposes 
into an isotropic and anisotropic part. 

• Care must be used when saying the "isotropic part" . What I am assert- 
ing is that this is the mass generated by an isotropic fluid sphere by the 
function g{r). I am not asserting that P^, (isotropic) = P,. (anisotropic). 

• The function (P^ — Pt) is a difficult thing to determine experimentally. 
Theoretically the best limit I can place on it is 2p by the DEC. 

• As our boundary condition we can choose lim^^o = f^rpc- We also 
choose the integrals to have lower limits of integration at zero, allowing 
evaluation of the constant of integration: 

m(r) 4 , q(r) ^ 

lim^-^ = -7rpe = lim^ + a 1.83 

With the expansion of g for small r given by we find that C = 

—4:7Tpc- The advantage of this boundary condition is we do not need 
to know anything about the anisotropy to calculate it. (For a look at 
other types of boundary conditions on isotropic stars see appendix lB|) . 

• The posit ivity of the mass and the isotropic mass allows us to (trivially) 
state 

-"^aniso < "^iso (1-84) 

for realistic stellar models. Note that this is not guaranteed by the algo- 
rithm, only by our prejudice for what constitutes physically reasonable 
solutions. 
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1.3.4 The tangential pressure solution 

Using the Gff equation in terms of g we have 

^-KPy 2rg{r) 2m{r) fl + 2g{r)r 



1 + rg{r) 1 + rg{r) r V -'- + / 
Using this we now have 



:i.85) 



m{r) '^'^'' + <!'^''^^:!"'^''^' /(r) dr + (1.86) 

/(r) Vi r\l + rg{r)) 2j y j 

where the integrating factor /(r) is given by the rather unsavoury expression 

7(r) = exp f 2 / ^rgjr) + 1 + r^r)^ + g'jry \ 
^ ^ \ J r{l + g{r)r) J ^ ' 

An integration by parts simphfies the integrating factor considerably: 

7(r) = r^(l + ry(r))^exp ^2 j dr) . (1.88) 

Yet another integration by parts will get rid of the derivative of g{r) in 
the main expression, as well as some algebra and actually doing one of the 
integrals (!) we finally obtain: 



r 



g{rYr 



mir) = - + - ---rrexp -2 / ' dr\ (1.89) 

^ ' 2 (l + r^(r))2 ^ V J ^ + rg{r) ) ^ ' 

X jcs -^T^ j Pt{l + 5(r)r)rexp (2 J ^^^^^ dr 

• This solution is fairly "tidy" 

• Finding isotropic solutions is messy (eg. substitute for Pt, use the Grr 
equation to find Pr and then use this for new Pt- Rinse, lather, repeat.) 

• Is Pt easily measurable? Can it be estimated reahstically? This is a 
difficult question, although I suspect that the answer is no. 

• To guarantee that 2m{r)/r < 1 the term in curly braces must be neg- 
ative. 



27 



• Boundary condition at r = 0: m(r) ~ ^pj'^ implies that C3 = — |. 
(All top level integrals run from zero to r, nested run from zero to 
integration variable). 

• Setting the lower bound of integration of the integrating factor at in- 
finity and using the boundary condition that m{r > R) — M and 
P{r > R) — gives 

M=- + CA -T^exp 2 / _._ — df (1.90) 



2 " [r-MY \ Jr f{f-2M){f-M) 
r r{r-2Mf ( {r-MY \ 

2^ ' {r-Mf \r{r-2M)) ^ > 

= ^ + C3(r-2M) (1.92) 

which gives the simple condition C3 = —\. 

• An easy inequality comes from requiring that 2m{r)/r < 1 throughout 
the fluid: 

Cs-SttJ Pt{l + g{r)r)rexp(2 J -^i^^^dr^^ < 0, (1.93) 
It is less useful than one might hope as it involves an integration. 



1.3.5 The 6 solution 

Consider the equation 

^ ^ 2g{r)r{l-5) {g{r)^r + g'{r)r - g{r))r 

^ ' l + g{r)r 1 + g{r)r ^ ' ' 

\r l + g{r)r r 1 + g{r)r 

(1 - 2S)g{r) + gjrfr + g'{r)r ^ 
1 + g{r)r 

'3 2 fg{rfr^ + {g{ryy + {1 - 6) - 26g{r)r 



r r \ 1 + g{r)r 



with S{r) = Pt/Pr 
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The integrating factor can be written down in closed form, but is long. 
The function / is introduced to shorten the notation 

g{r)r + 1 - 26{r) 1 - 5{r) 



exp 2 / /(r) dr = exp 2 



9{ry- 



1 + g{r)r r{l + g{r)r) 



dr 
(1.96) 



The solution for m(r) is given by 

g(r)r^(l + r^(r)/2) rV^J/^'- f f g irfil + 6 jr)) ^Jf^r,^ 
^ ' (l + ^(r)r)2 ^(i + ^(r)r)2 Vi r(l + ^(r)r) ^ 

(1.97) 

This is an algorithm, although not easy to use or find properties of. In 
particular, the most outstanding features of this algorithm are: 

• Problems at the surface as, although Pr vanishes Pt typically will not, 
implying that S{R) — oo. 

• The integrating factor is rather cumbersome as it involves both g and 
S. 

• However, by letting S{r) — 1 identically we easily recover the isotropic 
formula. 



1.3.6 Conclusions 

Algorithms for generating the solutions for both the isotropic and anisotropic 
cases have been developed. In the isotropic case we are able to get a solution 
by supplying g{r), while in the anisotropic case we have to supplement this 
with an additional function. If anything these algorithms work too well, as 
we can easily pick functions that lead to unphysical results. Specifically we 
would like to ensure 

• The energy density p is non-negative 

• The dominant energy conditions (p > \Pr\ and p > \Pt\) are satisfied. 



• The central pressure is finite. 
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• That the compactness x(r) < 1 so that a black hole is not sitting inside 
our star. 



With the exception of ensuring that the central pressure is finite, it is difficult 
to see how these conditions could be satisfied without picking g{r), calculat- 
ing the mass, density and pressure profiles and then explicitly checking that 
all of the above conditions are satisfied. In the next section we start looking 
at bounds which restrict the nature of "physical" solutions. The manner in 
which the bounds are obtained are reasonably general but have the disadvan- 
tage of not constructing valid solutions that saturate the bounds, or indeed 
demonstrate that the bounds can be saturated at all! Notice that we have 
to invoke energy conditions in the case Pt > Pr as an arbitrarily low central 
pressure could otherwise be obtained by letting Pf become high enough. 



It is interesting to note is that it is commonly claimed that the TOV 
equation ()1.46|) tells us that pressure contributes as a source of gravity, and 
so the pressure is higher than in Newtonian stars of the same density profile. 
This "regeneration of pressure" effect leads to a point where the gravitational 
attraction of the pressure is greater than the repulsion, hence the star is forced 
to collapse on itself. It is claimed that this mechanism is responsible for the 
Buchdahl-Bondi limit to be discussed in more detail in the next section. This 
argument is elegant and persuasive, but a glance at ()1.45|) shows us that it is 
seriously misleading. If the radial and the tangential pressures are allowed to 
be different we see that it is the radial pressure that acts as a source and that 
the tangential pressure (if positive) actually reduces the effect of gravity. We 
have also constructed an algorithm that uses Pt as free data, hence seeing 
pressure as contributing as an attractive source of the gravitational field is 
too naive. 
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1.4 Compactness bounds 

It has been well known for years that the compactness of a static spherically 
symmetric fluid x{R) = '^M/R is bounded above by 8/9, a result known as 
the Buchdahl-Bondi bound [HIE]- The result is significant as it shows that 
a star cannot get arbitrarily close to forming a black hole {x = !)• No such 
bound exists in the Newtonian case; although the concept of a black hole 
does not exist^ we can still require that the central pressure remains finite. 
The condition for hydrostatic equilibrium in a Newtonian star is 



dP / N / N / N rn{r) 4 

_ = _^(r)p(r), g{r) = — 
dr r 



—g{r)p{r), g{r) = — -— = -nrp, (Newtonian) (1.98) 



which is a linear differential equation for P{r). If the density is finite through- 
out the star then the average density and g{r) are also finite, hence for stars 
that have a finite radius R will also have finite central pressure. We can add 
a partial general relativistic correction onto g{r): 

dP / X / N / N "mir) + 47Tr^P(r) • n 

^ = -9ir)p{r), g{r) = ^2(1 _ 2m/r) ' (Post-Newtoman) 

(1.99) 

which still leaves a linear partial differential equation for P(r). The solution 
to this equation is 

P(r) = exp(-4./^-^^;iI^dr) (1.100) 



r(r — 2m(r)) \ J \ — 2m[x)/x 

It is now obvious that we cannot have xi"^) = 1 anywhere in the star, as the 
integrands will blow up. A calculation of a constant density star shows that 



^There are "dark stars" in Newtonian relativity where the escape velocity is at least 
the speed of light. However these stars are qualitatively very different; as the star can be 
a fluid free of singularities. While the light cannot escape to infinity, it will be able to 
propagate out to some finite distance Tc- Observers with r < Tc will see a perfectly normal 
star, and may be unaware that it is "dark" to distant observers. 
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the compactness can get arbitrarily close to one. The condition of constant 
density gives us the density and mass profiles easily: 

4 

p = po, rn{r) = -TrrVo- (1.101) 
The integrating factor can be expressed in terms of elementary functions: 

I J r^dxh.103) 



which implies 



= (I-X(r)f' (1.104) 

/ 3 \ 3/4 

= h--7rrVo] , (1.105) 



P(r) = (l-x(r))-^^/^|co-^^J^''r(l-^7rrVo)''^'dr 

= (1 - Xir))-'^' l^^o - f (1 - X)-'^' dxj (1.107) 
= Co(l-x(r))^/'-f (1.108) 

= ^o(l-^W^)'^'-f, (1.109) 

where the last equation makes reminds us that we cannot choose x(r) arbi- 
trarily as we are looking at a constant density solution. At the surface of the 
star P — and we have: 

p(i?) = o=(p. + |)(i-x)^/^-f (1.110) 

which implies that the compactness X — xi^) is given by 

\ 4/3 
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which shows that we may have the compactness as close to one as desired by 
allowing Pc/po to become large (but still finite). 

When we look at the full TOV equation with pressure as a "source" we 
obtain a non-linear differential equation. We will find that in the case of 
isotropic stars with a non- increasing density profile that the central pres- 
sure will blow up as the Buchdahl-Bondi limit is approached. This effect is 
sometimes called regeneration of pressure, although this term has to be used 
carefully (cf gTSU). 

1.4.1 Generalising Buchdahl— Bondi: > Pt 

In this section we consider a star subject to two restrictions: 

• The radial pressure Pr is never less than the tangential pressure Pt. 

• The averaged density defined by 

is a non-increasing function of r. 

The original derivation of the Buchdahl-Bondi limit stated that x < 8/9 for 
isotropic fiuids with non-increasing density profiles. That derivation is quite 
complicated and difficult to extend to more general results. More modern 
proofs, such as the proof given in Wald [79^, are very easy to generalise to 
Pr ^ Pt- As a smaller tangential pressure means that more radial pressure 
is required to achieve hydrostatic equilibrium, one would expect on physical 
grounds that the Buchdahl-Bondi bound would hold in this case as well. 

Rearranging expressions (fT"24b|l and (fT"^ and using the assumption that 
Pr > Pt we have: 



^<Pr-Pt) = --^\l^-—[-\ll-—C'\ +r[-) >0. (1.113) 
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The first term is bounded above: 



J ~ y rlrV t 



[I.IU) 



But here the left hand side is proportional to dp/dr which is non-positive 
by hypothesis. The bound ()1.114|) requires that the rhs is non-positive. But 
all the quantities outside of the derivative are positive, so we have the result 
that the derivative must be non-positive: 



d / 1 / 2m dC I ^ Q 
dr I r V r dr 



:i.ii5) 



As the quantity inside the brackets decreases monotonically with r, it is 
bounded below by its value at the surface of the star, r = R.^ i.e. 



R 



1 



2MdC 



R dr 



< 



r=R 



2m dC 
r dr 



:i.ii6) 



For all r > i? Birkhoff 's theorem tells us that the geometry is (locally) the 
Schwarzschild solution. The requirement that the first derivative is contin- 
uous allows us to evaluate the left hand side. After a rearrangement we 
obtain 



dr 



> 



M 



R^ 



2m 



CI 



r=R 



C\r=o > -Bi / > dr 

^ Jo , / 1 _ 2m 

C(o)< 



2M M 



R i?3 



=dr. 



2m 



:i.ii7) 



'1.1181 



;i.ii9) 



^In fact we know that if ri < r2 then ()i > ()2. But at the surface we know the 
geometry joins smoothly onto the Schwarzschild spacetime and can evaluate one of these 
bounds. 
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While we cannot perform the final integral as we do not know m{r), we can 
bound as we know the average density is non-increasing: 

m{r) > (1.120) 
which allows us to bound the integral in p.ll9j) 



Jl_2zn J, M-^-M^2 



The final integral can now be performed easily 



dr > / dr. (1.121) 



0<C(0)<^y/i-^-^. (1.122) 

The lack of horizons can only be satisfied if C > throughout the spacetime, 
which is impossible if x = IMjR > 8/9. When we consider the places where 
inequalities entered the argument we see that there are three places: initially 
Pr > Pt which is saturated by an isotropic fiuid, dp/dr < in obtaining 
()1.115|) which is saturated by p' = and the bound on the mass function 
m(r) > Mr^ / which is also saturated by p = 0. Thus the constant density 
or interior Schwarzschild solution saturates all these bounds and in this case 
we have 



N 3 / 2M 1 , , 

m = - IT - 2 

demonstrating that x < 8/9 is the best limit that we can achieve without 
invoking any energy conditions. 



To extend this result, we can integrate p.llTj) from r to R. This will give 
us information about ({r) anywhere inside the star. Repeating the above 
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calculations we find: 



C(r)< 



2M M 



R i?3 



-dr 



^ 2 m 

r 



2M 1 2M 



R 4 



1 / / 2M 



R 



u=l 



u=(r/RY 



2M /r \ 2 



1 - 



R \R 



(1.124) 

(1.125) 
(1.126) 



But the RHS is just C,{r) for a constant density star with p = SM/AttR'^. 
Denoting the constant density solution by C*('") we rewrite ()1.12(jj) as 



C(r) < C(r) 



;i.l27) 



for stars with non- increasing average density. This was done in i;54J, see 
appendix El A more significant generalisation would be to treat stars that 
have Pr < Pt everywhere, as is done in the ^1.4.2[ 



Improved bounds with the DEC. 

From equation ()1.117|) and equation p.l27|) we have 
C M r 



9[r) 



> 



> 



M 



C R^ C*{r)^l-2m/r '~ R^ C*{r)^l - x{r/RY C 



(1.128) 

Hence throughout the star the local acceleration due to gravity is always 
greater or equal to that of the constant density (interior Schwarzschild) so- 
lution. We can then use the field equation ()1.24b|) to show that 



lim87rPr(^) = hm 



2 , , 2m(r) 



r 



lim I -qir^ 

* 7 



Stt 



-pc 



(1.129) 



We note that we have g{r) > g*{r) and that pc > p* ■ Hence we have 

,2,.. Stt ,2 Stt 

lim87rP,,(r) = lini(-^(r)) - —p^ > \im{-g (r)) - — 

r^O r^O r 3 r^O T 3 



[1.130) 



^This last inequality follows from the fact that at the centre pc = Pc while p* = p{R). 
We then apply the condition that average density is decreasing: p{r) > p* . 
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or using ()1.69|) to expand g and g*: 

Pr(0)>P;(0) + ^(p*-pe). (1.131) 

This is actually a very weak bound. As shown in either appendix B or at 
the end of section 11.6.11 we can conclude that 

^'r(O) > P;(0). (1.132) 

Note that we do not have Prir) > P*{r) throughout the star, see appendix 
B for more information concerning this point. However, we see that if the 
interior Schwarzschild solution violates the dominant energy condition (DEC) 
P*(0) > p then all other decreasing density solutions will violate it too. The 
pressure profile for the interior Schwarzschild solution is given by 



P;(r) = - ^^^/^^l' - (1.133) 

3v/r^- V/1-X(r/P)2 



and in particular the central pressure is given by 



P;(0)=p*:i-=^^. (1.134) 
3vl - X - 1 

As the pressure decreases throughout the star, the maximum value is ob- 
tained at the centre. The dominant energy condition will be satisfied in the 
Schwarzschild interior solution iff 



^ ^ ^ <1, (1.135) 



3v/r^-i 

or in terms of a compactness, x ^ 3/4. This is not a maximum for all 

decreasing density stars as pc> p* ■ In general we need to satisfy 

p* 

— < 1. (1.136) 

Pc 

Using this condition and p.l33j) we have the result 

4[p./p1(2[p./p-] + l] p, 
^- (3[p,/p1 + l)2 ' p*-'- ^'-'^'^ 
Note that as the central pressure gets higher the Buchdahl-Bondi limit is 
obtained. This time we have only established that x ^ 8/9 is a bound for 
a solution that obeys the DEC, we have not provided an example which 
saturates the bound. With more work, better bounds may be available. 
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1.4.2 Bounds on compactness if P, < Pt 

When is less than Pt the tangential pressure is doing more to stop gravi- 
tational collapse than the radial pressure. Hence the required radial pressure 
is less than the isotropic case, and we would expect that even stars that 
have a non-increasing average density profile should be able to violate the 
Buchdahl-Bondi bound. By inspection of ()1.45|1 it is clear that we should be 
able to get the compactness as close to one as desired if we impose no upper 
limit on Pt. One limit we can place on the compactness will come about 
from insisting that the dominant energy condition holds: \Pt\ < p we 
are explicitly dealing with a type I stress-energy tensor of §4.3]. The bounds 
become difficult to saturate because the requirement of spherical symmetry 
means that Pt\r=o = -Pr|r=o and the condition Pt < p is in general difficult to 
satisfy. We make an attempt in 11.4.21 

To find more general bounds we look at the two measures of anisotropy: 

• Bounding the fractional change in the anisotropy, smoothed out by the 
energy density 



In the case that e is zero this reduces to the fractional change in 
anisotropy. 

• Bounding the anisotropy as a fraction of the energy density 



5p,,(e) := sup 



Pt + ep 
Pr + ep 



— 1 = sup 



Pt-Pr 

Pr + ep 



(1.138) 




Pt-Pr 



(1.139) 



= sup 



P 



Notice that because of spherical symmetry, we have Pj. = Pt a.t r = 0, 
implying that both 6p and Sp^{e) are non-negative. 
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Bounding by 5p^(e) 

If Sp^{e) exists and is finite, we know that 

as 5p,,(e) is the supremum. We can use the field equations ()1.24j] to find an 
expression bounding m: 

8n{Pt - Pr) < {Gff + eGii)6p,.{e). (1.141) 
Expanding out the pressures on the LHS: 

C) --0'<(^^^^ + ^-')^P.(^) (1-142) 



r / 2m / 1 A 2m 
(\ r I r V r 

which implies 



^ , ^ , m' m\ /m\' r 2m / 1 / 2m , , , 



The goal is now to attack each side independently, and make the LHS propor- 
tional to the average density. We are then able to get bounds for stars that 
have a non-increasing average density profile. The LHS is easy, if somewhat 
ugly: 

While this is not proportional to the derivative of the average density if 
e 7^ 1/3 it is easy to show via the product rule that if the average density is 
decreasing then m/r^+" is decreasing for positive n: 

m \' /m\' 1 /m\ n 

(1.145) 



Hence we can guarantee this term is negative ifO < e < 1/3 and if the 
average density is decreasing. 
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The RHS is a lot more messier, but it is fairly easy to show that 



By taking e < 1/3 we can make ()1 .14311 into the much nicer statement that 



The procedure is now very much the same as the > Pt case already 
discussed. By joining continuously onto the Schwarzschild metric we can 
obtain 



1 / 2m 



;i.i47) 



implying that 



1 / 2m , M _ X 

f V ^ ~ ~r ~ i?3+2%(e) = 2i?2+25p,(e) (1.148) 

As the density is decreasing we can again appeal to the fact that m(r) > Mj^ 
to be able to analytically bound this integral. Doing this we obtain 
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Here 2F1 is the Gauss hypergeometric function which takes on reasonably 
simple form when 5p^(e) is an even integer. I am also going to restrict atten- 



tion to the central value, in which case 



< C(o) < - 1 + 2 + • (1-153) 

Bounding by 5p 

As Sp is non-negative, we can use the Einstein equations to give: 



(1 I 2mdC\^^^ _^^d /m\ (1154) 
dr I r V r dr / ~ ** dr Vr^/ 



< (25p+ l)r3/(2^''+l)-2( ^^-3/(252+1) 



dr 



3^^-3/(2<5p+l)-l 

2Sp + l 



(1.155) 



V^i2Sp + ir<l. (1.156) 

Because in this case we know that the r exponent on the "density-like" term 
is less than three we cannot apply the simple result of ()1.145p to deduce that 
the LHS is decreasing. We can continue as the expression on the RHS is still 
simple enough to integrate by parts: 

R 



1 I 2m dC 
r V r dr 



<i / r'iv-D^fUL] dr (1.157) 
+ 3(1-7]) / —dr. (1.158) 

J r ^ 



m 

< — 



We can then find the bound for an arbitrary position throughout the star: 



2m dC I fm M, A 3(1 - n) m , 

-\1 - { - ^i^-v)] - ^dr. 1.159 

r V r dr T] \r-^ R'^ J r] Jr ''^ 

It does not seem clear from this that a simple bound can be placed on the 
compactness with knowledge of 6p, making it significantly less useful than 
the previous case. 
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Hypergeometric functions 

So far we have obtained the bound p.l5Hj) but this bound is not written in 
terms of elementary functions. The hypergeometric equation is defined to be 



, , X ab a(a + l)b(b + 1) 9 , ^ 

2Fi{a,b;c;x) = 1 + —x + ^ ^ / \, — -x^ 1-160 
c 2c[c + I) 

a(a+ l)(a + 2)6(6+ l)(6 + 2) 3 

^ 6c(c+l)(c + 2) 



a){a + l)...ia + i-l){b)...{b + i - 1) , 



c(c+l)...(c + i-l)i\ 

1=1 \ / \ / 



:i.i6r 



It is immediately apparent that the hypergeometric function is symmetric in 
its first two parameters a and b. The cases when either of these parameters 
are integers lead to tractable closed expressions for the Gauss hypergeometric 
function in terms of more elementary functions. In particular we have 



2Fi(^,l;2;x) = , , ~ (1.162) 



1 + v/r^ 



1 . „ \ 4 
3? 



2F1 ( 2; 3; X ) = — (2 - (2 + x) VI - X) (1-163) 



2F1 3; 4; = A(8 - (3^' + 4x + 8) v^T^) (1.164) 
2F1 Q, 4; 5; = ^(^28 - (40^^ + 48^' + 64^ + 128)^/1^). 



^1.165) 



These relations can be found by asking Maple, looking them up in mathe- 
matical tables, or using the relation 

2Fi(a, b- c- x) = (1 - Xf-^-' 2Fi{c -a,c-b;c;x), (1-166) 
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and then using the power series expansion in the new hypergeometric. We 
can also find the hmiting behaviour directly from the power series: 

/I \ 1 13 

lim 2F1 { -,l + n;2 + n;x] ^ I + -x + ^x^ + . . . (1.167) 
n— >oo y/ / 2 2! 

^ (1.168) 



These results can be put back into p.l53|) giving us bounds on the com- 
pactness. Let X* denote the value of x that saturates ()1.153|) in the sense 
that 



^ " ^* " 4(1+X,(e)) U' ^ ^ ^ ^ ^7 ^ °' ^^'^^^^ 

We can then find the bound on the compactness exactly for these few cases: 

*P,(e) = O^X-. = 5; (1-170) 

i.A') = 1 - X. = ^i^; (1.171) 

Me) ^ 2 ^ X. ^ 1 - (^i^^^i^l^^ (1.172) 



5 2 



2 



(4235 + 110 \/Tm)V3 3^ 

(5p,(e)=oo^X* = l- (1-173) 

A plot showing the points 5pr(e) = 0,1,2 as well a numerical solution of 
fll.l69j) and its tangent are shown in figure 11.21 Notice that if the anisotropy 
is unbounded then we cannot say anything about a star that is not a black 
hole, and this bound suggests that we can get arbitrarily close. Note that 
we do not currently have a way of obtaining a solution that saturates these 
bounds, so better bounds may exist. 



An attempt at saturation 

We see how close we can get to saturation of our bounds by looking again 
at a constant density star; this saturates the non-increasing condition. We 
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0.98 




Figure 1.2: The bound on the compactness X^^, as a function of the anisotropy 
parameter. 

would also like to let Pt = po, saturating the limit set by the dominant energy 
condition. This saturation would cause difficulties with the requirement Pr = 
Pt at the origin set by spherical symmetry, so we assume that the form of Pt 
is 



Note that the while the radial pressure profile may not be discontinuous there 
is no such requirement on the tangential pressure profile. Thus we can get 
arbitrarily close to saturating both our bounds. 

We begin with the anisotropic TOV equation p.45|) in the region r > e: 



The solution to this equation is given in terms of hypergeometric functions, 
but is not easy to interpret directly. We can learn quite a bit from inves- 
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tigating particular solutions. But first we rewrite (|1.176p in terms of the 
compactness: 

dPr (1 (mir) + ATTr^Pr)\ f m(r) + Anr^ Pr 2\ . , 

- —, ^ , XX + - n (1.177) 



dr \r r{r — 2m{r)) J \ r(r — 2rri(r)) r 
'2 {Mr^/R^ + 47rr^Pr)\ f 3x 



r r{r-2M{r/RY) J \87tR^ 
Mjr/Rf + Anr^Pr 2' 

r{r-2M{r/Rf) + r ' ^ ^ 

2 rix/R' + SnPr.)\ ( 3x 



r 2(1 -xr^/R^) J \8ttR^ 

,ix/R^ + S.P.)^2_ (1,179) 



2(l-xrVi?2) r 
At the surface of the star we know that Pr = 0. We see that 



dPr 



dr 



r=R 



R 2R{l-x) J \87tR'^ 
> if X < 4/5 
< if X > 4/5 



(1.180) 



Hence if x is less than 4/5 the gradient is positive and the radial pressure 
near the surface is negative. Gravity is not strong enough to hold the star 
together under the very strong imposed tangential pressures and the radial 
pressure must "help" . But then the second term in ()1.179|) is strictly positive, 
the only negative contribution begin the rx/i?^ piece in the numerator of the 
first term, which is maximised when r = R. Thus the pressure gradient is 
always positive and hence as we follow the profile inward the radial pressure 
gets more severely negative. 



If X > 4/5 then the radial pressure near the surface is positive and decreases 
to zero. For x ^ 0.92 the pressure eventually goes back to zero and then 
becomes negative, diverging as we approach the centre. For x ^ 0.92 the 
radial pressure is positive throughout, but diverges off to positive infinity 
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X = 0.9196 with P,=p 




Figure 1.3: The radial pressure of a star with compactness x — 0.9196. The second inter- 
cept is at « 0.07 and a constant density solution with compactness x* = 0.9196(0.07)^ w 
0.005 can be used in the region r < r*. 



as the centre is approached. To construct a solution that obeys the energy 
condition, one approach would be: 

• For 0.8 < X < 0.92 we have two places where Pr{r) = 0, denoted by 
r = (the interior solution) and r = R (the surface). 

• We use the solution Pt = p in the region r^, < r < R. 

• For r < we can generate an isotropic solution for a constant density 
star. It is equivalent to considering a constant density star of density 
Po and radius as Pr{ri,) = and we only require one initial condition 
for the solution. The compactness of this isotropic star would be 

This gives a legitimate solution provided that < 8/9. A solution is shown 
for X = 0.9196 in figure lOl 
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1.4.3 Theoretical limit on compactness: p' < 

We have shown that as Sp^{e) approaches infinity that x* approaches one 
()1.169j) . We have also shown that an exphcit example that obeys the energy 
conditions with x = 0.9196 in the last section. The question now is if the 
dominant energy condition places any bound on Sp^{e), hence ultimately a 
bound on x*- Note that it may not be possible to construct a stellar model 
with X = X*(sup[(5p^(e)]), but at least we know that x*(sup[(5p,,(e)]) will serve 
as an upper bound. 



From the definition of 5p^(e) we have 



Spje) 



sup 



Pr + ep 



where 



• 0<e<l/3 (Inequality condition) 

• \Pr\ < P) \Pt\ < P (Dominant energy condition) 

• p' < (Decreasing average density condition) 

• ^Pri^) ^ (Spherical symmetry requires Pr{0) = Pt{0)) 

As an extremely crude estimate, we may use the dominant energy condition 
to give us 



Spr (e) < sup 



Pt + Pr 



< 2 sup 



P 



\Pr + ep\ 



\Pr + ep| 

The denominator blows up for particular values of e if Pr is negative some- 
where in the star. Provided that Pr > — p/3 we can state that 



inf Sp(e) < oo, 

ee(0,l/3] 



;i.i82) 



and hence we can find some bound on the compactness. If we add the stronger 
hypothesis that > we have 

p-0 



inf 5p,.(e) < inf 

ee(0,l/3] <:e(0,l/3] + Cp 



3. 



;i.l83) 
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This allows us to conclude that 



X < X.(3) ~ 0.974. 



(1.184) 



It should be noted that to get this result we still have the restriction that the 
average density is decreasing outward and that the radial pressure is always 
positive. It should also be noted that our attempt with the constant density 
solution got us nowhere near this bound, suggesting that it is probably very 
weak. 

1.5 Pressure inside the star 

1.5.1 An upper bound on pressure if > Pt and < 

We can obtain both an upper and lower bound on the pressure inside at a 
radius r. As we are assuming that Pj. is greater than Pt throughout the star 
and that the average density is decreasing we can use ()1.115|1 




and integrating first from r to tq (r < tq) we obtain 




(1 



185) 



which in turn implies 




(1 



(1 



(1 



(1 



189) 



187) 



188) 



186) 
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Which gives a bound on ( at the centre of the star: 



C(0) < C(ro) 1 




2moC(ro)^ 




1 



2mo 



) 



) 



(1.190) 



The notation is by and large self explanatory, mo = m{ro) where tq is the 
arbitrary (but fixed!) radius. From ()1.35p we know 



so we find that the pressure is bounded above. On first sight this is not 
incredibly surprising; we have said that we have a given mass at a given 
radius, if the pressure became too large gravity would not be strong enough 
to hold the star together violating the stationary assumption we have already 
made. What is surprising is that this bound decreases as the compactness 
x(ro) = 2mo/ro increases and becomes zero as Xo 8/9. Physically the 
pressure increases as the mass of the star increases - what gives? 

We need to understand what this inequality is actually trying to tell us. If 
the pressure exceeds this value then we are going to have an infinite pressure 
by the time we reach the centre of the star. As x ~^ 8/9 then any radial 
pressure will result in infinite pressure at r = 0, hence why the bound is 
saturating. But we will also have to ensure that the pressure is large enough 
so that the star does not collapse, which would result in a lower bound. This 
bound does not, however give us any information about the central pressure. 

1.5.2 More bounds with the TOV equation 

We can use the TOV equation now to find a lower bound on pressure men- 
tioned in the previous section. This bound is more intuitive, as it comes 





Quite a but of algebra later we find that 
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about from requiring that the radial pressure is able to support the star 
rather than what pressure is required to stop a metric component changing 
sign. Of course, when the intersection of these bounds is empty we cannot 
hope to find solutions! 



dP m + Anr^P , , , , 

dr = - .^(i-2,„A) '''^^' 

mp mm' , „ , , 

< L = fi 194) 

- r2(l-2m/r) 47rr4(l - 2m/r) ' ^' ^ 



therefore 



PW - P(0) = -P(0) <-± r '^±(?^)' dr. (1.195) 



which implies 



^W^iE^W (1.196) 



>:^7 2M, — -— — +- — :/ —rrdr] 1.197) 

1 , / 1 M^+2\ 
> 7 2* 7^ 1.198 

1 1 (2M\ 

J— 

This holds for any spherically symmetric star, subject only to the constraints 
that the pressure and density are positive. We have not made any assump- 
tions about how the matter is distributed throughout the star. In the limit 
that the compactness approaches one, the pressure in the centre of the star 
becomes infinite to support the mass of the star from gravitational collapse. 
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1.6 Variations on a theme 

In deriving the generahsed Buchdahl-Bondi bound on the compactness we 
assumed that the density was non-increasing throughout the star so that we 
could bound an integral and that the pressure was isotropic. This bound 
happens to be saturated by the Schwarzschild interior solution giving an 
upper limit on the compactness of 8/9. 

By allowing an anisotropic pressure but still insisting that the density 
decreased as we went through the star it was found that the compactness 
could get arbitrarily close to one. It should be noted that in saying this I am 
assuming the existence of solutions of Einstein's equations that have constant 
Sp^{e) except for an arbitrarily small region near the origin, as Pr = Pt at the 
origin is required by spherical symmetry. By imposing the DEC we could 
claim that x* ^ 0.974. 

Kovetz jJU] concentrated his attention on isotropic fluids but no longer 
insisted that the average density decreased throughout the star. He looked 
at what would happen to the pressure profile as mass was moved around a 
perfect fluid without changing its total mass. This idea allows a comparison 
of the pressure profiles of two perfect fluid solutions with different density 
profiles. Do do this, Kovetz divides a fluid solution into "thin shells" and 
mass is taken from one shell at radius h + 5h and redistributed on a shell at 
radius a + (5a, where a < b. In the new solution the average density behaves 
like 



and from this Kovetz derives that the central pressure increases. Note that 
the total mass of the star remains unchanged so the variation in the (non- 
averaged) density cannot be non-negative or non-positive throughout the 
star. It is then argued that by successive use of thin shells that if we have 




>0 if a<r <b 



otherwise 
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any two solutions with pi > p2 everywhere but the same mass then the 
central pressure in star 1 is greater than the central pressure in star 2. One 
can strengthen this result: if pi > p2 everywhere then the central pressures 
obey the relation Pc,i > Pc,2- To do this we shall need to present Kovetz's 
argument in a slightly different way. 



Here the argument is presented in terms of a variation rather than a "thin 
shell" and Kovetz's result is rc-dcrivcd. The variational argument also allows 
for the change in average density to be non-monotonic, although no simple 
general statement about the central pressure can be given in this case. A 
simple argument for obtaining the central pressure increase result is also 
presented. 



1.6.1 A simple derivation for a "thin shell" 

As we are moving mass inward m(r) remains unchanged until we reach r — b, 
as m(r) is the mass within a radius r and does not depend on its distribution. 
A more formal derivation would be to use 5p is zero until r — b. The density 
p does not change except on the thin shells themselves, so it is certainly the 
same for r > b. Then the (isotropic) TOV equation for r > 6 

dPr ^ m(r) + Anr^Pr ^ 
dr r{r — 2m) ^ ' 

where there are no variations, so the rate of change of Pr is the same for r > b. 
As this is a first order differential equation it needs to be supplemented by 
an initial condition, a suitable one being that the radial pressure vanishes on 
the surface. But this is the same boundary condition for both solutions, so 
Pr is identical for r > b. 
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For definiteness's sake let us say that we are moving matter inward so 

mnew{r) > moid(r). Then 

/ s TTT'ne-w ~l~ 4:7rT Pr,new /-, rim \ 

'^'^ ^ r(r-2^_) ^'-'"'^ 

> ^oM + 4vrr^P,„ew_ .^ 202) 

r(r - 2moid) 

This translates into the following bound on the radial pressure gradient: 

< ^^;:''^r (i^,new + p). (1.203) 

dr r(r — 2moid) 

No effort has been made to distinguish the new and the old densities as they 
differ only on arbitrarily thin shells. Hence: 



dr ~ dr 



(1.204) 



and in particular 

Pr,new(&) " Pr,new(r) < P.,old(&) ' Pr,old{r). (1-205) 

But the pressure profiles are identical for r >b. Hence we have 

Pr-,new(r) > P.,oid(r) (1-206) 

with equality holding for all r > 6. To summarise: as matter is moved inward 
from a thin shell at r = 6 that 

• For r > b the pressure is unchanged. 

• For r < b the new pressure (and in particular the central pressure) is 
greater than or equal to the previous pressure. 

1.6.2 General variations 

We now generalise the results from Kovetz's thin shells. In Kovetz's orig- 
inal paper one mass profile was turned into another by moving thin shells 
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inward [50j. As a result more care was needed with the matching conditions 
so that one could find a variation in the pressure for each shell, and then 
integrate over the shells. Here we choose a slightly different approach, we 
allow for general small variations (but do not restrict ourselves to monotonic 
variations) of two solutions. 



We will write 

'nT'nir) = m{r) + 6m{r) (1.207) 

to mean the mass contained within a radius r after moving matter around is 
the same as that in our initial star plus some variation 5m(r). So far this is 
standard notation and we will prefix "5" in front of other quantities as well 
- no confusion should result. 



Variations in g[r) takes the following form: 

oF dm 



Aixr ^ 1 + STrr^P ^ 

■6P+^r-—- — --77(5m. 1.208 



1 — 2m/r r^(l — 2m/r)^ 

We remove references to the pressure (but not the variation in pressure) by 
expressing P = Pj. m. terms of g{r) using the Gff equation ()1.24bj) 

,- / N 47rr ^ 1 + 2q(r)r 

Sg(r) = SP + \ , 6m. 1.209 

^ 1 - 2m/r r2(l - 2m/r) ^ ^ 

The variation in the pressure is the expression that we are really interested 
in. By taking the first variations of the TOV equation: 

SP' = -{P + p)Sg - g6P - g6p (1.210) 
(P + p)) 6P - ^ ^] . {P + p)6m~g6P ~g6p. 



l — 2m/r J r2(l — 2m/r 

(1.211 
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We could eliminate pressure by substituting again from ()1.24b|) but it is easier 
to note 

m{r) + Anr^ Pr{r) 

9{r) 



r^(l — 2m{r)/r) 
4:7rrPr(r) 



+ ^ . — + 



l — 2m{r)/r l — 2m(r)/r \ 
inr'^Prir) + m{r)' 1 

r(l — 2m(r)/r) l — 2m{r)/r\ r 
4:7^r{Pr + p) 1 f m{r) 





(1.212) 


m(r)'\ 
r / 


(1.213) 


ir)\ 

.2 j 


(1.214) 




(1.215) 



1 — 2m(r)/r 1 — 2m(r)/r \ r 

We can drop the subscript "r" on the pressure as we are dealing with isotropic 
fluids. This allows elimination of P in ()1. 21111 



6P' = -(^2g-^ (ln(l - 2m/r))'^ 6P 



''^'"'^ ^9 + ,^^^]Sm-g6p. (1.216) 



Airr^ \ 1 — 2m /r 

This equation is a first order differential equation in 6P, so first we find 
the solution to the homogeneous equation to obtain the integrating factor. 
Setting 6p and 6m to zero in ()1.21()|1 leaves us with 

SP'h = -(^2g-^ (ln(l - 2m/r))'j 6Ph (1.217) 



a/1 — 2m/ r 



In I 1 1 SP„. (1.218) 



The homogeneous solution can be found by direct integration: 



(1. 219) 
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This allows us to get the general solution in the standard way: 

■5p + ■ ■ ■ 5m 



Jro 



SPh 

\/l - 2m/r / ((- 



2 



g{r) 5p{r) 



+ I 9{r) + /T./^ 1 '5m ^ dr. (1.220) 



y ^ yi - 2m/f 

1 +2g(f)f / w^y 

47rf3 V l-2m/r 



Here Tq is a constant of integration to be fixed by a boundary condition. By 
moving mass around it is possible that we change the location of the surface 
of the star. To avoid ambiguity I shall use R to denote the maximum radius 
of the two solutions. To evaluate Tq note that 5P{R) must be zero for both 
solutions as it is at the surface of one solution and either at the surface of 
the other solution or in vacuum, implying that Tq = is a suitable boundary 
condition. 



Let us take the first term in the integrand and integrate by parts: 



Sp{f)df=—I ' = {6m{f)) df (1.222) 

(5m(r) 



1 f arfgir) 



4:71 \ r^^l - 2m/r 



r'^^/l — 2m/r 



;i.223) 
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Hence 



_ g{r)Sm{r) a/1 - 2m/r / j C(0^(1 + 25'(r)r) 

The resulting integral is now in a form that we can work with, as we know 
the sign of Sm and can try to simplify the bracketed expression. Let us start 
with the differentiation and recall that C = 9 

(C7 = 2CC' = 2C'^?. (1.225) 
We can rewrite the final term in ()1.224|) as 

eg \ 2CV+CV 

r^i^l — 2m/r 





2m /r 


+ - 
V r 




) r2(l 




2m /r 


e 





(1.226) 



2/ + ^'-2^/r+f^ 



V r / 1 — 2ml r ^ 

(1.227) 

, , 2/ + o' - 2o/r - - ( ln(l - — ~ 

(1.228) 

As the pressure is isotropic we can use fll.H6|l to relate P and g' . 

2g' + g'-'^-^= ig' + + ^) + / - ^ (1.229) 
_ MP + {g + l/r)(m/r)' ^ ^2 _ 3^ 230) 



1 — 2m /r r 
2g{l — 2m/r) — 2m/r'^ + {gr + l){m/ry 2 3(7 
r(l — 2m/r) ^ r 

(1.231) 
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r ^ ^\ — 2m/r r^(l — 2m/r) 



'^\r) \ — 2m/r ' 1 — 2m/r 



+ 



4vr(p - p) 



or in terms of the original derivative: 
r^A/l — 2m I r j r'^^X — 2m/r 



(1.232) 
(1.233) 

47r (p - p) _ g_ 
r) \ — 2mjr ' l — 2m/r r 

(1.234) 



m\' 



'29 



+ 



Now we can go back to p.224|) 
C{rYg{r) 



C(r)2(l + 2(7(r)r) 



t^a/I — 2m/r 



g{r) + 



6m{r) dr. 



t^a/I — 2m /r 
Expanding out the derivative in the final term gives 

''-^^'^'^ W'""*'^' 4^C(r)= X rVl - 2m/f 



(m/r)' 
1 — 2m /r 

(1.235) 



X 



(1 + 2(?(f)r) f^^_^ ^ im/r)' 



r 



1 — 2m I r 



^2 I 2^(Ty I 4vr(p-p) , ^ 



5m(r) dr. (1.236) 



1 — 2m/r l — 2m/r 
Writing out the terms and cancelling allows for considerable simplification 



6P(r) + ^6m(r 
_ ^1 - 2m/r ''^ 



47rC(r)2 

a/1 — 2m/ r 
47rC(r)2 

a/1 — 2m/r 



R 



R 



C(r 




rVl- 


2m/r 


C(r 


? 


rVl- 


2m/r 


C(r 


? 


rVl- 


2m/r 



g{r 



7^^2 



47r(p - p) 



+ 



[m/r)' 



— Sm{r) dr 



l — 2m/r r{l — 2m/r) 

(girr - ^^^^ + ^I^iP^im 5m{r) dr 
\ l — 2m/r \ — 2m/r J 



g{r 



7^^2 



3(l-Tm/f) ,)^"^^^^^" 



(1.237) 
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Now we impose the condition that mass is moved inward or outward by 
requiring that 5m(r) be non-negative or non-positive respectively. Looking 
at the case where 6m{r) > we obtain: 

'nr) + ^^yra[r) ^^^^^^^ | r'^l-2m/r 

2 2m (r) /r^ 



> 0, (1.239) 

as all terms in the integrand are positive. In particular, as 6m{0) = (both 
solutions have no mass at r = 0) we obtain Kovetz's result 

5Pc > 0. (1.240) 



1.6.3 Variation in g{r) 

We can use the explicit expression for a variation in pressure express the 
change in g{r) by referring back to p.2()8|l : 

59 = 6P + J + 2,(r)r ^^ ^41) 

1 — 2m/r r^(l — 2m/r) 

-5m{r) + , „ , . — / o /. ^ ^{9 



r2(l-2m/r) - 2m/rC(r)2 Jr r'^^l - 2m/r 

^""^ U^(r)dr + ^^i^^^^MO (1-242) 



3(1 — 2m/r)j r2(l — 2m/r) 

l + ^(r)r ^ , r /"^ C(^) j 2 



r2(l-2m/r) y/l - 2m/rC(r)2 A " 2m/r 

+ ^-^^^^-^l-^^l^) dr. (1.243) 
3(1 - 2m/f) ] ^ ' ^ ' 

It is clear that if 5m is everywhere greater (less) than zero then 6g would also 
be greater (less) than zero. i.e. the local acceleration due to gravity increases 
as energy is moved toward the centre of the star. Note that we have no 
requirement that the average density decreases. If the average density does 
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decrease as we go out, then we know 6m{r) = m{r) — M{r/RY > and we 
recover ()1.128|1 as a special case. The integral formulae for 5P and 5g are of 
course still valid if 5m is not non-negative or non-positive corresponding to 
mass being moved around but not purely inward or outward. However there 
is very little one can say about this general situation. 

1.7 Conclusions 

Even within the confines of a spherically symmetric spacetime there is a lot 
of freedom for generating solutions. It is appropriate to reiterate a comment 
made at the beginning of Chapter 5 of that one can take any metric 
and define the energy momentum tensor so that the Einstein equations are 
satisfied. To restrict our attention to interesting solutions requires some 
prejudice on our part as to the nature of possible stress tensors. The attitude 
adopted in this chapter has been to ignore specific models of matter such as 
an equation of state but to make assumptions about the nature of any kind 
of matter. As noted in [T even the mild assumptions made in this chapter 
may be fundamentally flawed. It must also be admitted that such processes 
may become astrophysically relevant, as violation of the null energy condition 
is required for Hawking radiation to take place. ^ It is still hoped that the 
energy conditions describe bulk matter well, and outside of measurements of 
the cosmological constant no violations of any of the energy conditions have 
yet been found experimentally in hulk matter^. 

The goal of this chapter was to find some way of describing "reasonable" 
spherically symmetric solutions. The isotropic solutions were able to be 
generated by specification of g{r), the local acceleration due to gravity. By 

*The area increase theorem requires the nuU convergence condition: RabW°'W^ > for 
any null vector W. This is implied by Einstein equations {Rab oc Tab) and the null energy 
condition {TabW'^W'' > for null VF"). Hawking radiation requires a decrease in mass 
(and hence area) so must violate the NEC |41| . 

^The Casimir effect is one example of violations in the vacuum for small systems 
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the conditions of isotropy and spherical symmetry we only have one free 
function in the metric and the sceptical reader may ask why our method is 
of interest: after all does it not just formalise the idea of taking an arbitrary 
function, getting a spherically symmetric metric, generating a solution and 
then defining the energy-momentum tensor to be whatever comes out the 
other side in exactly the way mentioned in Hawking and Ellis? While in 
essence this is what happens we have the additional nice properties that 

• The gravity profile is physically measurable. 

• A class of physically unreasonable solutions can be discarded by re- 
quiring that g{r) > without discarding any physically reasonable 
solutions. 

• Another class of physically unreasonable solutions can be discarded by 
requiring that lim^-^o fi'('^)/''' = {Pc + Pc/'i)- 

• The solutions generated are guaranteed to be isotropic: Pj. = Pt- This 
is a good description of most stars throughout the bulk. 

The condition for a finite central pressure comes about from a single integral: 



and as the pressure profile decreases throughout an isotropic star we will 
then be guaranteed of a finite pressure everywhere. We can actually get a 
"quick and dirty" test of generating functions g{r) for finite pressure. We 
start by noting that 




(1.244) 





Using the lower bound, we have that ()1.244|1 is satisfied if 




exp I +2 / g{r) dr J dr < oo 



(1.246) 
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In contrast, the condition that p > and p > P are both exceedingly difficult 
to enforce. It is an unanswered question what the most general functions are 
that satisfy these constraints. It is unclear if such a project would even have 
a succinct answer. 

To this end other approaches have been developed, bounding the behaviour 
of the compactness throughout the star. The famous Buchdahl-Bondi result 
X < 8/9 is re-derived for stars that have a decreasing average density. To 
satisfy the DEC as well, we are able to show that the Schwarzschild interior 
solution requires that x < 3/4. A bound is established that limits to 8/9 
in the limit of the central density going to infinity for a star in which the 
average density decreases outward. However, no solutions that conform to 
the DEC have yet been found, and it is possible that better bounds exist. It 
is an unfortunate side effect of not imposing an equation of state that we are 
not able to obtain a maximum mass, only bound the compactness. One can 
see this by looking at the algorithmic solution for the mass ()1.78|1 : 
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{1 + rgy 



r(l + rg) 



We change to the dimensionless variables: 
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We begin by transforming d{r): 
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Hence our equation changes to 



Mm{x) 





x{l + Gx) 



G^ 




(1.247) 



m(a;) 



R / Gx^ f Gx 
M\{l + xGy V 



) 



(1.248) 





x{l + Gx) 



G^ 




The only term inside the braces that contains a mention of scale is an ar- 
bitrary constant of integration, and it becomes clear that we may scale our 
solutions to arbitrarily high masses. This is unfortunate, as masses are much 
easier to infer from the nearby orbiting bodies; the compactness requires the 
additional information about the size of the body. Hence we see immediately 
that general relativity does not impose a "maximum mass" of a star, and to 
find such information we do need to supplement our work with an equation 
of state. 

In the anisotropic case one has a lot more freedom. In this case the num- 
ber of important results seem rather limited. While one can construct an 
algorithm in this case as well it is far from clear that it is useful. The com- 
pactness bounds are harder to state as well: if one is prepared to ignore the 
DEC then it seems that the compactness can be arbitrarily close to one. If 
the DEC is taken into account and the radial pressure is assumed positive 
throughout and the average density decreases throughout the star then it 
has been shown that % < 0.974. An explicit example is constructed with 
X — 0.9196, and it is difficult to see how to improve on it, although not all 
constraints have been saturated in obtaining this profile. 



Perhaps the most important thing that we learn from anisotropy is to take 
the "regeneration of pressure" argument with a pinch of salt. 



Chapter 2 



Quasinormal modes, dirty black 
holes and Zerilli's equation 

2.1 Introduction 



Consider a stationary geometry such as a star or black hole in equilibrium 
and perturb it, say by infalling matter or a particle with a small impact 
parameter. For small perturbations we can linearise the geometry about the 
stable background solution, and the result is a dampened simple harmonic 
system. Hence the solutions are superpositions of the form: 



where g)J represents the value of Qah in the unperturbed geometry. The 
"resonant" frequencies a;„ are referred to as the quasinormal modes (QNMs) 
or ringing modes. 

In numerical relativity it is expected that one can model systems at late 
times accurately by using superpositions of a finite number of quasinormal 
modes over a spatially bounded region. Hence gravitational wave detectors 
are "tuned" to pick up frequencies of the order of the real part of the QNM 




(2.1) 
(2.2) 



= 9ab + {^9ah) exp(S5(u;„)t) cos(3?(a;„)t). 



{^{ujn) < 0), 
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with the smallest imaginary part, as they will decay the slowest. These 
are referred to in the literature as the "lowest lying QNMs". In looking 
at stellar structure, QNMs also determine the stability of the system under 
small perturbations. More detail can be found in the review article jlHl which 
discusses these applications. For all of these apphcations it is unsurprising 
that it is the slowly damped modes that are relevant. 

In contrast. Hod jlSI showed that if we take the real part of the highly 
damped modes of a Schwarzschild black hole to be a "transition frequency" , 
then it follows that the area of the black hole is quantised and some informa- 
tion on the area spectrum is available. Hod also showed that this argument 
could be used to determine the Barbero-Immirzi parameter 7, an arbitrary 
parameter in loop quantum gravity (LQG), making the theory essentially 
unique. For this conjecture to have validity we must find the same value 
of 7 for any black hole, not just a Schwarzschild black hole. The next two 
chapters are devoted to looking at the QNMs of spherically symmetric black 
holes surrounded by matter to see if the claim made by Hod is generic. 

Hod's QNM conjecture 

An investigation of the highly damped QNMs led to a proposal by Hod 
j43j : that the QNMs of a black hole given by classical general relativity 
tells us something about quantum gravity. The asymptotic spectrum for the 
Schwarzschild black hole was numerically found to obey |2l IS] : 



Recently this asymptotic form has been "proven" analytically by Motl and 
Neitzke |59t 160] . Hod's observation was that we could take the energy of a 
quanta emitted by a black hole to be: 




as n 



00. 
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via a "correspondence principle", where AM is the mass lost by the black 
hole and ip is the Planck length. But the horizon's area is determined by 
the mass with 

A^ie.(^)' ,2.5) 

which would imply that emission of such a quanta would be accompanied by 
a decrease in area: 

G^MAM , , 

AA = 327r . (2.6) 

Taking the change in mass from ()2.4|1 shows that the change in area is 

AA = 4£^ln3. (2.7) 

This suggests that the area of a black hole is quantised, and has an equi- 
spaced spectrum! More precisely, there is an area operator and its eigen- 
values are equally spaced. ^ The entropy of the black hole is given by the 
Bekenstein-Hawking formula: 

which has led to speculation that the entropy of a black hole is due to a 
"condensation" of 3-state systems at the horizon (such as spin-1 links of a 
spin network). 

This leads us to the QNM conjecture: 

The QNM conjecture is that the real part of the asymptotic form 
of the QNM of a black hole leads to a "transition frequency". 
This transition frequency is related to the dominant change in 
the area spectrum. For a generic black hole: 

OA hmun) 



^It is seen in i j2.2.2l that this conclusion is not automatic. It is possible to obtain a 
non-equi-spaced spectrum provided different assumptions are made. At the moment the 
arguments are heuristic only to provide motivation for the QNM conjecture. 



66 




Loop quantum gravity and the QNM conjecture 

The QNM mode conjecture is not a theory of quantum gravity, but an ob- 
servation. It says nothing about where the entropy comes from and does not 
give a formahsm for performing calculations. All that the QNM conjecture 
does is tie the QNM spectrum to the area spectrum of the black hole. LQG is 
a specific theory of quantum gravity, with one adjustable parameter 7 called 
the Barbero-Immirzi parameter. By making LQG compatible with the QNM 
conjecture, we arrive at a unique prediction for 7. This requires us to run 
the previous argument backward: LQG supplies the area spectrum, up to an 
arbitrary constant 7. The QNM conjecture fixes the absolute spacing in the 
spectrum and hence predicts a value for 7. 

An area spectrum {An} in LQG is given by 0: 



where rij is the number of spin-links with spin j. If we assume that almost 
all of the spin-links that pass through the horizon of a black hole have the 
same spin value jdom then the area is 



A spin j link on a spin network has 2j + 1 possible values. Hence we have a 
possible (2j + 1)^ states on the horizon itself, or an entropy of 



00 




(2.10) 



2j=0 




(2.11) 
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(2.12) 
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where the number of states has been approximated via ()2.1H) . This can 
be made compatible with the Bekenstein-Hawking entropy formula provided 
we choose 

7 = ^f''- + . (2.13) 

27rA/jdom(jdom + 1) 

There is some theoretical prejudice for jdom = 1/2, as fermions would fit 
neatly into such a theory. Instead, we now repeat Hod's argument: a loss 
of a spin jdom link would emit energy hM{uj). This would imply a change in 
area is given by ()2.11|) 



AA = 87r7£|Vjdom(jdom + 1). (2.14) 

Inverting the argument in ()2.4|) implies that the "transition frequency" should 
be given by 

_ AMc^ _ AAc^ _ c^AA 
^""^ ~ ~ 32nhGm ~ 327rGMiy ^^'^^^ 

Substituting in for the from (j2l4|l : 



87r7£|,v^Jdom(jdom + 1) 3 7Vjdom(jdom + 1) , 

= — 32^^Gm — " = — — ■ ^ ^ ^ 

Equation ()2.13|) constrains 7 in a LQG model to recover the correct expression 
for the entropy. We then find that 

ln(2jdom + 1) ^ „^ 

= 8nGM • ^2.17) 

This is the real part of the asymptotic spectrum ()2.3|1 if we choose jdom = 1' 
Hence the QNM conjecture and LQG are not only compatible, but forcing 
such compatibility uniquely determines 7. 



Issues with the QNM conjecture 

While the QNM conjecture is interesting, it relies on a numerical coincidence 
and it leaves many questions unanswered. We have seen that the QNM 
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conjecture relies on each spin link of spin j adding the same amount to the 
area (at least asymptotically), it does not rely on there being an equi-spaced 
area spectrum for different values of j. Indeed, the standard LQG spectrum 
provides an explicit counter-example! Provided that one spin dominates and 
each of those dominating spin links providing an equal amount to the area 
allows the black hole to decrease its area in multiples of a fixed amount, and 
satisfy the QNM conjecture. However, the prediction that the black hole 
loses area (and hence mass squared, in the Schwarzschild case) in constant 
discrete amounts is at odds with the emitted radiation (Hawking radiation) 
being purely blackbody. The prediction that the black hole should be a 
condensate of links with spin 1 also raises the question of how fermions are 
incorporated into LQG. Why does this correspondence only occur for the 
rapidly decaying modes (i.e. large n)? These are interesting questions and 
will be addressed in ^2.21 

The research in the next two chapters is more modest. All of the fore- 
going considerations were focused on a Schwarzschild black hole - that is a 
spherically symmetric black hole in an otherwise empty spacetime - but if 
the QNM correspondence is fundamental it should apply to all black hole 
spacetimes. As the general black hole entropy is formulated in terms of the 
surface gravity n, the QNMs should also rely on n. This is significant as k, is 
not completely determined by the black hole, but is affected by (non-black 
hole) matter outside the horizon. Hence we introduce the notion of a dirty 
black hole as a black hole surrounded by matter, and investigate if the QNM 
conjecture is still valid. In terms of LQG specifically, we know that 7 is a pa- 
rameter that can only be fixed once; that is, two different geometries should 
not give different values for 7. If one enforces LQG only to be consistent 
with the Bekenstein-Hawking entropy then different geometries do give the 
same value for 7 fl2.13|l . although jdom is still undetermined (see §B of j^l). 
Our task will be to see if the QNM conjecture will predict jdom = 1 for all 
static, spherically symmetric black hole geometries. 
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The next section of this thesis outhnes in more detail some of the arguments 
for and against the QNM conjecture. Readers famihar with the arguments or 
wanting to see concrete resuhs would be best to look at §2.H[ where QNMs are 
described as a scattering problem, and then at ^2. 41 where the wave equations 
for dirty black hole spacetimes are derived. In chapterElthe different methods 
of finding quasinormal modes are discussed. 

2.2 The first window: Alice's curiosity 
2.2.1 Black hole entropy 

It has long been known that one can associate a temperature and entropy 
with black holes, as was initially advocated by Bekenstein jHj. One argument 
is that black holes have many features reminiscent of thermodynamic quan- 
tities. In particular, the area of a black hole horizon plays an important role, 
as can be seen by considering the following results: 

• Reversible processes leave the area unchanged, a result due to Christo- 
doulou 120]. 

• The area theorem, proved by Hawking in 1971, showed that subject 
to the null energy condition the area of a black hole never decreases. 

mun]. 

The area theorem trivially contains Christodoulou's result, as the initial and 
final areas are the same for a reversible process and the area never decreases. 
Despite this redundancy, Christodoulou's result nicely emphasises the idea 
that the area of a black hole is an adiabatic invariant. Stationary black holes 
can be described by a limited number of parameters - mass, charge and 
angular momentum - and all of these can be changed by reversible processes 
|2()j . Hence if a black hole has entropy, it is plausible that: 



^-BH = S{A). 



(2.18) 
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At this stage we only know that S{A) is a monotonically increasing function 
of area. To obtain the functional form of the entropy, Bekenstein argued that 
truly reversible processes were impossible: black holes had to increase their 
area and that the amount of increase was independent of the parameters of 
the black hole! To see this requires Christodoulou's result that the increase 
in the area of a horizon, when a black hole captures a particle of mass m 
with a turning point at proper distance b away from the horizon, is given by 
[UliSI: 

Tfib 

(AAU^ = SnG —. (2.19) 

The process is then reversible if and only if 6 = 0; that is, the turning point 
occurs on the horizon. Bekenstein argued that because of the uncertainty 
principle, the turning point for a real particle could not be determined ex- 
actly. Taking a particle that classically had a turning point at the horizon 
would be subject to some uncertainty in both the momentum and position, 
leading to 



> 47r^ = ATril. (2.21) 

This minimum area increase is completely independent of the characteristics 
of the black hole! This led Bekenstein to the idea that the entropy was 
proportional to the area and that the area was discrete with equi-spaced 
eigenvalues. 

Bekenstein's original suggestion of an entropy associated with black holes 
was motivated by a very different consideration: ensuring that the validity of 
the second law of thermodynamics was not violated. In fact the suggestion 
pre-dated both the area theorem and Christodoulou's result! If black holes 
did not have an entropy, then by lowering objects into a black hole we could 
"lose" the entropy associated with that object and hence violate the second 
law of thermodynamics! However, while this argument suggests that an 
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entropy is associated with a black hole, it does not tell us what form it 
takes. The result that the entropy is proportional to the area leads to an 
interesting problem - that the entropy of an object is bounded by the amount 
of area it would have if made into a black hole. This is the idea behind 
the holographic principle. Unfortunately, naiVe attempts to create a bound 
on the amount of entropy an object can have that is compatible with the 
holographic principle are plagued with difficulties, and it is an ongoing task 
to find a satisfactory statement. However, if we insist that the total entropy, 
or so-called generalised entropy, never decreases, we will require some sort of 
holographic principle. 

Notice that neither of these arguments suggest any sort of "microstates" of 
a black hole that the entropy is associated with. This is another unresolved 
problem in black hole thermodynamics which will be re-addressed in chapter 
il 

2.2.2 Discussion concerning the QNM conjecture 

There are many unanswered questions about the QNM conjecture that leave 
it somewhat unsatisfactory. The most important such question is whether 
or not the conjecture is generic, or just a numerical coincidence for the 
Schwarzschild black hole. Even if the result turns out to be generic, we 
still have to explain how the discrete spectrum from the QNM gives rise to 
blackbody Hawking radiation. We also have to fit fermions into a theory that 
seems to be dominated by states with a three fold area degeneracy^. Many 
of these questions have been raised by Polychronakos |(i2j . 

Let us first introduce the issue of an the logical relationship between the 
QNM conjecture, LQG and quantum gravity. The QNM conjecture is non- 
generic as it requires the area of a black hole horizon to be quantised. This 

^These are "spin-1 links" in LQG, but the QNM conjecture allows for any three mi- 
crostate subsystem. 
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quantisation is not required by a quantum theory of gravity; there exist 
"perfectly good" string theories that do not quantise area.'^ In LQG the basic 
variables are holonomies which describe how to parallel transport spinors, 
and the fundamental excitations of the geometry is analogous to electric field 
lines, or quantum vortices. In the classical limit the geometries approach 
the standard area (see for a review, or jMl IHH IHH] for calculations of 
the area spectrum and the relationship of these "vortices" to spin network 
states). Each vortex line penetrating a surface then contributes to its area. 
The vortex lines are formally analogous to spin network states and can be 
labelled by the half-integers. Although the QNM conjecture is most often 
used to determine the Barbero-Immirzi parameter in LQG, it should be noted 
that the conjecture and LQG are logically independent. In particular: 

• If the QNM conjecture is false then it does not tell us anything about 
the nature of quantum gravity, and in particular about LQG. 

• If LQG is false, then we have discarded a specific model that quantises 
area on the horizon. The QNM conjecture could still be compatible 
with another theory of quantum gravity that quantises area. 

• If the area is not quantised, then of course LQG and the QNM conjec- 
ture are both false. 

Once we have assumed that the area is quantised, we must decide on an 

appropriate area spectrum. As seen in the introduction, the QNM conjecture 

implies that (for a non-rotating black hole) the emission of any quanta gives 

the same change in area. In a generic spectrum this would imply that one 

particular transition dominates. In an equi-spaced area spectrum transitions 

would could occur between the n and the n — c states for any n and some 

constant c. In principle it is also possible for different transitions in a non- 

equi-spaced spectrum to produce the same change in area via selection rules. 

■^Here perfectly good string theories means that string theories that do not quantise 
area suffer no more technical problems than those that do. 
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In LQG there are two spectra that are frequently discussed. The first is the 
standard spectrum (c.f. ()2.1()j) ): 



„2 <^ I ^ 



where Oj is the area from a spin hnk with spin i/2. The other comes from 
applying "radiative corrections" on ()2.10|) (see .IJ) 

ai = A'K-ffp{i + l). (2.22b) 



One of the arguments presented in Polychronakos's paper joz] was that 
fermions could be included in LQG, while still having the dominant contri- 
bution to the area due to links with spin one, thus remaining compatible 
with the QNM conjecture. More precisely, a theory of LQG with gauge 
group SU{2) is still compatible with the QNM conjecture. Polychronakos 
demonstrated that the dominance of spin one links could occur due to ther- 
modynamic suppression of spin half particles. 

To see how this comes about, consider a horizon of total area A. Denote 
the number of spin links that pass through the horizon with spin i/2 by rij. 
Then it follows that 

oo 

^ = ^n,ai, (2.23) 

where is the area added by each link with spin i/2. The details of the 
area spectrum are left open at this point in the argument, allowing us to 
compare (|2.22a,j) to ()2.22b;i later on.^ Each set {rii} satisfying fTT^ is called 
a partition of A. A link with spin i/2 has a degeneracy given by 

gi = t + l. (2.24) 



■^The most general spectrum would be a function of all the n^, and would not require 
that each spin link of the same value would add the same amount of area. In the case 
where the horizon area is some arbitrary function we cannot relate the QNM conjecture 
and LQG. However, the purpose of this argument is simply to show that with the inclusion 
of fermions, LQG and the QNM conjecture are not mutually inconsistent. 
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Here it is assumed that the area added by a hnk is given by the total spin j 
and all of the 2j + 1 internal states are independent. 

Now a partition function for the area is constructed [HS IBBj 

oo 

Z{P) = J2Il^' + 1)"' exp(-/5m,a,). (2.25) 

{nii} i=0 

Explicitly we are summing over all sequences of positive integers - not only 
partitions of the area A. Here we are looking at a grand canonical ensem- 
ble, rather than a micro-canonical ensemble. Note the absence of factorials, 
implying that the spin-links with the same value of j are not indistinguish- 
able. Note also that this distinguishability (or partial distinguishability p2]) 
is an assumption that allows fermions to be incorporated, and has not been 
"proven". For standard thermodynamic arguments to be valid it is neces- 
sary to assume that the black hole is enclosed and that the total area of the 
system is conserved. As pointed out by Gour and Suneeta there is no 
natural analogue of a heat bath in this system. Instead we rather recklessly 
assume that standard thermodynamic arguments hold.^ 

The rest of the argument is nicely given by Gour and Suneeta, given the 
assumption that we can use thermodynamic arguments. If we have (3 less 
than a critical value (3c then higher spins dominate the partition function 
and the partition function does not converge. In particular, the contribution 
from a given spin i/2 will not converge unless (3 > jSi, where (3i is given by 

A = = i^^^. (2.26) 

To see this, note that the sequences of positive integers has the natural num- 
bers as a subsequence. Thus the partition function has the geometric series 

^Note that the assumption that "area thermodynamics" is vahd is weak, but this has 
nothing to do with standard black hole thermodynamics which are provable statements in 
classical general relativity. 
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as a subseries, and ()2.26|) ensures converges of this particular subsequence. 
To ensure all the subsequences converge it is necessary that 

/3 > /3c = maxA. (2.27) 

i 

Provided [3 > I3c^e have 



=0 



As the area is the analogue of the energy for this "thermodynamic" system 
we have 

_ dln(Z) _^ (i + l)a^exp(-/ja^) 
^ ' d/3 Z^i-(, + i)exp(-/?a,)' ^ ' ' 

The numerator vanishes rapidly, provided that the denominator is well be- 
haved. Thus large areas are possible if /? ~ /3c- 



Nothing more can be said until an area spectrum has been assumed. Using 
the standard area spectrum ()2.22a|) we have 

lim/?,=0 (2.30a) 

In 2 0.0318 

A = = — ^ 2.30b 

In 3 0.0309 

= 7= = — ^ 2.30c 

In 4 0.0285 
= /— = — (2.30d) 



All subsequent spins have a lower /3i, thus the standard spectrum predicts 
domination by spin 1/2 particles. The prediction for 7 does not agree with 
the QNM conjecture. 
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However, if one uses the equi-spaced spectrum ()2.22b|) : 



Inl 

In 2 0.0276 



Po = j—^ = (2.31a) 



In 3 0.0291 
lii4 0.0276 

with all subsequent /5jS lower. Then the dominant contribution is from i = 2, 
or spin links with spin 1. Note that the fermions are still described by the 
theory but are suppressed. 



Although there are a number of aspects of this treatment that are at best 
unsatisfactory, it provides a concrete example of how LQG with fermions 
and the QNM conjecture can be compatible. Before taking the argument 
too seriously, it would need to explain why the area can be treated as a 
conserved quantity. More seriously, objects that are not black holes tend to 
have larger areas than a black hole of the same mass; so the above argument 
would suggest that everything macroscopic was a collection of spin-1 links! 

Polychronakos has also looked at the possibility that black hole states 
are fully distinguishable. Assuming the QNM conjecture, the Bekenstein- 
Hawking entropy cannot be (exactly) recovered in LQG without changing 
the gauge group to SO (3) and hence excluding coupling to fermions. 

A separate issue facing the QNM conjecture is Hawking radiation. While 
the conjecture is not sensitive to the details of the area spectrum it does tell us 
that the radiation from the black hole will be a discrete spectrum rather than 
a continuous (blackbody) spectrum predicted by Hawking. Polychronakos 
[n2] states that this may be analogous to an atom in thermal equilibrium 
with a heat bath which would radiate its own line spectrum if placed in 
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isolation. The problem with this analogy is that blackbody radiation of 
atoms is derived when many "independent oscillators" are present, whereas 
black hole temperature calculations assume the metric of a single black hole. 
Only at high frequencies does the Hawking radiation "look" like a black body 
spectrum. 

The condition that the QNM conjecture is generic is the most serious and 
also the hardest to check. The goal of the next two chapters is to investigate 
dirty black hole spacetimes, and try to determine if the "correspondence 
principle" that applies to the Schwarzschild black hole applies to these more 
general black hole spacetimes. 

2.3 An introduction to quasinormal modes 

The idea of a quasinormal mode was presented very briefly in the introduction 
of this chapter. Here, we extend that idea, and show how QNMs are defined 
and what they are used for. Of direct relevance to the investigation of the 
QNM of black holes will be the interpretation of QNM as a scattering problem 
in quantum mechanics. This motivates the "monodromy" approach of Motl 
and Neitzke jHOj and the "Born approximation" of Visser, Medved and Martin 
|56j . In addition, a brief overview will be given on how quasinormal modes 
are used in other of astrophysics. 

2.3.1 Defining quasinormal modes 

The following is a modified discussion from Kokkotas and Schmidt [4Sj. Let 
us consider a wave equation with a potential function that vanishes provided 
that > Xq: 



+ V{x)x{t,x) = Q. 



(2.32) 
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When normal modes are studied we restrict attention to systems of finite 
energy. In this case the operator H, defined by 



is self-adjoint and has a continuous spectrum if the domain of x and t is 
infinite. Because of the continuous spectrum, we no longer have a square- 
integrable eigenfunction |1H]. By introducing distributions we can still write 
superpositions of solutions as new solutions. 

Normal modes are useful in studying steady oscillations of a system and 
non-localised finite energy solutions. They are not suitable for defining the 
evolution of data of compact support or systems that radiate energy off to 
infinity. Notice that these problems only occur in infinite systems where the 
spectrum is continuous. Even finite systems can be dissipative, and here 
normal modes do not work either. In modelling perturbations, we want to 
treat some initial perturbation that is generically of compact support and 
that will radiate away to infinity; for this we use quasinormal modes. To get 
decay of the solutions we require that operator H is no longer self-adjoint, so 
that we may have complex eigenvalues. Hence we cannot look at the space of 
finite energy solutions! The physical reason for this is that the wave equation 
is reversible, so if we have a finite energy solution at t = that decays at 
late times it must grow at early times. This is not a problem physically as 
something would cause the initial perturbation; it would not be a decay of 
an infinite "perturbation" at an infinite time in the past! 

We are assuming that the potential vanishes if |x| > xq, and is always 
positive. Then the solution xi^i ^) is always bounded above. Define the 
Laplace transform of x{^i^)'- 



H 



+ V{x), 



(2.33) 




(2.34) 
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which makes sense for all 3fJ(s) > 0. The wave equation ()2.32|) becomes 

^'XL -^ + VXL = sxiO, x) + ^(0, x) (2.35) 

where the right hand side contains all the initial data. One approach to 
solving such a problem is to use Green's functions. To find a Green's function 
one looks at solutions that give a "delta function" contribution. A well known 
method is to find two linearly independent solutions to the homogeneous 
equation /+ and /_. A Green's function can then be explicitly found 



^, . ^ , f-is,Xi) U{S,X2) iXi<X2) 

G(s,Xi,X2) = X <^ (2.36) 

W^(^) /_(s,X2)/+(.,Xi) (X1>X2) 



where W{s) is the Wronskian of the two solutions: 

W{s) - 



f- u 

9f- df+ 



(2.37) 



dx dx 

After finding a Green's function, the original differential equation is solved, 
albeit in integral form: 

X{t,x) = j G'(s,xi,a;2)|sx(0,a;) + ^(0,x)| ds. (2.38) 

In selecting /+ and /_ it is convenient to pick them such that they are 
bounded as x becomes large. As they are solutions of the homogeneous 
equation and the potential is assumed to vanish, we have 

/ = exp(Tsa;), (|a;| > xq) (2.39) 

Hence we pick so that it decays for large x and /_ so that it decays at 
large negative x: 

f±{x) = exp(=Fsx), (±x > Xq). (2.40) 

Although the large distance behaviour of /+ is exp(— sa;), toward small values 
of X, /+ is represented by a linear combination of exponential terms. This is 
similar to the tracking of WKB solutions used in the monodromy approach 
to QNM. 
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Im(s) 



Re(s) 
► 



Figure 2.1: Taking the inverse Laplace transform. The poles in the Green's 
function are denoted by "x" and the contour of integration is the dotted line. 



We finally have enough background to define the QNM themselves! Sup- 
pose that there is a (complex) frequency Un such that: 

U{x) = c{s)f^{x) (2.41) 

where c(s) is a complex number. This would imply that f+{x) tends to one at 
large positive x, and to c{s) at large negative x. In particular the Wronskian 
of the two solutions vanishes and the Green's function becomes singular. 
The fact that the solution is bounded means that the Green's function exists 
and is analytic for 3fJ(s) > 0, so that any poles if they occur must occur for 
3fJ(s) < 0. We call such values of s quasinormal and they form a countable 
set for positive potentials of compact support OEHI, hence we label them 
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We can then take the inverse Laplace transform or Bromwich integral |64j : 



where a > for convergence reasons. By closing the contour as shown in 
figure 12.11 and using the Cauchy residue theorem, we see that it is the poles 
of the Green's function that contribute to the line integral. The solution in 
a finite region of space (xq < x < a;i) and at late times (t > ti) can be 
approximated by the finite sum: 



where the s„ are quasinormal frequencies ordered by their real parts: 3?(a;„) < 
3fJ(s„+i). The coefficients dn depend only on the initial data. 

So how good an approximation is ()2.43|) ? Bacholet and Motet-Bacholet [H] 
showed that for positive potential of compact support, 



Xit,x) -XNit,x)\ < C(xo,Xi;ti;e)exp((-|3fJ(s^+i)| +e)t) (2.44) 



where e is greater than zero. At late times the approximation ()2.43j) becomes 
very close to the exact numerical solutions, making QNM a very effective way 
of studying perturbations. 

In models of black holes the "potential function" V{x) does not have com- 
pact support and this introduces additional complications. The functions 
/+ and /_ may not exhibit exponential decay as large values of x are ap- 
proached. The longer tails of the "quasinormal eigenfunctions" means that 
the decay of the error ()2.44|) is also less than exponential. As pointed out by 
Kokkotas and Schmidt |1H|, if one restricts attention to a finite interval of 
time Iq < t < ti then the wave equation cannot depend on the values of the 
potential for large x and exponential decay of the error is recovered. 




(2.42) 



N 




(2.43) 



n=l 
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2.3.2 Relating QNM to scattering 

The solution tplx) to the time independent Schrodinger equation 

-IrS + ''<^W' = E4, (2.45) 
will also solve the homogeneous Laplace transform of the wave equation 



+ V{x)i) = -s-^ip 



(2.46) 



if we identify E = —s^ and m = ^/h/2. Under this identification, the two 
systems are equivalent. From the previous discussion we know that QNMs 
are restricted to 3?(s) < if V{x) is positive and of compact support. In 
terms of E, this means that any value of E is allowed except for values lying 
on the negative real axis.® 



We can show that negative real E is not allowed directly from ()2.45j) . Let 
us multiply ()2.45|) by the complex conjugate ip* and integrate over the real 
line: 



OX'^ 



V{x)i)*i) dx = E dx = E 



(2.47) 



where the wavefunction is taken to be normalised. Applying integration by 
parts to the first term gives 



dx = ip* — — ■ 

dx'^ dx 



dip* dip 
dx dx 



dx 



POO 


dip 


1 —oo 


dx 



dx. 



Hence the Schrodinger equation becomes 



dih ^ /"°° 
/ -- dx+ / V{x)\ip\^dx = E 

J —CO J —OO 



(2.48) 
(2.49) 

(2.50) 



^Such energies are a requirement for bound states. Note that a negative energy eigen- 
value does not necessarily imply a bound state as the state must also be normalisable to 
qualify as a bound state. 
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V(r) 



V(r) 



Incident wave 





Reflected wave 



Transmitted wave 



A. 







Figure 2.2: If particles are only incident from one side, then a clear separa- 
tion between the reflected and transmitted particles can be made (left hand 
figure). If we have particles incident from both sides, then the ingoing am- 
plitudes are still incident waves (Ain and -Bin), but the outgoing waves are a 
superposition of transmitted and reflected waves. 

Positivity of the first term implies that 



which contradicts the assumption that V{x) is never negative. Note that 
values of E on the negative real axis are permitted if V{x) < anywhere, 
although these values can only be associated with finite energy on a point 
spectrum. If a negative eigenvalue En does exist, then solutions can grow 
exponentially in time. 

As the wave equation is a linear differential equation in ip (or x) then the 
solutions linearly superpose. We can place purely left travelling waves off 
at infinity, by interacting with a potential these will turn into a mixture of 
transmitted and rcficctcd amplitudes. As the wave equation (or Schrodinger's 
equation) is second order we have two linearly independent solutions. Let 
us choose one solution that is left travelling when V{x) = denoted iJjl, the 
other purely right travelling denoted ipR. Then we have in general 




(2.51) 




(2.52) 



in regions where the potential vanishes. Let us analyse waves coming in 
from the left only and hitting a potential. Toward the left the solution will 
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be a superposition of reflected waves and the incident wave. The boundary 
condition of no incoming radiation from the right imphes that the solution 
will be purely right travelling. Only one boundary condition has been fixed 
so far, let us fix the incoming amplitude as A^^. Let us assume that the 
potential vanishes once > xq. Then the solution can be found: 

, , , J Ai^^L{x) + A^^ti^Rix) X < -Xo 
^(x) = < . (2.53) 

[BoutipR{x) x> Xo 

This is not equivalent to solving the problem as it does not tell us about the 
solution inside the potential. It is often the case that we are only interested 
in what happens far away from the potential. We can find the behaviour 
a large distance from the potential by finding the so-called reflection and 
transmission amplitudes. 

(Incoming from the left only). (2-54) 

Bout — Ta Ain 

If the "energy" is real (or s is purely imaginary) then we have the "conser- 
vation law" : 

\Ra\' + \Ta\''^1. (2.55) 

For obvious reasons, Ra and Ta are called the (left) refiection and transmis- 
sion amplitudes. Analogous relations exist if there are only incident waves 
from the right: 

^ '° (Incoming from the right only) . (2.56) 

Bout — Tb -Bin 

It is more difficult to give the interpretation of refiection and transmission as 
any outgoing radiation will be a superposition of refiected and transmitted 
parts. However, the principle of superposition makes finding the solution 
very easy: 

(2.57) 

The matrix here is the familiar S-matrix from quantum mechanics. 
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We know that as x becomes large (small) ^(x) f+ix) f-{x)) as 

it is the well-behaved (decaying) solution in this regime. At a quasinormal 
mode the functions f+{x) and f-{x) become linearly dependent, and we can 
have 

^{x) cx f+{x) (2.58) 

with well behaved solutions at both large positive and negative x. But these 
solutions are purely outgoing: A-^^ = B-^^ = 0! In terms of the S- matrix we 
have 

Bout J \Ta Rb ) \ 

Hence at a QNM the S-matrix becomes singular! We can also tell that QNM 
modes occur only for purely outgoing waves at infinity and the horizon. 

This similarity with scattering problems in quantum mechanics means that 
we can try and adapt techniques in quantum mechanics to finding QNM. The 
techniques that are particularly useful are modified WKB approaches and 
Born approximation approaches as they are particularly sensitive to poles in 
the S-matrix, and hence in the Green's function. 

Some care is warranted as it not only QNMs that correspond to poles in the 
S-matrix. Bound states of potentials where nothing propagates out to infinity 
also cause poles in the S matrix. For purely positive potentials this is not an 
issue, but bound states can develop if the potential is negative anywhere. We 
have also made a lot of use of the fact that the potential vanishes, so that we 
can talk about "out" states without having to be too careful about whether 
or not a sensible meaning can be given. The cases dealt with in black hole 
physics are not of compact support, but decay exponentially. It is assumed 
that "outgoing" states still make sense. 
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2.4 Wave equations for dirty black hole 
spacetimes 

Instead of a restricting attention to a Schwarzschild black hole, the spacetime 
is chosen to be a dirty black hole: 

,2 _ o^^»,\\ / 1 ^(^) \ , 1 J^2 , ^2 102 



ds' = - exp(-2$(r)) 1 - ^ dr + —dr' + r^dn\ (2.60) 

\ r J 1 — b/r 

This black hole spacetime is manifestly spherically symmetric and static, 
however it is not required that the exterior region is in vacuum. A quick 
computation shows that the energy density is given by 

Sttt-^ dr 

This metric is capable of describing the special cases often discussed in the 
literature, such as Reissner-Nordstrom and Schwarzschild-de Sitter, as well 
as more generic black holes. 



The restriction to spherical symmetry is one that is used to make the calcu- 
lations more tractable. The condition that the metric is stationary (i.e. time 
independent) is a reasonable one as it is difficult to define QNM in a generic 
time dependent background. The additional restriction that the spacetime 
is also static is not necessary but simplifies calculations. In particular, we 
know that grr and Qu must go to infinity and zero respectively at the same 
locations to preserve the signature of the metric. As the spacetime is station- 
ary, the place where gu vanishes is an infinite redshift surface as seen from 
infinity. Thus the spacetime ()2.60|) has a horizon where 

b{r) = r. (2.62) 

Note that while gu and 1/grr must share the same zeros, they may be different 
in other locations. The function $(r) allows for this (TH]. In the special cases 
of Schwarzschild, Ressnier-Nordstrom and Schwarzschild-de Sitter $ = 
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{gtt = l/grr), but this is not required by the conditions that the spacetime 
is spherically symmetric and static. So that gu and l/grr still preserve the 
same zeros $ will be assumed to be finite at the location of any horizons. 

2.4.1 The generalised tortoise coordinate 

When defining wave equations, it is useful to make the t-r piece of the metric 
conformally flat: 

ds' = exp(-2$(r)) ( 1 - -] ( -dt' + ^^^i^^l^^dr^^ + (2.63) 
\ '"/V (1 — o/r)^ / 

= exp(-2$(r))(l - b/r) (-dt^ + drl) + r^dQ"^. (2.64) 

The last equation introduces the tortoise coordinate r^, up to an arbitrary 
constant by 

drj, exp($(r)) 



dr 1 — b/r 



(2.65) 



It is intuitively clear that the "fiattening" of the space to produce will 
mean will approach infinity as a horizon is approached. To prove this, 
expand b{r) about a horizon Tq: 



d6 

^(^) = ^0 + d^: 



, 1 d% 
(---o) + ^^ 

ro 



(r - ro)' + . . . (2.66) 



where b{ro) = ro has been used. This implies that 
dr^ rexp($(r)) 



dr r — b 

r exp($(r)) 



(r-ro)(l-6') 
exp($(ro)) f r 



(2.67) 

dr (2.68) 
dr (2.69) 



1 — 6' J r — tq 
= (r - ro + ro In |r - ro|) . (2.70) 

1 — 

The special case b' = 1, corresponding to an extremal horizon, has been 
excluded and is discussed further in jTHI. As r — > ro the generalised tortoise 
coordinate tends to infinity. 
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Near a horizon the and r are related to the surface gravity of the horizon. 
For a spherically symmetric static black hole the surface gravity is given by 



1 

K = — 

2 



dr 



(2.71) 



The surface gravity will be discussed more carefully in chapter ^ In a dirty 
black hole spacetime we have 

exp(-2$(ro)) 



2ro 

The equation ()2.70p then becomes 

1 / r — Tq 



:i - b') 



2Kexp(<l>) \ Tq 

In |r — ro|. 



+ In |r — ro 



2Kexp(<l>) 
Inverting this relation we have 

r — ro ~ exp(±2fi;exp($(ro))r^ 
as the horizon is approached. 



(2.72) 



(2.73) 
(2.74) 

(2.75) 



The potentials that describe dirty black hole spacetimes generally have 
polynomial falloff as a horizon is approached. As seen from ()2.75j) this implies 
that the falloff in will be exponential. If the region of interest lies between 
the outermost horizon and infinity, however, then polynomial decay occurs 
for both r and as infinity is approached.'' 

The next task is to find the wave equations for different types of fields. In 
slightly more detail, the function that plays the role of the potential relies 
not only on the spacetime, but on the type of field under consideration and 
its spin. In practice, for the static case we are perturbing the geometry but 
requiring that the linearised Ricci tensor is unchanged. Thus the QNMs of a 

^As the outermost horizon is approached the decay is stiU exponential in r^,. 
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field are determined by how that field contributes to the stress energy tensor. 
However, in the literature there are comments about the potential for spin j. 
More appropriately, these are the gravitational, electromagnetic and scalar 
potentials. In vacuum spacetimes the difference between spin j and a specific 
field is not great, as the field is not in the spacetime and at linear order basic 
fields look the same. We are also spoilt in vacuum spacetimes by the fact 
that different polarisations of the same field are isospectral, so while different 
polarisations are technically associated with different potentials the spectra 
are the same and people sloppily talk about "the" potential for a spin j field. 
In the non-vacuum case different polarisations are not necessarily isospec- 
tral, for example Victor Cardoso showed in his PhD thesis ^H] that the two 
gravitational wave polarisations are not isospectral. The motto is that the 
effective potential in a generic spacetime needs to consider the polarisation 
and the specifics of the field. 

Although I disagree with the interpretation of the "potential of spin j" , it 
is a commonly used potential in the literature. For this reason, I repeat the 
standard formula below. The spin j effective potential in a Schwarzschild 
spacetime is claimed to be |1H| 



2.4.2 Scalar wave equations 

Even in the case of massless scalar waves propagating on a background there 
is some cause for ambiguity. The simplest generalisation of the wave equation 
for scalar particles is 



□0 - + C - - U = ^'^*VaVfe0 - m^(P +{(--] R(P = 0. (2.77) 



Here m is the mass of the excitations of the field and ( is the conformal 
coupling. If C, is zero and m = then ()2.77|) is conformally invariant. The 
case C = 1/6 eliminates the final term and the resulting equation is called 




(2.76) 
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the minimally coupled Klein-Gordon equation. Note that the mass automat- 
ically breaks conformal invariance, so for massive fields ( = 1/6 is almost 
universally used. However, the following is no more complicated by keeping 
( general. 

Let us rewrite the term involving the Laplacian. We have 



where (f)^ is a one-form. By exploiting the metric connexion we have 



□0 = ^7'"'V«Vfc0 



(2.78) 
(2.79) 



(2.80) 



(2.81) 



We also exploit the well known formula [TH] 




(2.82) 



to obtain 




(2.83) 
(2.84) 





(2.85) 



(2.86) 



We now have to apply ()2.86|) to the dirty black hole spacetime: 



\/\g\ = exp(— $(r))r^ sin 9. 



(2.87) 
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The diagonal metric simplifies terms immensely: 



E [9'"'V\9\^.a) ^ (2.88) 



e-*('-V2sin^^ {g''dU + 9''^dli 
+ e-^^'^^r'^g^^de sin Odect) + sin Odr (^'^'^e'^^'^V^a.^) (2.89) 



9\ { 9''dU + -t^AA + ^-2^8 sin Ode 



+ sin (^"'^e-*("V29,0) (2.90) 
: ^/\i\g''^^^ + sin^a, (/'•e-^^'^V^a,. 



+ ^ ( -^S.sin^a, + -^aj ) 0. (2.91) 
r"' V sm 6 sm ^ 



The final term in parenthesis is the Laplacian on the unit sphere! Let us 
denote 

= -^de sin 9dg + (2.92) 
Our four dimensional Laplacian then takes the form 

^ ab 



□0=^((?»Vl^7l0,O (2.93) 
9\ 



gttQ2^ + '-^dr (^7-e-*Wr2a,0) + (2.94) 



Mr) 

We now separate this equation into angular and non-angular parts by making 
the ansatz 

0(r,^,(^,t) = rV(r,t)r(^,<^) (2.95) 
The Klein-Gordon equation ()2.77j) separates: 



^2 tt *(r-) V72y 

3 .qIp + "L^Q^ (^-e-*(^V2a,rV) - + (C + 1/6)^?] - ' 



r 



(2.96) 

It is worthwhile to note that m is independent of Q and y?, as it is a scalar 
in a spherically symmetric geometry. In the case of a conformal field we are 
using = i?/6 which may depend on r. 
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To find regular solutions of the angular equation we require that 

V^F = + l)Y. (2.97) 

Hence the angular pieces are the Legendre polynomials. The non-angular 
pieces are less trivial: 

- [m^ + (C + 1/6)^] yp) = 0. (2.98) 



The equation ()2.98|) differs from a simple wave equation in the appearance 
of first derivative terms. We may simplify the term with radial derivatives 
in ()2.98|1 by transforming to the generalised tortoise coordinate ()2.()5|1 : 

d^grr^-Hr)r'^g^ = dy{l - b/ r)e-'^^'^ Or (2.99) 

= ^rr'^i- = dy^. (2.100) 

dr^ or OTi, 

To eliminate the first derivative, define the function P{r,t) via 

r'^p{r,t) =r'P(r,t). (2.101) 

The relevant term is then 

a,/^e-*(''V25,(rV) = a,/^e-*("V2a,(r^P) (2.102) 

d 

= dy—r'P (2.103) 

= a, ( r2+-^ + sr^+i^P^ (2.104) 
V dr^ dr^ J 

- dr^ I ^i+s dP ^ 

dr I dr^, dr^ dri 
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Expanding out the derivative and collecting terms gives 

,dP 
dr 



9,/'-e-*('^V29,(rV) = 2(1 + sy+'^ + r 



dr 

+ s(s + l)r'—P + sr'+ , , , 
dr^, dr \ dr^ 



dr drl 
^d f dr 



P. 



(2.106) 



The unwanted first derivative is eliminated by the choice s = — 1. The radial 
term simplifies dramatically: 



dr^ d^P d / dr 



dr 9r2 dr \ dr 



P 



(2.107) 



2 +(<l''(l-6/r) + (6/r)')e-*Mp. 



(1 — h/r) drl 

(2.109) 

In terms of the tortoise coordinate and the Ansatzsdnderung ()2.10H) . equa- 
tion ()2.98|) becomes 

'rexp(2$(r)) d'^P 



2<I.(r) 1 

.dlr-'P + 4 



1 — 6/r 



r^ \ (1 — h/r) dr 



2 +(<l''(l-6/r) + (6/r)')P 



- [m' + (C + l/6)i?] (r-ip) - ^i^-tll(r-ip) = 0. (2.110) 

By multiplying out by the coefficient of the time derivatives we get the equa- 
tion for a scalar field in a dirty black hole in the form of a wave equation: 

Q2p Q2p 

dt"^ drl 
The "potential" V{r) is given by 



+ 1/(r)P = 0. 



V{r) 



r^ D 



$'(1 -6/r) + (6/r)' 



(2.111) 

{X-hjr) e-2*W. 

(2.112) 



Notice that the differential equation is expressed in terms of r^, while the 
potential is expressed as a simple function of r. The problem of finding QNMs 
of scalar fields reduces to finding the purely outgoing harmonic solutions to 
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2.4.3 Gravitational perturbations — a wave equation 

The gravitational case is somewhat more comphcated. By virtue of spherical 
symmetry we can take any perturbation to be axisymmetric (see §4.1 of jl8j). 
Let us take the metric in the general axisymmetric form: 

ds^ = - exp(2[z/°(r) + 6t{r, t)])dt^ + exp(2[/i°(r) + 5,(r, t)])dr^ 
+ exp(2[^°(r) + 5,(r,t)])d^2 ^exp(2[/ij(r,^) + 5^(r,^,t)]) 
X (d<^ - g,(r, e, t)dr - g,(r, 6, t) - qt{r, 6, t)f . (2.113) 

The functions i/°(r), /i°(r), /ig(r) and ^^{r,9) are arbitrary at this point 
and represent the background spacetime. A dirty black hole spacetime is 
recovered when we take 

exp(2z/°(r)) = exp(-2$(r))(l - b/r) (2.114a) 

exp(2/^°(r)) = (2.114b) 
I — o/r 

exp(2/i^(r)) = (2.114c) 

exp(2/ij(r, 9)) = sin^ 9. (2.114d) 

The functions 6i and qi are then considered to be perturbations on the back- 
ground spacetime. Gravitational wave equations are found by insisting that 
the perturbations be "massless" ; more precisely we insist that the variation 
in the linearised Ricci tensor vanishes. The problem decouples into two cases: 

• 6i = 0,qi^ 0: 

These terms introduce cross products g^i into the metric at linear order, 
and thus contribute to frame-dragging effects. Reversing the sign of (p 
will reverse the direction of these effects, hence they are called axial 
perturbations. 

• Si^O,qi = 0: 

These terms induce no frame dragging and are unaffected by reversal 
of (p. These perturbations are referred to as polar. 
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The terminology polar and axial terminology is misleading as, in a spheri- 
cally symmetric system, the names suggest that by a rotation one could be 
transformed into the other. This is certainly not the case! 

In order to obtain the perturbations we define 

z/(r,t) = z/°(r) + 5t(r,t) (2.115) 
/ii(r, e, t) = /i°(r, 9) + 5,(r, 9, t). (2.116) 

We could take the metric, calculate the Ricci tensor directly and then discard 
higher order terms. Instead I shall find the Ricci tensor in an orthonormal 
frame and then discard higher order terms. The metric ()2.11H|1 suggests the 
following orthonormal covariant basis: 

6(0), = (-e^ 0,0,0). (2.117) 

e(l). = (0,e^^O,0), (2.118) 

e(2), = (0,0,e'^«,0), (2.119) 

e(3), = (-?*, -gr,-ge,l),e'^^ (2.120) 

The corresponding contravariant basis vectors are given by 



e 



(0) 



(e-^0,0,g^e-'^)* (2.121) 

eji) = {Q,e-^\Q,qre-^^Y (2.122) 

ej2) = (0,0,e-^«,g,e-^«)* (2.123) 

ejg) = (0,0,0,6-'^^)'. (2.124) 

It is easy enough to check that these do form an orthonormal basis. The 
definition of a general tetrad is ([ZHj chapter 3): 

where t]^,^ is a constant matrix and rj'^'' is its matrix inverse. In the case of an 
orthonormal basis 1] is the special relativity metric. To construct the Ricci 
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tensor we require the Ricci rotation coefficients defined by 

i^fiua = e(^)*Vae(,,)^ (2.125) 

Notice that the first two indices are simply labels, and the 24 Ricci rotation 
coefficients are really 6 one- forms Ufj^i,. The notation here differs from Wald 
|79j . but his notation tends to hide the geometric nature of the rotation 
coefficients. It should be noted that uj^i^a is antisymmetric in the first two 
indices. We rewrite the coefficients in terms of directions of the tetrad, rather 
than t, r, 9 and ip: 

^ijiva= ^^i,vae''"\. (2.126) 

We can then evaluate the spin coefficients by combinations of Lie brackets: 
uJ^^ua = ^ (e(a)fc [e(/^),eH] '' + 6(^)6 [e(„),e(^)]* + e(^)^ [e(^), e(„)]^) . (2.127) 

In a torsion free theory, the Lie bracket is constructed using partial deriva- 
tives and the tetrad basis vectors. For example 

[e(o),e(i)]'' = e^o)%e1,^ - e^i^%el,) (2.128) 

= {e-'dt + qte-'d^) (O, e-^''^^'*), 0, g,(r, ^, t)e-^'-('-'*))" 

- (e-'^^ dr + qre-'^d^) (e'^^^'^'*), 0, 0, gt(r, ^, t)e-^('-'*))" 

(2.129) 

= exp{-fir - I') {i',T , -/i,t , , 2q\r,t\ + qti^,ry ■ (2.130) 

In evaluating the Lie brackets, it is worthwhile noting that the tetrad vectors 
are independent of ip, due to the fact that the spacetime is axisymmetric. 
Any time derivative is automatically at least first order, as the background 
is independent of t. The other Lie brackets are given in table ITTl 

The procedure of getting the spin coefficients is now straightforward, if 
somewhat tedious. The result of doing this calculation is summarised in 
table 12.21 We can now use 

Rpafiu = '2e[{p)\"'daUJp„\a] — 2?/"^ {uJl3fi[p^\au\a] + ^^f3[ap]^fMiya) (2.131) 



97 



[^(0), 

L V^/ 


6(1)] 


= exp(- 


-n - 




[^(0), 

L 


6(2)] 


= exp(- 


-V - 


-fie) (0,0, -fie,t , 2g[6),t] - qelJ^e,t) 




6(3)] 


= exp(- 




-u) (0,0,0, 




6(2)] 


= cxp(- 




- fig) (O ,0 , —fle,r , 2g[6l,r] " /We,r?e) 




6(3)] 


= cxp(- 


-fir 


- (0 , , , -fl^^r) 


[e(2), 


6(3)] 


= exp(' 




- /i^) (0 , , , -fi^,e) 



Table 2.1: The six independent Lie derivatives of the tetrad basis 



or in less mathematical notation 

Rpap.0 — e[(p)\°'daUJ^,y\a] — V°''^ {^^Pfi[pi^\au\a] + ^0[ap]^pua) 

— (everything on the first line, but a ^ p). 

The orthonormal Ricci tensor is then obtained by contracting over the first 
and third indices. 



Notice that up until this point everything has been exact; no use of the 
fact that qi and 5i are of first order has been made. Let us denote the Ricci 
tensor by Then the change in the Ricci tensor must be at least second 
order: 

5R,i^R,i-R%^0{{q,,5,f) (2.132) 
The linearised perturbations fall into two separate categories: 

• Three are axial: these terms vanish in the background and all the 
corrections are linear or higher in qj_. Hence any term containing 5i is 
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=-i/,rexp(-//) 




'^oii 


=-//,rexp(-i/) 




^013 


=exp(//^ - 1/ - llr)qir,t] 




^022 


=-//e,texp(-i/) 




^023 


=exp(/i^ - i/ - i^e)q[e,t] 




<^03i 




<^013 




=exp(//^ - i/ - 


''^023 


•^033 


=-/^v,texp(-i/) 




'^122 


=-/^e,rexp(-//^) 




<^123 


^ex.p{ii^ - He - lJ'r)q[e,r] 




'^136 


=- exp(//^ - 1/ - 


= -<^03i 


'^132 


=- exp(//^ - //r - /^e)?[r,e] =+^^123 


'^133 


=-//<^,rexp(-//r) 




''^236 


=exp(//^ - i/ - /X0)g[t,e] = 


''^032 


<^23i 


=exp(//<^ - /^r - /^e)?[r-,e]= 


'-<^123 


'*^233 


=-//<^,eexp(-//0) 





Table 2.2: Spin coefficients for the given tetrad basis. They have been 
grouped by the ffist two parameters (which are antisymmetric) and the last 
index is the direction of the one-form. All other components are zero. 
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of second order. 



5i?o3 = exp(-2/x° - /i^) I exp(-/i°) {q[t,r] exp(3/xj - i/° - /i^ + /i^)) 

+ exp(-//° - 1/°) {q[t,e] exp(3/xj)) ^ | (2.133a) 

(5i?i3 = exp(-2/x° - - 1/°) {qir,e] exp(3/i° - - + u^)) 

+ exp(/xj - 2j/° - /x°)g[t,,],t (2.133b) 

+ exp(-2z/0 + /.J - /i°)g[t,,],t (2.133c) 



The remaining seven are axial. These terms do not vanish in the back- 
ground spacetime. The Qi only enter at second order, so are irrelevant 
for calculating the linearised perturbations. The 5i do enter at first 
order. At linear order the Qi and 5i are decoupled, in the sense that 
the axial perturbations rely only on and the polar perturbations rely 
only on S^. As the gj's can be ignored for the polar perturbations, this 
problem is spherically symmetric. This symmetry reduces the number 
of independent equations from seven to five. 

In writing down the axial equations the values of //^ and //° for the 
background spacetime have been substituted: 



SR^^ = ?^^L^^^ [^S, + 6^),e + {S^ - So) cot 9] 



dt 



5i?6i = exp(-z.» -/.;!) - 



-Sr - {So + S^)^r 



+ [K--]{So + s^) 



(2.134a) 



(2.134b) 



SR 



exp(-//°) 



12 



{St + (5^),r6» + {S<p - So),r cot 6 



J^,r+-\ Sr.e +{{l'j 



(2.134c) 
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'1 

r 



+ {25^ + 5t-5e)fiCoie\ +e 



(2.134d) 



5R 



33 



-e 5^^tt + ^ + + (5^ - 5e + 25^ + 5t)fi cot 1 



r 



- [St - Sr + 5e + 



2 1 

(2.134e) 



Axial perturbations 

Let us now restrict our attention to tlie axial perturbations. To simplify the 
equations, we define the variable Q by 



Q(r, e, t) = q[r,e] exp(3/ij - - /ig + z/° 



(2.135) 



Replacing the and z/° with their values from the background spacetime 

h{r) 



Q = r'sin^(^) 1 - 



~^(r) 



(2.136) 



We can now rewrite the equations governing the perturbations ()2.133b|) and 
fl2.133cj) in the following form: 

r^sin^(6') 



Q 



e-*(l - h/r) 
Qfi = -r^sin^(6')e*g[i,r],t 



(2.137) 
(2.138) 



We now assume that the perturbations are harmonic i.e. the time dependence 
is of the form exp{iujt). The perturbations are then governed by 
2e-/^^(l - b/r) 



r2 sin'^(^) 



Q^r = 2iu;g[t,e] = iujqt,e + 



—lujqt^r — ui qr- 



(2.139) 
(2.140) 
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The qt term can be eliminated by virtue of the fact that partial derivatives 
commute: 



2e-M^(l - b/r) 



sin^(^) 



Q,r] +{ -T^T77^Q,e ) = ^\Qe,r - qr,e) (2.141) 



sin^(6') 



Q. 



r2 sin3(0)(l - 6/r)e-* 
This leads to an equation in separable form: 



(2.142) 



Q,, +sin^(0) 



sin3(^) 



(1 - h/r) 



We take as our ansatz: 



Q = 0. 
(2.143) 



Q(r, t) = rR{r)^{e) exp(icjt). 



(2.144) 



which separates out ()2.143|1 : 

'e-*(l-6/r) 



4 <i> 

re 



ri?) + r e 



:i - h/r 



-rR - j^rR = 



(2.145) 
(2.146) 



where the separation constant has been denoted 7^. The angular equation 
p.l46|l is solved in terms of Legendre polynomials and Q^: 



^{9) = sin^ 9 



\/9+472-l 



\/9+47^ 



(cos^) + Sg2 ' (cos( 



(2.147) 



To ensure the solutions are real, we restrict the values of i to integers. This 
requires that 7 is related to i by 



=e + i-2 = (£ + 2)(£-l). 



(2.148) 
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The aim is to rewrite the radial equation ()2.145|) in terms of a wave equa- 
tion over the real line. Changing variables to the tortoise coordinate (cf 



dr^ 1 



dr exp(-$)(l - 6/r) 
the radial equation becomes: 



(2.149) 



R ^ dr d r dr\R o f V ^ | -^i? cxp( 'I') - 

'^*^* dr^ dr \dr^J r \dr^J dr^ 

(2.150) 



This can be written in "Schrodinger form" 

R^r.r. + u^R = V{r)R (2. 151) 

with the "potential" V{r) defined as 

^/ N _dr/ d f dr\l ^^f dr\ I ^-i"^ ^ ^ \ 
dr^ \ dr \dr^/ r \dr^^) r^ J 

l_i!)('i.('l_^),^k,!(i±i)_») ,2,53) 

where in the second term we have substituted back ()2.65|) and written 7^ in 
terms of the angular momentum ^. 



So far we have not used the 5Rq^ equation. This equation can be used 
to complete the solution - that is, solved to find qt and g(r,e)- However the 
form of qt is not of direct interest to us and it is assumed that all of these 
equations can be simultaneously satisfied. 



Polar perturbations 

The equations governing polar perturbations ()2.134|) can be simplified by 
assuming an exp(zco't) time dependence for all 5j. In this case we can simplify 
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mm- 

= 
= 

= 



{Sr + S^)^0 + {S^ - Se) cot 6 

{St + 5^),re + {S^ - 6e)^r cot 9 - (^U^r + Srfi + ^ 



(2.154a) 
(2.154b) 



- 1 

r 



(2.154c) 



-UJ' 



exp(-2z/°) (5^ + 5e)-— [{5t + (5^),ee + 25^ 



(2(5,^ + (5t - 5e),e cot 



■(45,1/° -25t,,) 



1 \ 2 

- + (5e + 5^),r + 



(2.154d) 



= u'e ^"6^ + — {d^^ee + 25^ + (5r - 5e + 25^ + 6t)fi cot i 



z/° - /i° + - 

r 



2 1 



(2.154e) 



In the vacuum (Schwarzschild) case, the above equations are considerably 
simphfied by = — z/°. In this case, the above equations are identical to 
equations (31)-(35) in chapter 4 of Chandrasekhar [THj. Let us introduce the 
Ansatz (due to Friedman [21]) to separate the equations: 



St{r,e) = N{r)P\cos9) 
6r{r,e) = L{r)P\cose) 

6e{r,e) = T{r)P\cose) + A{r 
(5^(r, 6) = T{r)P\cos 9) + A{r) 



d^P^(cosi 

' d^2 

dP^(cose 



cot 6*. 



(2.155) 
(2.156) 

(2.157) 
(2.158) 



Here P^(cos^) is a Legendre polynomial, which will be abbreviated to P^ in 
future. The first two equations are standard separations, whereas the last 
pair have been chosen separate out the radial and angular parts. 
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The 5Rq2 ()2.154a|) equation shows that not all the radial functions are 
independent: 



= 5rfi + + (5^ - 8e)fi cot 6 
= {L{r)P'),e + {T{r)-A{r))P^, 
= {T + L- A)Pl. 



(2.159) 
(2.160) 
(2.161) 



Following Chandrasekhar we choose to eliminate T(r) in favour of L{r) and 
A{r). 



Turning attention now to (j2.154bj) we have 
2 



= -LP' - [2{A -L)- i{i + l)A]^rP 

+ (^z/J, - [2(A - L) - i{i + 1)A]P' 

i-^'r + l) [2L+{i + 2){i-l)A] + h. 



(2.162) 
(2.163) 



Equation ()2.154c|l becomes: 

= {[N ~ L]P^ + A[P' + P^g cot 0]) . + A^r [Pfi cot e - P^ss) cot 9 



re 

-[v.r + -1 LP' + (v r - - 1 NP', 



^2 



(2.164) 



[A^ - LlrPfi + Ar[Pfi{^ + cot' 6) + (P;, cot ^),, - P^.g COt I 



yr + -\ LP 



r 



(2.165) 



— + z/^-- iV- — + + + 
dr r / \ dr r 



sin'^ 



+ (cot0),0 



,9 

(2.166) 



which implies that 



— + I/,, - - 1 AT 
dr r 



+ Z^,r + - 1 

dr r 



(2.167) 
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Up to this point the equations derived agree exactly with Chandrasekhar's 
treatment of the Schwarzschild black hole. This is not true of the final two 
equations, where the dirt surrounding the black hole makes a non-trivial 
difference. Taking ()2.154d|) : 







-LU' 



exp(-2z/°) (2T - £{£ + 1)A) P' - - [NiP% + P^e cot ( 



- ii + 2){i-l)TP^] +6-^"" 



-{2Lv%-N,. 



P'. 



By using the defining differential equation 

Pfi6 + P'eCoie = -(i{l + l)P'. 
for Legendre polynomial we find 

,0 



(2.168) 



(2.169) 



^ (2L + (^ - + 2)A) + [iV£(£ + 1) 



dr 



+ (£ + 2)(£-l)(A-L)] +e 



r \ ' r I r 



0. (2.170) 



The final equation ()2.154e|) is also the most complicated: 



= uj^e-^" 6^ + — {6^,ee + 26e + (5r - Se + 25^ + 5t),e cot ( 



+ 



-)5r 



+ -[5t- 5r + 5e + 



(2.171) 



= u^e-^"' [TP^ + AP^g cot e] + \ {-{£ + 2){£ + 1)TP^ + {N + L)P^g cot 9) 



,.0 



2 
r 



d ■ 

dr 



{TrP^ + A^rPlcote) 



+ I - ^° -^Alp' + \[N-L + 2T- £{l + l)A\^ P' 



(2.172) 
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We now collect terms proportional to and Pj, separately: 
= {^2g-2.0y _ + 2)(£ + 1)T + ^6-"-° (^fil - ul - L+ 



[N -L + 2T + I) A] ^ + e- 



{-V-"A+l(iV + L) + e- 



d 
dr 



2/.S! 



d ■ 

dr 



A^r \PnCOte 



(2.173) 



(A - L) - (£ + 2)(£ + 1)(A - L) 
+ ^^-2^0 _ ^0 - i") L + ^ [iv - 3L - (£ + 2)(£ - l)A]^ 



+ e 



- III + 



+ 



dr 

2/.J! 



(A-L),,}P^ 



+ 



dr 



A,,|p^gCot^. 

(2.174) 



The final expression must vanish for all Q. But the only Q dependence is in 
the linearly independent factors P^ and P^ cot Q. For ()2.174|1 to be satisfied 
for all Q we require 

^2g-2.0(^ _ ^) _ (£ + 2)(£+ 1)(A - L) + ^e-2/^° - vl - L+ 



-2mSJ 



[A^-SL- (£ + 2)(£- l)A]_, + e- 



-2m? 



dr 



{A - L),r = 0, 
(2.175a) 



and 



u'e-^''°A + \iN + L) + e-^''" 



dr 



A,^ = 0. (2.175b) 



Deriving the Zerilli equation 

In the previous section we have taken the equations for polar perturbations, 
although finding the actual equations we need to satisfy gets lost in the 
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working. The equations governing polar perturbations ()2.154a|) - ()2.154e|) and 
rewritten with the variable (3 = {£ + 2){i ~ l)/2: 



-] |L + PA] + -L = 

r r 



dr ~'~ r / I dr ~'~ ~'~ r 



+ 2(3{A - L)l + e 



r 



dr 



- [2^1 + -]L- -N, 
r \ ' r I r 







^^e-''^ (A - L) - 2/3(A - L) + -e- 

r 



[AT - 3L - 2/3 A] , + e 



a; e 



yl + 4(Ar + L) + e- 



-2/.S! 



/i° + 



dr 



+ 



dr 



(2.176a) 
(2.176b) 



(2.176c) 







(2.176d) 



A^r = 0. (2.176e) 



This is five equations in three variables, namely A^, L and A. The first 
three equations ()2.176ap - ()2.176cp form three first order equations in the three 
variables. As such, we can solve for the first derivatives of A^, L and A 
in terms of these same three functions. The simplest way of solving these 
equations is to introduce functions a(r), 6(r) and c(r) such that 



A^,^ = dN + bL + d(3A, 13 



(£ + 2)(£-l) 



(2.177a) 



Using equations ()2.176b|) and ()2.176cp easily yield 



-\L + c(3A 
r 



(3A^r 



a + .°-i)Ar+(2.°-6-i)L+(.? 



1 



(2.177b) 



~c]pA. 



(2.177c) 
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By substituting these expressions into the longer and more messy fl2.176c|) 
we can then find these three functions. The results are 

a = ^^^i^ exp(2/x°) = ^ (2.178a) 
2r r — 

6 = -- - ^ + T^^e^* + (z/°)V + vl (2.178b) 



r r — 6 {r — b) 
1 1 r^u;^ 
r r — 6 (r — b) 



+ ^ + T^^e^* + (^°)'^' (2- 178c) 



where the values for z/° and have been substituted back into the dirty black 
hole functions $(r) and m(r) via ()2.114a|l . The derivative v% is complicated, 
and the expressions are simpler not expanding it. 



The task now is to find a variable that can be expanded into a one dimen- 
sional wave equation. We try the function 



Z{r) = f{r)pA{r) + g{r)[L{r) + pA{r)]. 



(2.179) 



Intuitively the three functions A^, L and A span a three dimensional vector 
space, and the above combination has no "component" . If we can express 
the Zerilli equation as a one-dimensional wave equation then I require to 
have Z^r^,r^, pointing in the same direction, so in particular the N component 
must cancel. The L and A components have been chosen as candidates for 
a one dimensional wave equation as ()2.17(iaj) tells us that at least [L + l3A)^r 
transforms into a component independent of N. We have 



1 ry^ 

— = f\r(3A + ff3A,r + g,r{L + PA) + g{L + /5A), 



f,rPA + fPA,.+ 
lr/3A + f 



+ 



9 ^. 



1 

a 

r 

n 1 

r 



9,r+9[ - 

N + 



[L + PA) - 

r 



2Z/'; - 6 - - 
r 



L 



(2.180) 
(2.181) 

--~^PA 



{L + PA) - ^L. 

r 



(2.182) 
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Before taking the second derivative, we regroup the expression into pieces 
with (L + (3A) 



f 

+ 

f 
+ 



(^l-a- z/° ^ N + (^u^ - 1 - (L + + (a + 2i/°)L 



9,r + 9[ T^,r - 



- a - V, 



(L + pA)-^L + lrPA 



(2.183) 



N+(a + 2i^Z)L 



9,r + 9[ ^. 



1 



(L + /5A). 



(2.184) 



The procedure to obtain the ZeriUi equation is now clear: we differentiate 
this expression again and use the fact that 

e^'^°+^*([2^J + $,]Z, + Z,,). 



(2.185) 
(2.186) 



By ehminating the term in N we will get one differential equation in / and 
g. Another equation is obtained by setting the ratio of the / and g terms to 
be the same as the ratio of the coefficients of (3 A and (L + (3 A) in 2'^^^^^ . This 
is two equations in two unknowns and hence in principle we can determine / 
and g. If the solution turns out to be a linear equation then we trivially have 
a wave equation for polar perturbations. With some regret, I must report 
that I have not done this yet. The difficulties are purely mechanical, for 
example the second derivative is given by 



j,2 ^,rr 



f 



+ 



a-c + 3.°-l)/,+ (2.o,-ai-^ 



N 
3z/0 



— H \ l-v 



(2^° - ^) 9,r + 9, 



{L + PA) 



-9,r 

r 
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c 1-6 
- + } 



/ 



+ 



a — V 

r 



fNr + 



/(a + 2z.°)-^ 
r 



L,r + UP A. 



(2.187) 



We would need to expand the derivatives on the last line and substitute using 
fj2.177a|) - ()2.177c|) . then collect terms and place into ()2.186|) . This has not yet 
been completed. 



The Zerilli equation for a Schwarzschild black hole 

I shall not derive it, but I quote the significant results for the special case of 
a Schwarzschild black hole from Chandrasekhar's book jTHj. In this case the 
one dimensional variable that obeys the wave equation is given by 



The Zerilli equation is given by 



+ LJ^Z = V,Z (2.189) 



where the effective potential is given by 

At some point I would like to find the dirty black hole equivalent and show 
that it reduces when b = 2M and $ = 0. 



2.4.4 Isospectrality and the NP formalism 

Even though the solution has not been obtained for the Zerilli equation, it 
is at least clear how we could obtain it given enough time and patience. 
The method does seem devoid of any sort of intuitive interpretation, other 
than we have realised that there are two decoupled modes. The rest is just 
going through rather tedious algebra. The NP formalism seems to allow 
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for a cleaner separation of variables, and makes various appeals to the fact 
that stationary black hole spacetimes are of Petrov type D to considerably 
simplify matters. Even after finding the principal null directions in static 
dirty black hole spacetimes the spin coefficients did not separate into two 
similar sets of equations. To some degree, one would expect that the dirtiness 
of the spacetimes would complicate matters somewhat, as in the standard 
approach in vacuum spacetimes the NP formalism a natural way of seeing 
that the spectra are isospectral. As we know of examples where the spectra 
are not isospectral (namely AdS) then it is possible that the NP formalism 
is not of great utility. In any case, I have investigated the possibility and 
have found the non-zero $o, and $2 have prevented me from finding any 
obvious separation of variables. 

2.5 Conclusions 

This chapter has hopefully served as a brief introduction to the idea of QNM 
and how they are defined, as well as some of their analytic properties. The 
QNM conjecture has also been outlined, which has served as the motivation 
for looking at the highly damped QNMs. After all, they are not likely to 
contribute in any significant manner to any astrophysically relevant process! 
The method to obtain the wave equation for a QNM is outlined, and the 
equations are found explicitly for axial gravitational fields and for scalar 
fields. In the next chapter methods of probing the asymptotic structure of 
the QNM spectra are discussed. 

It has also been pointed out that finding the polar perturbations is non- 
trivial, and while it has been indicated how one could carry it out in principle 
the task still remains incomplete. But a more serious problem has developed: 
while the QNM conjecture looked promising it was ultimately based on a 
numerical coincidence. As such, it is a promising direction for research and 
investigation, but we should not argue vehemently that it is true until we 
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have solid arguments leading us to the QNM conjecture from well established 
basic principles. If one adopts this philosophy, a recent paper by Domagala 
and Lewandowski has shown the QNM conjecture to be false. Recall that 
we anticipate being able to set the Barbero-Immirzi parameter as follows 

7 = -^^ 0.123637... (2.191) 
2v 27r 

More generally, if we allowed the black hole to be dominated by links with 
spin j we found that to get the correct relationship between area and entropy 
required 

2vrv/jU + 1) 

The choice j = 1 is made by the connection to QNM. If the QNM conjecture 
is not assumed then there is no reason to assume j = 1, any half- integer value 
in ()2.192j) will do. The Lewandowski and Domagala paper shows that in the 
case of j = 1/2 that ()2.192j) does not hold. It is shown that to reproduce the 
area of the horizon you need 70 to be the solution of 

00 ^ 
J2 exp(27r7i/2Vm(m + l)) = - (2.193) 

m=l/2,l,3/2.... 

which leads to 

7^=1/2 = 0.2375329... (2.194) 
whereas the solution predicted by ()2.192|1 is 

^.^y^ = }^ ^ 0.127384 . . . (2.195) 
7rv3 

This is because in considering ()2.192j) we did not take into account various 
constraints, such as the sum of all the spin links intersecting the horizon 
being zero (as we are considering a Schwarzschild black hole). The value 
'jj=i/2 = 0.2375 . . . works over dirty black holes, dilaton black holes, rotating 
and charged black holes, and hence it seems that it is robust. However, a 
similar result will be true for the j = 1 case as well, as the same constraints 
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are still there. Thus 7 will be different from what is predicted in ()2.192j) and 
we will no longer be able to make a correspondence with the QNMs while 
retaining the semi-classical relationship between area and entropy. 

Is this too dim and dark a conclusion? Is it not possible to fix up the QNM 
conjecture in some manner to overcome this difficulty? Not an obvious one 
that I can see. Of course one can try and repair it, but such a task seems 
fruitless. The QNM conjecture had a nice intuitive appeal on the basis of 
particles and some [admittedly odd] use of a correspondence principle. The 
only evidence for it, however, was not a solid argument but an apparent 
numerical coincidence. Now even that has gone. As Feynman remarked 
about proving the mass of the photon is zero, you can only prove that it is 
less massive than some given upper mass limit. Thus, in principle, if the 
photon is truly massless you could spend forever getting better and better 
bounds on its mass. At some point we just accept it is massless until we have 
a reason to believe otherwise. With the QNM conjecture a similar principle 
should hold: until we get so desperate for ideas on how to explain aspects of 
quantum gravity, the QNM conjecture should be discarded. 
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Chapter 3 



Finding the QNM 



In the previous chapter the definition of QNM was given, and the problem 
of finding QNMs was related to finding solutions of the Helmholtz equation 



with the boundary conditions that the wave is outgoing at both the horizon 
and infinity, and that the time dependence of the perturbation is of the form 
exp{—iujt) . These boundary conditions can be rewritten explicitly in terms 
of as: 



if we take '^{uj) < 0. Solutions with '^{uj) positive are not allowed as 
they would correspond to solutions that grow exponentially in time, whereas 
QNMs are the decaying modes of a spacetime.^ It is easily seen that these 
boundary conditions ensure that the solutions are outgoing at both the hori- 
zon and infinity. 

^In the language of ^2. 3. li the exponential solutions correspond to /+ and which are 
the well behaved solutions to the homogeneous equation when either infinity is approached. 



dr^ 



+ ij'R{n) = V{n)R{r,) 



(3.1) 




(3.2) 
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The previous chapter also showed how to find V{r^) for specific fields. 
Together with information of the initial perturbation we can then find how 
various fields are perturbed in the spacetime. The goal of this chapter is to 
outline methods that may be used to calculate the quasinormal frequencies. 
The problem with solving (jH.lj) subject to the boundary conditions are: 

• The generic problem of solving a partial differential equation. 

• While V{r) is a rational function of r, there is no (simple) relation in 
terms oir^,{r). For example, in the Schwarzschild geometry we need the 
Lambert W function. Hence V{r^) is generally not a simple function. 

• The boundary condition is numerically difficult. In general we have a 
superposition of an exponentially growing term and an exponentially 
decaying term as — > cxo at fixed t. The QNM boundary condition 
()H.2|) requires that the exponentially decaying piece vanishes. In a prac- 
tical numerical solution, the exponentially decaying piece would be lost 
in the numerics, as the exponentially growing term would dominate. 

To find the quasinormal modes we look at four different methods: 

• Analytic approximation methods: By replacing the function V{r) by 
an approximate function V(r) we are able to solve the new Helmholtz 
equation with the required boundary condition ()H.2j) for the frequen- 
cies Un analytically. The underlying assumption is that if V is well 
approximated by V, then Un will be well approximated by Un- 

• Series expansion methods: By making a series Antsatz for the form of 
the solution a recurrence relation or continued fraction method can be 
found for the QNM. Also known as the Leaver method jSH], and also 
used by Motl 

• WKB or Monodromy methods: The Helmholtz equation (jH.lj) is for- 
mally analogous to the Schrodinger equation. By treating the QNM 
problem as a pure scattering problem, an analogue of the standard 
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Bohr-Sommerfeld rule is found. This technique was used by Motl and 
Neitzke in their recent paper [HO] . 

• Born series methods: By treating the QNM problem as a scattering 
problem, it is hoped that the Born series would give accurate results; 
especially for highly damped modes. 

With the exception of Leaver's method, all of the above methods make an 
approximation at the level of the differential equation, rather than finding an 
exact solution to the exact differential equation. The above list is only a sur- 
vey of possible methods and is by no means exhaustive. In the review article 
by Kokkotas and Schmidt ^4^, the last section outlines several techniques 
for finding QNMs. In particular, they mention techniques such as resonance 
methods, direct integration methods and variation principles. For details, 
see ||4Sj and references therein. 

It is hoped that we could approximate the QNMs of any spherically sym- 
metric static black hole. This would include the special cases of Schwarzschild, 
Schwarzschild-de Sitter, Reissner-Nordstrom, as well as the more generic 
cases where the black hole was surrounded by spherically symmetric matter 
fields.^ Some methods (notably Leaver's) do not require spherical symme- 
try and can be applied even more generally, but the others apply to one 
dimensional scattering processes and it is difficult to see how to generalise 
them. 

Even within the restrictions of a spherically symmetric black hole things are 

still unclear. If we allow arbitrary matter fields, then there is a possibility of 

a cosmological horizon or even a "multi-horizon" structure. It is argued that 

QNMs are a property of the spacetime between two horizons, as this is the 

^The reader is reminded that while axisymmetric matter fields are more astrophysi- 
cally relevant, we are not performing astrophysically relevant measurements here. Rather, 
the goal is to look at highly damped modes and investigate the robustness of the QNM 
conjecture, as highly damped QNM are not of astrophysical relevance in any case. 



118 



only region accessible to a given observer, under the additional assumption 
that neither of the horizons are extremal. It is tacitly assumed that the 
horizons are boundaries of an inward trapped and outward trapped region 
respectively, so that we may apply outward going boundary conditions at 
each horizon. By introducing a generalised tortoise coordinate (as was 
done in Chapter |2)) that maps the region between the horizons to the entire 
real line, the boundary conditions of the differential equation ()3.ip remain 
unchanged. 

Much more subtle is the possibility that the space becomes asymptotically 
anti-de Sitter (AdS), in which case the space is closed (in the absence of 
topological obstructions). Thus there is no hope that the second bound- 
ary is an outward trapped region and the above prescription fails. A QNM 
prescription is still needed, as energy in the exterior region can propagate 
through the horizon and become lost, making the use of normal modes inap- 
plicable. To define the QNMs precisely we need to fix a boundary condition 
at the outer boundary of the spacetime. Conventionally this boundary con- 
dition is taken to be the vanishing of fields at the outer boundary that are 
damped in time. Using this convention the QNMs still correspond to poles 
in the Green's function, as there is no incoming radiation. A different con- 
vention (not adopted here) is that the outer boundary is perfectly reflecting. 
For a discussion on boundary conditions defining QNMs in AdS spacetime 
see §3.3.1 and references therein of Cardoso's thesis [T3] and the paper by 
Horowitz jH]. 

Finally the reader is warned that the QNMs found in this thesis are ob- 
tained by assuming that the topology is the simplest one compatible with the 
local geometry. While topological effects may make a significant difference 
to the QNM structure, I have assumed throughout that the local geometry 
suffices to specify the topology. 
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3.1 Born series method 



The Bom series method is the least rehable of all the different ways of finding 
QNMs and yields very little information. However, it is very simple to apply 
to wide classes of black holes. The Born approximation is described in detail 
by many books on quantum mechanics, for example [221 IHl] • A very brief 
discussion is given here. 



The Born approximation splits the function -R(rj,) into a series 

R{r,) = + R^'\r,) + R^'\r,) + ■■■ (3-3) 

As we are interested in the limit of highly damped modes, we are really 
looking at the modes for which |u;„| ^ sup provided that sup |V^(r)| 

is finite. In this limit, the potential can be ignored in the lowest order (see 
()2.1|1 to recall the definition of Un)- Hence the lowest order solution R^'^^ is 
a plane wave solution. From R^'^^ , the n*'^ Born approximation can be found 
via 



+ cj2/?W = V(r^)i?('^-i). (3.4) 



By finding the Green's function G{r) for the Helmholtz equation we convert 
this into 

POO 

/2(")(rJ= / G(r,-r;)r(r;)i?("-i)dr;. (3.5) 



But for the one dimensional Helmholtz equation the Green's function is well 
known 

G'(r^) = 2 exp(ico'r^) (3.6) 
and hence the Born series becomes 

/oo 
exp{-iujn)V{n)R^''-^\Q dn. (3.7) 
-oo 



120 



3.1.1 Finding the QNM 

We start by utilising the first order Born approximation 

/oo 
exp{—iujfi,)V{fi,) exp(— iturv) d?\. 
-oo 

(3.8) 

The transmitted part is then in the first term, while the second term describes 
particles that are scattered back. We want to concentrate on when the back- 
scattering blows up, so that we can have a case with no incoming particles 
and only outgoing ones - the boundary conditions of a QNM! 

To do this, let us take the second part and label it a(a;). This is essentially 
a Fourier transform of the scattered part. Then we have 



a[uj) oc 



/oo 
exp{-2iuf,)V{f,)df,. (3.9) 
-oo 



Undergo a change of variables back to r. This means that the lower limit is 
Tb, the black hole horizon and the upper limit of integration is Vc, the cos- 
mological horizon. If a cosmological horizon does not exist then we formally 
take Tc to infinity. 

r" ~ dr* 

a(a;) oc / exp{—2iur-),{r))V{r)-r^dr (3.10) 

~ / exp(— 2ia;f,t(r))y(r) 



1 1 

+ TT—, T + a 



2Kb{r - n) 2Kc{r - Tc) 



dr. 
(3.11) 

Here we have approximated 

dr. 11 



dr 2Kb{r - n) 2Kc{r - Tc 



+ «, (3.12) 



where ni, and the surface gravities of the black hole and cosmological 

horizons respectively. The fact that the expression takes this form near a 
horizon follows from earlier comments about the exponential decay of r^^ 
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as a horizon is approached.^ If a cosmological horizon does not exist, the 
requirement of taking Vc to infinity aheady kills off the term that contains 
it. Here we have assumed that the horizon is non-extremal so that the sum 
of the two terms is a good approximation to the derivative of the tortoise 
coordinate, at least near the two horizons. The other large assumption is 
that it is the near horizon geometry that is important. 



Let us assume that r^, <^ Tc. Then the poles that come about from the 
near black hole horizon geometry Uf, are those for which 

1 1 



a{ujb) oc / exp{-2iuJbri,{r))V{r) 



+ 



exp{—2iLJbri,{r))V{r) 



2Kb{r - Tb) 2Kc{r - 



+ a I dr. 



+ a 



2Kb{r - n) 



But near the black hole horizon we also have 



a{r - n) + — In 

2Kb 



r — Tb 

Tb 



+ 0{[r-rbf 



dr 
(3.13) 
(3.14) 

(3.15) 



and hence 

exp(— 2iti;r^) ~ exp(— 2ici;Q;(r — r^)) 
~ exY>{+2iujarb))e' 



r — Tb 

Tb 



r 

1 

rb 



-iuj I K-^ 



+ (exp(2icu[r-rf,]Vft;6)) 
(3.16) 

(3.17) 



where z = —2iuja{r — r^).^ Placing all of this back into the integral ()3.14|) 
we obtain (neglecting factors of proportionality) and making z the variable 
of integration: 

iz 







2uja 



+ n 



(2^ + ' 



(3.18) 



■^Or, as a horizon is a place where 2m(r) — r, we have a pole in the defining expression 
for the tortoise coordinate. If the horizon is non-extremal by definition this pole is simple 
and the result follows. 

*This corrects an error in our paper, where we claimed that z = —2iujar. 
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We can then expand V{rh + iz/2uja) as a power series in z: 



2uja 

Hence we have 



a U<j oc > a. 



oc 



oc 

71 



«n e-^^-'"/'^^-^+" + dz (3.20) 

/■oo 

V&n / 6-"^-*"/"''-^+" dz (3.21) 
/■oo 

V6„ / e-"z-^'^/^''-^+" d^ (3.22) 
■>'o 



oc 

n 



5^6„r(-^+n). (3.23) 



Now the Gamma function blows up at all negative integers. The argument 
of one of the Gamma functions will be a negative integer provided that 

uo = —nbini, (3.24) 

for some integer m. Note that my sign convention is opposite to the one 
used in our paper (appendix Ej). Similarly, by looking at well separated 
cosmological and black hole horizons we find analogously 

uo = —KcfT^i. (3.25) 

In the case where the horizons are close, or more precisely (r^ — rh)/rh is of 
order one or lower, then treating the two horizons as individual supports for 
the first Born approximation fails. See appendix O for details on this case. 



3.2 Analytical approximation schemes 

Analytical approximation schemes are based around the premise that a "small" 
change in the potential will result in a small change in the QNM frequencies. 
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It is not clear that this is always the case, or how one should meaningfully 
quantify "small" . One candidate for a measure of smallness is the integrated 
difference in the absolute value of the actual potential V{r^) and the analytic 
approximation V(r^) 

/oo 
\V{n)-V{n)\dr. (3.26) 
■oo 

One possible way of addressing this question is by looking at generic one- 
dimensional potentials and bounding the reflection and transmission coeffi- 
cients, as is done by Visser [77j and Froman and Froman j28j . 



So far analytic approximation schemes have been applied to four cases 
in the literature. The first paper by Ferrari and Mashhoon jEH] covered 
the Schwarzschild, Reissner-Nordstrom and slowly rotating Kerr black holes. 
Much later the Schwarzschild-de Sitter black hole was investigated by Suneeta 
|73j . If one considers only the near horizon geometry to be important in dic- 
tating the QNM structure (and there is evidence in the Born series method 
which suggests this may be the case) Suneeta's results may be extended to 
generic dirty black holes as well. 



3.2.1 Schwarzschild black holes 

Matching the potentials 



V(r.) = (3.27) 



The potential of choice for modelling black holes has been the Poschl-Teller 
potential 

cosh^(Q;(r^ — r^)) 

where Vq is the value of the potential at the peak = r°. The "matching 
conditions" taken by Ferrari and Mashhoon were 

• That the maxima agreed: 



nWo) = V((r,)o). 



(3.28a) 
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a was determined by matching the second derivative at the peak: 



d^y 



dr^ 



d^V 



(r-*)o 



dr? 



(3.28b) 



(r-*)o 



These conditions emphasise the role of the peak of the potential. Instead of 
()3.28bj) . I adopt a convention that is easier to fit into multi- horizon space- 
times: 



That the rate of decay at the horizons are equal: 



r*->-oo r+ 



(3.28c) 



This condition determines a. 



Equation ()3.28c|) is inspired by Suneeta's approach to a SdS spacetime. At 
the end of this section my variation on the theme of the "Suneeta" approach 
will be compared with the results of Ferrari and Mashhoon, and with the 
numerical results. 




Figure 3.1: A comparison of the Regge- Wheeler potential (solid) and the 
Poschl- Teller potential (dotted) when the matching conditions ()3.28p are 
used. These graphs are, from the left, for £ = 0, 1 and 6 respectively. 
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Before starting a discussion on surface gravity, we find some concrete re- 
sults for the Schwarzschild black hole. The potential for a Schwarzschild 
black hole is (for a scalar field, cf ()2.112j) ) 

= (l - — j (— + j ■ (3-29) 

We find the location of the maximum: 

-4(^(^+l)r^ + 3M(l-£(£ + l))r-8M2) = (3.30) 
dr 

3[£(f + 1) - 1] + v/9[£(£ + 1) - 1]^ + 32f (£ + 1) 

= WTY) ^^-^'^ 

It is now simple to find Vq = V{ro), where tq is the location of the maximum; 
but such an expression is not illuminating. 



To impose the second condition we need to investigate the near horizon 
behaviour of the potential. Introducing the tortoise coordinate (cf ()2.65|) ): 

dr+ 1 



dr 1 - 2M/r ^^'^^^ 
= r -2M + 2M\n{r -2M) = X + 2M\n{x), x = r-2M. (3.33) 

Near the horizon, x is small. As x ^ ( or r^, — — oo) 

X ~ exp(r^/2M). (3.34) 

From equation ()3.29|) the near horizon form of the potential is found: 



X / X \ f ^ ?>x 



n^) = F77? 1-7T77 l-TTI7+^(^+l) 1 " 77 77^ , (3.35) 



8M3 V 2M/ V 2M 
or in terms of r^,: 



X 



M 



( X 



VM6 



Hence the near horizon decay rate is 1/2M, or more formally 

lim M]:>^ = ^. (3.37) 

r*^-oo 2M 
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The Poschl- Teller potential has exponential decay as one heads toward the 
horizon as well: 

lim V(n) = 7 ^ ^TT^ (3.38) 

r.~*~oo ^ *^ exp(-2a[r^ -r0])/2 ^ ^ 

= [2yoexp(-2ar°)] exp(2ar^). (3.39) 
The decay rates will agree near the horizon if and only if 

a = -i- = (3.40) 

AM ^ ^ 



where Kh is the surface gravity of a Schwarzschild black hole. 
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Finding the QNM 

Our attention is now turned to finding the QNM frequencies of the approxi- 
mate Poschl- Teller potential. That is, we wish to find the QNM given 



dr^ cosh^ far 



+ u'R = 27 T^' (3-41) 



where the axis has been translated so that r° = 0. The exact solution is 
given in terms of hypergeometric functions: 



R{r) = Ci{2cosh{2ar^) +2) 



l/i+G/a 



^ f-2iuj + a + AG a + 2iuj + AG AG I , 1 

^ \ Ta ' A-a ' ' + 2^' 2 ^^^^^'"^^^ + 2 

+ C2(2cosh(2ar,) + 2)1/^-^/" 

f2iuj + a-AG a-2iuj -AG AG 1 , , , 1 \ 

^ [ 4^ ' 4^ ' ^ - 2^' 2 + 2 j ' 

(3.42) 

G = ^^^Ja^-AV^. (3.43) 



^To find the surface gravity for any stationary black hole see (for example) Wald |79| . 
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We are wanting to look at the behaviour as goes to infinity. The argument 
of the hypergeometric is dominated by the exponentially growing piece: 

lim R{r) = 

r*— >itoo 

^ n , Ar^M fa-2iuj + 4G a + 2iuj + AG ^ , AG 1 ^.,\ 
C, expda + 4G||,,|/2),f, , ; 1 + -; -e-'-l j 

+ Cexpda - 4G||.r.|/2).f, (^i^^. ^^^^^1; 1 - ie-'-l) 

(3.44) 

To express this in terms of exponential functions (and hence ingoing and 
outgoing waves) near infinity the following result (see equation (e.6) of [52j ) 

2Fi(a, 6; c; z) = IjfJ^f — ^(-^)""2i^i ( a,a + 1 - c;a + 1 - b;- 
T[b)T{c — a) \ z 

To find the asymptotic form at infinity the series expansion for the hyperge- 
ometric function is employed: 

, , , ah z a(a + l)b(b + 1) z"^ , ^ 

.F.(a.t;c;.) = l + --+ J^\^ ' -+... (3.46) 

This allows us to replace the hypergeometric near infinity with unity. 
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The asymptotic form of ()3.42|) is then (up to an irrelevant constant fac- 



tor): 



hm R{r) oc CiF 1 + 

r— >±oo 



2a 



T{iuj/a) 



+ 



r{[a + 2iu + ^a^ -AVo]/Aa)r{[3a + Vi^^ - ^Vo + 2iuj]/4a) 

T{—iuj/a) 

T{[a + y/o^^^^Wo - 2itu]/4a)r([3a + V"^ - ^Vq - 2iuj]/Aaf 



-iuj\rt\ 



+ Car 1 



2a 



T{—iuj/a) 



-iuj\ri\ 



+ 



T{[a - 2iuj - V/t^ -4\/o]/4a)r([3a - 2iuj - ^ a^ - 4K)]/4 

V{iijj ja) 



r([a + 2icj - Va^^4K)]/4a)r([3a + 2iuj - ^a^-WQ\l^a) 



(3.47) 



The QNM boundary condition means that the incoming waves (the part 
proportional to exp(— zct;|r^|)) must vanish. This implies that 



= r - 



a 



CiF 1 



y/g^ - 4^0 
2a 



X 



4a 



4a 



^ a - 2it<; - ^/a^ ^ fia - 2iuo - ^a^ - 4Vo 



4a 



4a 



-1 



-iaj|r*| 



(3.48) 



It is well known that the gamma function has no zeros in the complex plane. 
Candidates for QNMs then come from the fact that the Gamma function 
diverges at the negative integers. Specifically, letting n be a positive integer, 
candidate QNM modes occur at 
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Ci = 0: Then we have one of the two remaining Gamma functions 
encountering a pole, i.e.: 



a 



2iu - ^Ja^ - 41^0 = -4na (3.49a) 



or 



3a - 2iuj - - = -Ana (3.49b) 

which imphes that 



o^n=\/^-(|) -tln + -]a. (3.49c) 



• 6*2 = 0: The only change is the sign of the square root changes. 
We can finally summarise the QNM for the Poschl- Teller potential: 



^„ = (±) ^IVo-i^y + z[n + ^ja. (3.50) 
3.2.2 Comparison of results 

In comparison, the results of Ferrari and Mashhoon [2^] are given by {i ^ 1) 

" ^±2(2£+l)+i fn + i^y (3.52) 



3v^ V V 2 

These results are presumably given for gravitational waves, although no such 
statement is made by the authors (electromagnetic QNM are independent of 
i, and the paper focuses on gravitational and electromagnetic QNM). 



It should be noted that the type of field and the angular momentum i 
determine the potential completely. In particular, they determine Vq, and it 
follows that the real part of Un will depend on i, irrespective of the matching 
condition chosen. 
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The values of n occur from the imposition of the boundary conditions on 
the Poschl-Teller potential. Hence the imaginary part will be evenly spaced 
if we match to the Poschl-Teller potential irrespective of the details of the 
matching conditions. The value of the spacing between the levels does depend 
on the matching condition. More specifically 

• Matching decay rate: Real part depends on I via Vq. The matching 
condition gives 

a — Kh (3.53) 
which is independent of i, hence the imaginary spacing is constant. 

• Ferrari and Mashhoon: Now a depends on £ as well: 

2M\ 1 



1 







dr 


r=ro 



(3.54) 



This condition is £ dependent, so the imaginary spacing would rely on 
both n and £. The real part has £ dependence from both Vq and a. 

It is worth emphasising again that the imaginary spacing of the real QNMs for 
a given £ does not have to be constant; this is a property of the approximate 
Poschl-Teller potential. 

The tabulation of the QNM mode frequencies show that, at least for low 
values of £, the decay rate of the potential allows much more accurate mod- 
elling than matching at the peak. It should be noted that the decay toward 
infinity is actually polynomial, as 

lim V(r) = ^ii±ll (3.55) 

1 — >oo T 

independent of the "spin" . At large r there is little difference between r and 



Mujn for a Schwarzschild black hole: 
Numerical and analytic approximations 



\ ^ 

n \ 


1 


2 


3 


4 


5 


6 







?+?z 


0.374 + 0.089i 


0.599+0.093Z 


0.809 + 0.0942 


? 


? 





1 


?+?z 


0.347 + 0.274i 


0.583+0.281Z 


0.797 + 0.2832 


? 


? 


0.375 


2 


?+?z 


0.301 + 0.478i 


0.552+0.4792 


0.772 + 0.480i 


? 


? 


0.625 


3 


?+?z 


0.252 + 0.705i 


0.512+0.690Z 


0.740 + 0.6842 


? 


? 


0.875 


4 


?+?z 


0.208 + 0.947i 


0.470+0.9162 


0.702 + 0.8982 


? 


? 


1.125 


5 


?+?z 


0.169 + 1.196i 


? + ?2 


? 


? 


? 


1.375 


200 


?+?z 


0.061 + 49.860i 


0.078 + 49.843i 


0.085 + 49.8372 


? 


0.241 + 49.397Z 


50.125 


400 


?+?z 


0.056 + 99.864^ 


0.069 + 99.852Z 


0.083 + 99.8372 


? 


0.129 + 99.678Z 


100.125 


600 


?+?z 


0.054 + 149.8662 


0.065 + 149.8562 


0.078 + 149.8432 


? 


0.0360 + 149.548i 


150.125 


800 


?+?z 


0.052 + 199.867i 


0.062 + 199.8582 


0.074 + 199.8472 


? 


0.080 + 199.817 


200.125 


1000 


?+?z 


0.051 + 249.866Z 


0.060 + 249.860i 


0.072 + 249.8502 


? 


0.086 + 249.824i 


250.125 


3000 


?+?z 


0.048 + 749.8712 


0.053 + 749.8662 


0.060 + 749.86O2 


? 


0.077 + 749.844i 


750.125 


6000 


?+?z 


0.0469 + 1499.872i 


0.051 + 1499.8682 


0.056 + 1499.864i 


? 


0.068 + 1499.852i 


1500.125 


10000 


?+?z 


0.0462 + 2499.873i 


0.049+ 2499.870i 


0.053 + 2499.867i 


? 


0.063 + 2499.8571 


2500.125 




0.060 
0.108 


0.368 
0.378 


0.597 
0.602 


0.807 
0.811 


1.011 
1.014 


1.211 
1.214 





Table 3.1: The modes n = 0,1,2,3,4,5 were obtained from Anderson and Glampedakis [SI]. The higher 
modes were taken from work kindly provided by Emanuele Berti (private communication). Here Us are the 
frequencies by imposing fl3.28c|) (the "Suneeta" condition), while Uf comes from imposing ()3.28bjl (the "Ferrari" 
condition). Question marks denote frequencies that I have not found numerical results for. The real parts of 
both Poschl- Teller approximations are independent of n. The imaginary part of Ug is independent of i as it 
is related to the surface gravity; but as Uf relies on the curvature of the (£ dependent) potential the '^{ujf) 
changes as i does. 



Mun for a Schwarzchild black hole: 
Using the Mashhoon & Ferrari matching condition 





1 


2 


3 


4 


5 


6 


100 


500 


n = 


0.054z 


0.189i 


0.301z 


0.406z 


0.507z 


0.607i 


9.669i 


48.160? 


1 


0.163i 


0.567i 


0.904z 


1.217z 


1.521i 


1.820i 


29.008Z 


144.481? 


2 


0.2717: 


0.946i 


1.506i 


2.029i 


2.535i 


3.034i 


48.347?; 


240.802? 


3 


0.379i 


1.324j 


2.108? 


2.840i 


3.549i 


4.247i 


67.686? 


337.123 


4 


0.488i 


1.702i 


2.711i 


3.651i 


4.564i 


5.461i 


87.025? 


■1 O O AAA' 

433.444? 


5 


0.59m 


o no 1 ■ 

2.081i 


6.6l6i 


4.4d3« 


5.578i 




106.364? 


529.764? 


200 


21.738i 


75.844i 


120.785i 


162.692i 


203.330i 


243.315i 


3877i 


19312ii 


400 


43.422i 


151.498ii 


241.268i 


324.978i 


406.154ii 


486.023ii 


7745i 


38576ii 


600 


65.105i 


227.153i 


361.751i 


487.264i 


608.977i 


728.731i 


11612? 


57841? 


800 


86.789i 


302.808Z 


482.235Z 


649.55n 


811.80n 


971.440Z 


15481? 


77105? 


1000 


108.473Z 


378.462Z 


602.718Z 


811.837i 


1014.624Z 


1214.148Z 


19349? 


96369? 


3000 


325.31H 


1135.008Z 


1807.552Z 


2434.699i 


3042.858Z 


3641.230Z 


58026? 


289011? 


6000 


650.568Z 


2269.827i 


3614.803Z 


4868.992i 


6085.209Z 


7281.853Z 


116042? 


577973? 


10000 


1084.244Z 


3782.919Z 


6024.471Z 


8114.716? 


10141.677i 


12136.017Z 


193398? 


963256? 




0.108 


0.378 


0.602 


0.811 


1.014 


1.214 


19.339 


96.321 



Table 3.2: This table shows the QNM calculated using Ferrari and Mashhoon's method. The real part depends 
only on and is displayed at the bottom of the table to prevent needless duplication. At constant ^ the spacing 
between adjacent n levels is constant. 



Chapter 4 



Near horizon conformal field 
theory 

When entropy was first introduced into classical physics, it was as a state 
function and bore no relation to "microscopic degrees of freedom" until the 
pioneering work of Boltzmann. It is then perhaps not completely surprising 
that the explanation and calculation for black hole entropy yields an answer 
but does not enlighten us as to the microscopic origins of black hole en- 
tropy. What is surprising is that the answer is definite; in classical physics 
the entropy was determined only up to an additive constant, whereas the 
prescription of the logarithm of the number of microstates yields an exact 
value. However the calculations of black hole entropy do not yield any clue 
to the nature, location or even existence of black hole microstates.^ 

The black hole temperature and entropy calculations rely on very basic 
physics, and can easily be extended into other dimensions. This has led 
to the hypothesis of universality:'^ any explanation of black hole entropy 

^It may be argued that the existence of entropy is proof of the number of microstates. 
It is possible that black holes are fundamentally different from other systems and possess 
entropy without microstates but this is a fairly radical and desperate idea. 

^Universality is overused in the context of black hole entropy. Carlip |15j uses univer- 
sality to mean that many different theories of quantum gravity reproduce these results 
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must also apply in any number of spacetime dimensions. This may seem 
puzzling at first - we only live in four spacetime dimensions so the fact that 
the calculation may not work in other dimensions seems largely irrelevant. 
Having the argument work in only four dimensions could actually be seen 
as a positive sign if there was some fundamental reason for it, as it may 
be an "explanation" for why the world is four dimensional. The problem is 
that general relativity is written in the language of differential geometry, a 
language that is very flexible about the dimensionality and type of manifold 
that it is placed on. Consequently we can take our semi-classical entropy and 
temperature calculations and apply them to mathematically legitimate (but 
presumably unphysical) black holes. Given that these same semi-classical 
calculations provide the only justification for us to look for microstates of 
black holes it would be odd if the unphysical spacetimes did not posses similar 
microstates; hence the universality principle. It is only a hypothesis, and it 
may just be (however unlikely) that quantum gravity links quantum field 
theory and general relativity in such a way that the only consistent theory is 
in four dimensions; hence in other dimensions the classical (and hence semi- 
classical) approximations break down. Until evidence of such a claim can be 
provided (or indeed suggested) by any calculation, the universality principle 
seems like a good guiding principle. 

We also have a large number of quantum gravity candidates that are able 
to reproduce the correct temperature and area of a black hole such as string 
theory, D-branes, loop quantum gravity (LQG) and Sakharov induced grav- 
ity models. In black holes with negative cosmological constant (anti-de Sit- 
ter space or AdS) there is a correspondence with conformal field theories 
(CFTs) that gives the correct entropy and temperature as well. An excellent 
and recent review of calculations of black hole entropy that shall be used 
extensively throughout this chapter is Fursaev [21]. None of these models 

hence they are universal. Fursaev |29j calls a method universal if it gives the correct result 
for any spacetime. 
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is entirely satisfactory. The LQG calculations were outlined in chapter |2l 
where the microstates were spin-links that penetrated the horizon. String 
theory calculations are often done in a spacetime that is dual to the original 
spacetime. As such the microstates are not even located in the right space! It 
is not clear that the string theory calculations can be carried out on a generic 
black hole; see Fursaev [SHj for details. The fact that the entropy and tem- 
perature seems generic across many different theories seems to suggest one 
of the three possible alternatives: 

• It is just a numerical coincidence. Any quantum gravity calculation 
that does not give the Hawking temperature or Bekenstein-Hawking 
entropy maybe discarded. It is just an odd curiosity that so many 
theories remain. 

• These theories are not really different. They are the same underlying 
theory written in a different way, and they are formally equivalent once 
correct identifications are made. 

• Something fundamental is being assumed in all of the models. One 
such suggestion is that a symmetry in the classical theory is being 
preserved under all these quantisations making the temperature and 
entropy calculations generic. 

To some degree the first idea must be correct, as it is possible to concoct the- 
ories which do not give rise to the correct entropy or temperature relations 
(e.g. the completely indistinguishable spin-links in Polychronakos [HI])- To 
claim that there is no reason for the existence of so many completely dif- 
ferent theories that reproduce these results seems highly unlikely. A literal 
interpretation of the second idea is demonstrably false; among string theo- 
ries alone many are known to be inequivalent! In the spirit of M-theory it is 
possible that many of the theories are different ways of writing one theory in 
different languages and in different limits. 
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This chapter is written as part of an investigation into the third possibihty, 
that a near horizon classical conformal symmetry is preserved on quantisa- 
tion and fixes the density of states and hence the entropy. This idea started 
when Brown and Henneaux ^U] noted in classical AdSs that the generators 
of the group of diffeomorphisms at large r is becomes (two copies of) the gen- 
erators of CFTs, known as Virasoro algebras. Brown and Henneaux where 
able to conclude that any theory of quantum gravity on AdSa that was diffeo- 
morphism invariant must be a conformal field theory! In two dimensions the 
group of conformal transformations is infinite dimensional (any homomor- 
phic function is a conformal transformation) and the conformal symmetry 
is a powerful constraint; it determines the density of states up to a con- 
stant known as the central charge. In classical theories without boundaries 
the central charge vanishes, and hence it is normally associated with quan- 
tum theories. However Brown and Henneaux showed that even classical 
field theories could acquire a central change if they had a boundary. In the 
AdSs calculation this boundary was taken to be the t — (f plane at infinite 
radial distance. Strominger [21] showed that it was possible to obtain the 
Hawking-Bekenstein entropy of a negative cosmological constant black hole 
from the conformal field theories. In higher dimensional spaces the CFTs 
are finite dimensional and do not fix the density of states, and it was unclear 
how to perform direct generalisations. Even in three dimensions the CFT 
was unable to attribute individual black holes with an entropy as the CFT 
correspondence only applied at large r. I think that this problem is some- 
what overstated as entropy cannot be localised; entangled particles can fly 
arbitrarily far apart in a "pure state" without increasing the entropy. When 
I allow one particle to interact with a heat bath I change the entropy of the 
entangled pair. Naturally one would like to be able to associate entropy with 
(essentially) non-interacting subsystems of the universe, but if black holes 
are essentially non-interacting then the asymptotic symmetry should be "es- 
sentially" obtained far enough away from each black hole that each one could 
be assigned an entropy. It is not a forgone conclusion that it is possible to 
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associate entropy with individual black holes if they are strongly interacting 
(such as in a binary orbit). A more serious objection is that the method of 
using a CFT at infinity cannot distinguish between a spacetime that contains 
a star and a spacetime that contains a black hole. 

To have a CFT account for the black hole microstates in a generic space- 
time there must be a boundary, and a 2D surface on which the generators of 
diffeomorphisms look like (multiple copies of) Virasoro algebras. In generic 
black hole spacetimes one such candidate for the boundary is a horizon [T^. 
As the horizon has a dimensionality one lower than the spacetime any bound- 
ary conditions specified on the horizon eliminate gauge freedom there. The 
question of what type of horizon needs to be used is not resolved, but for 
the stationary case where a Killing field is available a Killing horizon is ap- 
propriate. In the stationary case the two plane spanned by the normal to 
the (spacelike section of the) horizon and the Killing field possess a 2D CFT 
as the horizon is approached |14^ 1701 ITKj ITT)] . As the CFT is on the horizon, 
the microstates can be associated with an individual black hole. It should 
be emphasised that the black hole horizon is not being treated as a phys- 
ical boundary, but as a convenient place to set boundary conditions that 
ensure the existence of a black hole (and hence in the stationary 2D 
CFT) ^3]. As a consequence the "near horizon" CFTs can distinguish be- 
tween stationary black holes (which posses 2D CFTs) and stationary stars 
(which generically do not). 
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Chapter 5 

Near horizon stress tensors 



5.1 Stationary metrics 

In this chapter we investigate stationary metrics, meaning that there is an 
isometry of the spacetime, the orbits of which are timehke curves. We have 
two possibihties for a spacetime containing a black hole: 

• The spacetime is static; that is the timelike (Killing) vector field is 
orthogonal to a family of spacelike surfaces of "constant time" . This can 
be shown to be equivalent to the statement that there exist coordinates 
in which all mixed space/time metric coefficients of the form gu vanish 
(i = 1,2,3). See for example Wald, page 119 

• The spacetime is stationary and axisymmetric. An axisymmetric space- 
time is a spacetime with an isometry that has closed spacelike curves. 
For a spacetime to be both stationary and axisymmetric the action of 
the isometries must also commute. 

There is actually a third case that is less often discussed: 

• The spacetime is static and non- axisymmetric. For this to occur, an 
ergosphere would have to exist and it could not intersect the horizon 
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There are no known solutions of the third type, and good reasons to suppose 
that they do not exist. Hawking and Elhs (41J point out that a "Penrose 
process" allows the existence of negative energy particles (as measured by 
an observer at infinity) as long as there is an ergosphere.^ Hawking and 
Ellis then argue that if the horizon and the ergosphere have no point in 
common then the negative energy particles cannot pass through the event 
horizon but must remain in the ergosphere. If it was possible to continue 
this process indefinitely we could mine an infinite amount of energy from the 
black hole. The ergosphere is constrained from spontaneously vanishing, as 
we require a place for the negative energy particles. The other possibility is 
that the event horizon and the ergosphere eventually touch, allowing particles 
to escape from the ergosphere into the singularity. Once there are no negative 
energy particles the ergosphere can be reduced to zero, as any negative energy 
particles can be forced into the black hole. 

The above arguments come from §9.3 of Hawking and Ellis, but do not 
completely preclude the idea that the ergosphere is disconnected from the 
event horizon. For example, the Penrose process could be limited by the 
gravitational effects of the numerous test particles fired in to extract energy! 
There is a stronger statement by Hajicek [HH! that insists that the ergosphere 
and event horizon must coincide. Hajicek's result requires the mild technical 
assumption that the orbits of the timelike Killing field are complete. 

Once we have discarded the possibility that the ergosphere and the event 
horizon do not intersect, we have to then show that a non-axisymmetric 
metric does not exist. 

"'^To see this, remember that if the four momentum of a particle is k"' and the asymptot- 
ically timelike Killing vector is then the energy of the particle measured by an observer 
at infinity is E — —k°-C,a- When becomes spacelike, E may be less than zero. But C" 
spacelike defines the position of the ergosphere, hence negative energy particles can only 
exist in the ergosphere. 
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5.2 The ADM decomposition 

The Arnowitt-Deser-Misner (ADM) decomposition jl] is the process of break- 
ing a four-dimensional spacetime into a time dependent three-dimensional 
geometry. To do this the spacetime (or a region of the spacetime) is foliated 
with three dimensional spacelike surfaces, labelled by a continuous param- 
eter t. The idea is that the surface represents the geometry of space at 
"time" t. There is no natural way of doing this even in special relativity, 
as even if one observer sees two spatially separated events as simultaneous 
other observers do not necessarily agree! The decomposition of spacetime 
into space and time (or more precisely, a series of spacelike surfaces labelled 
by a parameter t) is thus highly non-unique. For obvious reasons the ADM 
decomposition is also referred to as a (3 + 1) decomposition. While there 
are many references available on the ADM decomposition, such as the living 
review article or I have found an elementary treatment lacking. Thus 
the work in this section is meant to be clarifying and does not contain any 
new results.^ 

What do points on different surfaces Sj^ and have to do with one 
another? By introducing a timelike vector field representing the "flow of 
time" we are able to identify points on with by insisting that if two 
points are on the same orbit as r'* then they are the same point (see figure 
15. 2j) . Instead of thinking of spacetime as a foliation of spacelike hypersurfaces 
we can instead consider space to be the surface Sq with a time dependent 
metric. Notice that the vector field is highly non-unique, and different 
choices lead to different decompositions. 

To actually carry out the ADM decomposition we need to make three 
(non-canonical) choices: 

^One of the examiners (David Wiltshire) has pointed out some elementary resources 
that I had not found. He suggests the references |37[ I72L IST| . 
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Figure 5.1: Two figures showing the dragging of the points of Eq along 
integral curves of the fiow r". The time function t is different on the LHS 
and RHS. The level surface t — T (i.e. Et) on the LHS is spacelike and 
well-behaved, if this occurs for all T then i is a vahd time function. On the 
RHS the surface E^- is not spacelike so t is not a valid time function. 

• An initial surface Eq, representing time zero. 

• A time function t : M ^ that assigns a time to each point in the 
spacetime. To be consistent with the interpretation of Sq we require 
that Vp e So : t{p) = 0. 

• A flow of time r". 

A change in any of these will lead to a different decomposition of the same 
spacetime. Notice that the choices are not all independent or arbitrary. For 
example, the time function t is constrained so that the surfaces of constant 
t (i.e. Ef) are spacelike. A function t that satisfies this constraint is called a 
valid time function. Normally the normalisation condition 



is insisted upon. I shall not insist on it, as becomes null as the horizon is 
approached. 



T«V„i = 1 



(5.1) 



Once we have defined the time function t we have a future pointing timelike 
vector field, Vt. Note that there is no a priori reason to pick so that it 
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points in the same direction as Vt. In general, we can decompose the flow 
into a timelike normal part and a spacelike part. First, we introduce a 
timelike normal 

and then decompose the flow of time in the following manner: 

= Nn"" + AT" (5.3) 

where is a scalar function called the lapse and A^'* is a spacelike vector 
called the shift. In an intuitive sense N°- lies in S^. Of course, this cannot be 
taken literally as and A^'^ both lie in the tangent space of a point in M, 
not in M itself which is where is deflned! However, this intuitive idea is 
captured by the orthogonality requirement of the normal and the shift: 

N'^Ha = 0. (5.4) 

Figure 15.21 shows this decomposition. Getting the lapse is easy, as one only 
needs to flnd the component of along n": 

A^ = -r"na, (5.5) 
or by using the normalisation condition ()5.ip 

N = +^= ^ (5.6) 

The shift can be obtained by a rearrangement of ()5.3p . However, it is nicer to 
introduce a projector that will allow us to flnd the component of an arbitrary 
vector orthogonal to n". First we note that the component of parallel 
to is 

"^along n = ~'>^aS"' (5.7a) 

and hence the orthogonal component is 

^Ito n = 5" - <iong nn" = + n,s'n\ (5.7b) 
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Figure 5.2: The decomposition of into the shift and lapse. Picture inspired 
by Wald 



Finally we note the identity 

= g\s'' (5.7c) 

and hence 

si,, ^ = i9\ + n^n,) s\ (5.7d) 

The quantity in brackets projects an arbitrary vector to the spacelike surface 
orthogonal to n", so we define 

h\ = g\ + n''nk. (5.7e) 

In particular, we can apply the projector to the flow vector fleld to obtain 
the shift: 

= /i^Vfo. (5.8) 



Now we give some insight into the meaning of the shift, lapse and projector. 
Let us consider a timelike observer moving along an integral curve of n". The 
distance measured by such an observer is given by 

ds^ = -N^dt\ (5.9) 

Using A as the proper time of such an observer we have 

dt = ^. (5.10) 



145 



Thus a small lapse means that a observers travelling along integral curves 
of 71°' ("normal observers") experience a small "lapse" in time between two 
constant time surfaces. For example, if t is the Killing time coordinate it 
takes an infinite amount of Killing time to cross an event horizon, while the 
proper time of a "stationary" observer slows as the horizon is approached^. 
When using Killing time as the time function we then when an event horizon 
exists we have iV = 0. 

Note: There is an error in the paper reproduced in AppendixKj\ and 
the published version. We strongly imply that the statements "the 
horizon", "N = 0" and "infinite red-shift surface" can be used inter- 
changeably. This is not the case: the ergosphere is an infinite redshift 
surface (with respect to Killing time) but has no special causal prop- 
erties. 

What is the interpretation for the shift? Let {t = ) be the 

coordinates in a patch of the four dimensional manifold M. Let {y^,y'^,y^) 
be the three coordinates on a patch of Sj. Now consider and St+dt- By 
construction, we are identifying the point y E T,t with the point in T,t+dt if 
they lie on the same integral curve of r". In fact, we can see the stronger 
statement 



With a coordinate change, we can find the metric in terms of the coordinates 
{t = y^, y^, y"^, y^). We have 



ds^ = gtt + 2gudt {Ay' + iVMt) + g,,{dy' + N'dt){dy^ + iVMt) (5.13) 



■^I really want an apparent horizon here, rather than an event horizon. But as I have 
used a Killing time I have assumed the spacetime is stationary and hence the apparent 
horizon and event horizons coincide. 



dy'' = dx° - N°dt. 




ds^ = gabdx''dx'' = gabdy°dy' 



(5.12) 



or, in terms of the y coordinate system 
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where i and j run over spacelike indices only. But the second term is zero, 
as the t direction is orthogonal to the spacelike directions by construction. 
The last term has sums over dy"' and A^**, both of which are orthogonal to 
n"'. Thus there is no different between using gij or hij and we finally obtain 

ds^ = Qu dt^ + hijidy' + N'dt){dy^ + N^t) (5.14) 
= {gtt + KjN'N^) df + 2{KjN')dy' dt + h,,dy' dyK (5.15) 

Thus we see that the projector hab is really the time dependent metric on 
the three dimensional geometry S^. To identify the lapse we look again at 
the normal observers 

ds^ = -N'^dt^ = {gtt + hijN'N^) dt^ (5.16) 

or that 

g,, = -{N^ + h,,N'N^). (5.17) 
Hence we can write the metric in the terms of the shift and lapse by 





/ 


-(N^ + hijN'N^) 


Ni 


9ab — 




iV, 


hij 






V 







5.2.1 Modified ADM for stationary spacetimes with 
horizons 

A glance at eq. (G.IO) would lead one to think naively that the lapse is given 
by and the shift hj uj dip, based on a simple comparison with the form of 
the metric in ()5.18|) . It is true, but quite difficult to prove, that = is the 
horizon. The difficulty stems from showing that the geodesies are trapped, 
rather than just having a perverse coordinate system that "shifts" close to 
the shifts at the speed of light. 
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To see this, we note that we cannot use the Kilhng vector to generate 
time translations as ( is spacehke in the vicinity of a rotating black holes 
horizon. In fact, this is the definition of the ergosphere! As the vector field 

is not timelike this also implies that I cannot use the Killing time t as the 
parameter labelling different surfaces. To distinguish the Killing coordinate t 
and the time function I shall denote the latter by T. There is a Killing vector 
field that is null on the horizon and timelike in the vicinity of the horizon, 
formed by a linear combination of the time translation Killing vector field 
and the axisymmetric Killing field i/j"". 

x'' = C + n,r- (5.19) 

As a linear combination of Killing vector fields, is also a Killing vector 
field. It is this Killing field that generates the horizon. The constant is 
the rotational rate of the black hole, defined by taking the locally measured 
rotation 

uj = (5.20) 

in the limit as the horizon is approached. The result that the rate of rotation 
is constant along the horizon is the rigidity theorem due to Carter [T7] . 



We could use x"^ as our flow of time vector field, but that is not timelike 
at large distances. Instead, we use 

r" = C + UJ t/j". (5.21) 

Notice that is not a Killing vector field. As uj is the locally measured rate 
of rotation, it does have the physically nice property of identifying a point 
on one time slice with a point on another time slice that has "drifted" along 
with the space. We also choose the time function so that r = VT. Hence the 
flow and the normal are pointing in the same direction and the shift vanishes. 
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The lapse is given by 



= gtt + 2u gt^ + u^g^^ 



a 



(5.22) 
(5.23) 




(5.24) 



(5.25) 



5.3 Near horizon geometry for stationary Killing 
horizons 

5.3.1 Metric for stationary non-static Killing horizons 

In this section we finally tie together the loose threads developed in previous 
sections. First we use the assumption of an stationary (but not static) metric 
with a horizon; then from ^5.11 we find the spacetime is also axisymmetric. 
Hence we may introduce coordinates (t, xi, X2) so that the metric functions 
g^i, are independent of t and y?. Then the condition that the spacetime 
is time reversible is imposed; that is the spacetime is invariant under the 
simultaneous changes 



As the metric components themselves are independent of t and but the 
intervals are not, distances will be preserved iff 



t 



-t, 



(5.26) 



0, 



i = 1,2. 



(5.27) 



Conditions that ensure the existence of a time reversal symmetry are dis- 
cussed in ^5.3.11 Equivalently, time reversal symmetry can be assumed as a 
"physically reasonable condition" . Together these conditions ensure that the 
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metric can be taken locally as 



















(5.28) 



v 















gxiXi{Xi,X2) gxiX2{Xl,X2) 
gx2Xi{Xi,X2) gx2X2{Xl,X2) J 



As the spacetime decomposes into the direct sum of two 2-d manifolds, an 
index convention that respects this is now adopted. Greek indices a,(3, . . . 
will represent a general spacetime index; that is it may take values from to 
3 representing the coordinates (t, ip, xi, X2). Upper case latin indices A,B,... 
take the values or 1 and represent the t — (f submanifold. Lower case latin 
indices a,b, . . . take the value 2 or 3 and represent the coordinates Xi or X2 
respectively. This convention is not used in earlier parts of this thesis. Such 
a convention allows us to write ()5.28|) in the much more compact form 



The ADM-like decomposition described at the end of ^5.21 is performed on 
the t — (fi submanifold. That is, 



where the horizon is located at = 0. The parameter u defines how much 
the horizon "rotates" by as one travels along the orbits of x", the vector field 
that is null on the horizon. We have the earlier result that u on the horizon 
is Qh- 

So far no use has been made of the coordinate freedom in xi and X2- Near 
the horizon we can introduce a normal coordinate system. The coordinate n 




(5.29) 




(5.30) 



(5.31) 
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measures the distance normal from the horizon, with the horizon occurring 
at n = 0. The coordinate z is a spacehke vector that hes on the horizon. 
The decomposition is then 



9ab 




(5.32) 



Notice that the (now overused) word normal in this context means a spacelike 
vector that lies in the X1-X2 plane that is orthogonal to the horizon. For 
example, in the Schwarzschild geometry the radial coordinate r would be a 
normal coordinate in this sense of the word. The full metric can now be 
given 
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(5.33) 



Note that this coordinate system only exists locally, it will fail as soon as 
two of the geodesies leaving the horizon in the n direction intersect. As we 
are only interested in the near horizon results, the local construction is good 
enough. 



5.3.2 Static metrics 

If the metric is static then we no longer have the result that it is also ax- 
isymmetric. Physically this is because there is no rotation, and hence no 
tidal friction. As we know there that the KiUing vector field is hypersurface 
orthogonal we have the result 
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(5.34) 
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In the neighbourhood of any point on the horizon we may decompose the 
coordinates into an on horizon part (99, z) and a normal coordinate n. The 
metric then takes the form 
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(5.35) 



This exhausts the four coordinate choices we have available (choosing Killing 
time t, picking the normal n, and choosing z and ip so that = gmp = 
0). Without the introduction of more symmetries no further reduction is 
possible. 



5.4 Out, damned singularity! 

Once we have adopted appropriate coordinates (either for a static or station- 
ary non-static spacetime) near the horizon, we can insist that the curvature 
scalars remain finite as the horizon is approached. Specifically, we require 
that 

_ pab r? 

• -ft -Kob 

• and RabcdR""^ 

all remain finite as the horizon is approached. (Any linear combination with 
constant coefficients will also suffice, as a computational matter it is some- 
times more convenient to use the traceless versions of the Ricci tensor and 
Riemann tensor). If any of these scalars becomes infinite as the horizon 
is approached then a genuine curvature singularity exists, as a scalar must 
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be independent of the coordinates used. It is possible to have a directional 
curvature singularity, where any tetrad will have directions distorted as the 
singularity is approached. However, it can happen that the curvature in 
two or more independent directions can contribute opposite amounts to the 
curvature scalars, rendering the limit finite. Our finiteness condition is thus 
necessary but it is not sufficient. 

In order to investigate the near-horizon geometry we take n = as the 
horizon and expand the lapse N and the rotation rate cu: 



Naturally, uj is trivial in the case of static black holes! We have used the 
additional information that n = 0, the horizon, corresponds to iV = and 
that the rigidity theorem implies that a;(n = 0) = 

The rest of the details are simply a case of calculating the Ricci scalar, 
Ricci tensor squared and the Weyl tensor squared. 1 summarise the case for 
the stationary non-static case below. The static case is similar, and in both 
cases more details can be found in the appendices 

5.4.1 The curvature invariants for stationary non-static 



The Ricci Scalar 

The Ricci scalar for a stationary non-static black hole at the horizon is given 



N{n, z) = Kh{z) Ki{z)n —K2{z)v? ^ ... 
uj{n, z) ^flh + uji{z) + ^uj2{z)n^ + ... 



(5.36) 



(5.37) 



horizons 



by 




(5.38) 



Assume: 
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• Non-extremal black hole: K.h{z) ^ 0. 

• The horizon geometry is regular: g(<^(^ 7^ and is finite 
Hence we must have lo\{z) = so that R remains finite as n ^ 0. 



The traceless Ricci tensor 



After using the fact that 0Jx{z) = we find: 



ah 



RabR' 



:R' 



dln^r^^ 



gzz{^,z)Kh{z) \ dz 



n=0 



n=0 



d In g^^. 



„=o 

4 + 0^1. 



n 



(5.39) 



As this is a sum of squares, we may conclude 



• An alternative derivation of the zeroth law: the surface Kh is constant. 

• The next term in the expansion of the lapse, k,2{z), vanishes. 

• Both dnQzz and dnQip^ vanish on the horizon. 



The other constraints 

We have yet to apply the constraint that the Weyl tensor squared yields 
a finite result, or use the terms in the Ricci scalar and square of the 
traceless Ricci tensor. The six constraints found above render the square of 
the Weyl tensor finite at the horizon upon substitution, and also ensure the 
rT^ pieces of the Ricci scalar and traceless Ricci tensor vanish as well. 
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5.5 Near horizon geometry 

In all the stationary black hole cases, when we construct a tetrad from 
(x, n, (/?, C) we have the result that the Einstein block diagonahses 
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(5.40) 
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The exact form of the different components depends if you are looking at 
the static case or the stationary non-static case. However, the general form 
of the stress tensor remains in tact. It is at this point, and this point only, 
that we can invoke the Einstein equations to tell us that a Killing horizon 
requires that 

P = -Pn (5.41) 

if p is the energy density in the timelike direction In particular, this 
conclusion is reached without any consideration for the energy conditions. 



Chapter 6 
Conclusions 



With a range of problems so diverse, it is hard to think of a concise way of 
summarising all of them. Indeed, to a large extent I already have summarised 
them at the end of each section. Instead of explaining the problems again 
and how far through completion they are, I would like to leave with a list 
of things that I believe are easily do-able but have not yet been done. I am 
afraid that to the reader that this may read as a checklist of things to do, 
and that is indeed how I intend to use it. But after providing a 250+ page 
thesis, I hope that the reader will allow me this indulgence. 

Algorithmic constructions 

• I believe that looking at spherically symmetric but anisotropic stars is 
going to be a hard problem, and likely one of little physical relevance. 
While I believe that very compact objects such as quark stars or neu- 
tron stars are probably going to be very anisotropic, especially at the 
surface where crusts are expected to form, the assumption that they 
are static will almost certainly be false. Even a very slowly spinning 
star will have a large angular momentum, resulting in a very quickly 
rotating neutron star that has formed as the result of gravitational 
collapse. 
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• Finding a bound on the compactness for isotropic stars that also obeys 
the DEC looks like a simple problem. The goal would be to find a family 
of stars that get arbitrarily close to saturating the 8/9ths bound, or to 
improve the bound and saturate that. 

• Improving the variational changes method in the gravity profile to allow 
for arbitrarily large changes in the profile, so that more comparison 
bounds can be done. As can be seen from appendix El many bounds 
were obtained by comparison with the Schwarzschild solution. 

• The problems of ensuring the positivity of integrals, such as the ex- 
pression for p{r), seem too difficult. A more modest problem of finding 
large classes of functions g{r) that ensure the energy density is positive 
everywhere may, however, be achievable. 

Quasinormal modes 

This section is somewhat depressing for two reasons. The first is that the 
motivation for looking at highly damped QNMs does not exist anymore, as 
the recent paper by Domagala and Lewandowski showed that the QNM con- 
jecture was almost certainly false j21j. The other problem is that everything 
is very messy; one may consider that the adjective "dirty" could refer to 
finding the potentials as well as or better than to the spacetime! Finding an 
approach that is more intuitive or concrete would make me feel better than 
just going through screeds of algebra. This may, of course, be impossible. 
After all, the nice geometric results tend to span classes of spacetimes and 
tend to be existence theorems [for example, singularity theorems]. Here we 
really are interested in numerical answers for QNM, not existence answers to 
existence questions such there an infinite number of QNMs for every 

black hole spacetime"? And as a numerical question, there is some expec- 
tation that the answer will come about from arduous working for different 
spacetimes. There is one avenue that seems worthwhile pursuing, namely 
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• Can we express the problem in a clean and separable way in the Newman- 
Penrose formalism? If so, can we assign a similar interpretation as in 
the Schwarzschild case? 

The NP formalism allows separation of Reissner-Nordstrom, Kerr, Kerr- 
Newman and of course Schwarzschild. It is easier to see the relationship 
between the spectrum of the axial and polar cases when in the NP formalism, 
and indeed it provides a clear way to show that the two cases are isospectral 
in these examples. As we know that general dirty black hole spacetimes are 
not isospectral, we can ask the related question 

• Is the utility of the Newman-Penrose formalism only significant if the 
spectra under consideration are isospectral? 

A slightly different, but related question, would be 

• How useful is the NP formalism in a situation where the standard 
massless fields are isospectral (as, excluding the scalar field, massless 
fields only possess two polarisations) at describing a massive field? Do 
all polarisations occur independently? Are all polarisations isospectral? 

It would also be interesting and feasible, although algebraically tedious, to 
extend the potentials for the standard fields to the stationary dirty black 
holes used in the conformal symmetry section rather than just static and 
dirty black holes. 

Although it is not something that I plan to spend much time on, I would 

be interested in at least seeing if the claims about the real part of highly 

damped QNM changes now that the QNM conjecture is almost certainly 

false. Do these analytical proofs hold up to closer scrutiny, or are physicists 

just as prone to the bandwagon effect as sociologists^? 

^Physicist Alan Sokal managed to get a paper entitled "Transgressing the bound- 
aries: Toward a Transformative Hermeneutics of Quantum Gravity" published in 
Social Text #46/47 pp. 217-252 (1996). It is still available onUne from 
http://www.physics.nyu.edu/faculty/sokal. 
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Conformal symmetry 

There are two things that I would hke to work through in more detail when it 
comes to the work that we have done on conformal transformations. The first 
is to see if I can find a different coordinate system in which both the central 
charge c and the eigenvalue of Lq remain finite as the horizon is approached. 
This would suggest that it is possible to have a well behaved CFT on the 
horizon, rather than having a theory which one is only allowed to answer 
certain questions for which well defined limits exist. To do this I would want 
to have a better understanding of the role that a central charge plays. Steve 
Car lip has suggested looking at a (1 + 1) Dilaton black hole and looking at 
null sheets to see if a well defined central charge can be obtained in this toy 
model. 

The other thing that would be interesting does not relate to the CFT, but 
rather simply the form of the geometry near the horizon. We have a null 
surface, and by hypothesis of stationarity the fields on that surface cannot 
evolve in time. One can imagine that even if this null surface was not a 
Killing horizon that the stress tensor would be constrained - after all massive 
particles cannot travel along a null surface! So while intuitively we may 
expect that on a null surface only "zero mass" fields would be able to exist, 
as massive "quanta" would be forced to travel at the speed of fight to remain 
on this surface, I face two difficulties. The first is that it is difficult to express 
in general terms what is meant by a field with zero mass quanta in terms of 
allowed stress energies. At the moment I am not playing with a Lagrangian, 
only the locally defined stress energy tensor. The second point is slightly 
more problematic: the energy conditions have not been imposed anywhere! 
A timelike particle is still required to have a non-null four momentum, and so 
it is still possible that the only way to make such quanta have null momenta 
is to have singularities. This is an interesting problem, as it indicates that 
the restrictions on the stress energy we are seeing may have no bearing on 
the question of the existence of a CFT near a Killing horizon! Still, as an 
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idea it is still too vague for me to use to prove anything with. 

Finally, as the near-horizon question is more one of existence 1 would hope 
that a nice clean argument could be presented that gives the same result. A 
candidate for such a result is the bifurcate Killing horizon argument outlined 
in appendix O 
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Appendix A 



Algorithmic construction of 
static perfect fluid spheres 

Damien Martin and Matt Visser 

Perfect fluid spheres, either Newtonian or relativistic, are the flrst step in 
developing reahstic steUar models (or models for fluid planets) . Despite the 
importance of these models, explicit and fully general solutions of the perfect 
fluid constraint in general relativity have only very recently been developed. 
In this Brief Report we present a variant of Lake's algorithm wherein: (1) 
we re-cast the algorithm in terms of variables with a clear physical meaning 
— the average density and the locally measured acceleration due to gravity, 
(2) we present explicit and fully general formulae for the mass proflle and 
pressure proflle, and (3) we present an explicit closed-form expression for 
the central pressure. Furthermore we can then use the formalism to easily 
understand the pattern of inter-relationships among many of the previously 
known exact solutions, and generate several new exact solutions. 
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Appendix B 



Bounds on the interior 
geometry and pressure profile 
of static fluid spheres 

Damien Martin and Matt Visser 

It is a famous result of relativistic stellar structure that (under mild technical conditions) 
a static fluid sphere satisfies the Buchdahl-Bondi bound 2M/R < 8/9; the surprise here 
being that the boimd is not 2M/R < 1. In this article we provide further generalizations 
of this bound by placing a number of constraints on the interior geometry (the metric 
components), on the local acceleration due to gravity, on various combinations of the 
internal density and pressure profiles, and on the internal compactness 2m{r)/r of static 
fluid spheres. We do this by adapting the standard tool of comparing the generic fluid 
sphere with a Schwarzschild interior geometry of the same mass and radius — in particular 
wc obtain several results for the pressure profile (not merely the central pressure) that are 
considerably more subtle than might first be expected. 
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Appendix 121 Bounds on perfect fluid spheres 



Appendix C 



Dirty black holes: Quasinormal 
modes 

A.J.M. Medved, Damien Martin and Matt Visser 

In this paper, we investigate the asymptotic nature of the quasinormal modes for "dirty" 
black holes — generic static and spherically symmetric spacetimes for which a central 
black hole is surrounded by arbitrary "matter" fields. We demonstrate that, to the leading 
asymptotic order, the [imaginary] spacing between modes is precisely equal to the surface 
gravity, independent of the specifics of the black hole system. 

Our analytical method is based on locating the complex poles in the first Born approx- 
imation for the scattering amplitude. We first verify that our formalism agrees, asymp- 
totically, with previous studies on the Schwarzschild black hole. The analysis is then 
generalized to more exotic black hole geometries. We also extend considerations to space- 
times with two horizons and briefly discuss the degenerate- horizon scenario. 
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AppendixlH; Dirty black holes: QNM 



Appendix D 



Dirty black holes: Quasinormal 
modes for "squeezed" horizons 

A.J.M. Medved, Damien Martin and Matt Visser 

We consider the quasinormal modes for a class of black hole spacetimes that, informally 
speaking, contain a closely "squeezed" pair of horizons. (This scenario, where the relevant 
observer is presumed to bo "trapped" between the horizons, is operationally distinct from 
near-extremal black holes with an external observer.) It is shown, by analytical means, 
that the spacing of the quasinormal frequencies equals the surface gravity at the squeezed 
horizons. Moreover, we can calculate the real part of these frequencies provided that the 
horizons are sufficiently close together (but not necessarily degenerate or even "nearly 
degenerate"). The novelty of our analysis (which extends a model-specific treatment by 
Cardoso and Lemos) is that we consider "dirty" black holes; that is, the observable portion 
of the (static and spherically symmetric) spacetime is allowed to contain an arbitrary 
distribution of matter. 
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Appendix\D[ DBH: QNM for "squeezed" horizons 



Appendix E 



A note on quasinormal modes: 
A tale of two treatments 

A.J.M. Medved and Damien Martin 

There is an apparent discrepancy in the hterature with regard to the quasinormal mode 
frequencies of Schwarzschild-de Sitter black holes in the degenerate- horizon limit. On the 
one hand, a Poschl- Teller-inspired method predicts that the real part of the frequencies 
will depend strongly on the orbital angular momentum of the perturbation field whereas, 
on the other hand, the degenerate limit of a monodromy-based calculation suggests there 
should be no such dependence (at least, for the highly damped modes). In the current 
paper, we provide a possible resolution by critically re-assessing the limiting procedure 
used in the monodromy analysis. 
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Appendix\^ A tale of two treatments 



Appendix F 



Dirty black holes: Spacetime 
geometry and near-horizon 

symmetries 

A J M Medved, Damien Martin and Matt Visser 

We consider the spacetime geometry of a static but otherwise generic black hole (that 
is, the horizon geometry and topology are not necessarily spherically symmetric). It 
is demonstrated, by purely geometrical techniques, that the curvature tensors, and the 
Einstein tensor in particular, exhibit a very high degree of symmetry as the horizon is 
approached. Consequently, the stress-energy tensor will be highly constrained near any 
static Killing horizon. More specifically, it is shown that — at the horizon — the stress- 
energy tensor block-diagonalizes into "transverse" and "parallel" blocks, the transverse 
components of this tensor are proportional to the transverse metric, and these properties 
remain invariant under static conformal deformations. Moreover, we speculate that this 
geometric symmetry underlies Carlip's notion of an asymptotic near-horizon conformal 
symmetry controlling the entropy of a black hole. 
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Appendix\F\ DBH: Spacetime geometry and near horizon . 



Appendix G 



Dirty black holes: Symmetries 
at stationary non-static 
horizons 

A J M Medved, Damien Martin and Matt Visser 

We establish that the Einstein tensor takes on a highly symmetric form near the Killing 
horizon of any stationary but non-static (and non-cxtremal) black hole spacetimc. [This 
follows up on a recent article by the current authors, gr-qc/0402069, which considered 
static black holes.] Specifically, at any such Killing horizon — irrespective of the horizon 
geometry — the Einstein tensor block-diagonalizes into "transverse" and "parallel" blocks, 
and its transverse components are proportional to the transverse metric. Our findings are 
supported by two independent procedures; one based on the regularity of the on-horizon 
geometry and another that directly utilizes the elegant nature of a bifurcate Killing horizon. 
It is then argued that geometrical symmetries will severely constrain the matter near any 
Killing horizon. We also speculate on how this may be relevant to certain calculations of 
the black hole entropy. 
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Appendix\^ DBH: Symmetries at stationary horizons 



